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ABSTRACT 


Critical phenomena in the Arnold Engineering 
Development Center (AEDC) High Performance Demonstration 
Experiment (HPDE) and the U.S. U-25 Experiment, are analyzed. 
Also analyzed are the performance of a NASA-specified 500 
UW(th) flow train and computations carried out by STD Research 
under Contract AC-01-79ET15501 concerning critical issues for 
the scale-up of UHD Generators. 

The HPDE is characterized by computational simulations 
of both the nominal conditions and the conditions during the 
experimental runs. The steady-state performance is discussed 
along with the Hall voltage overshoots during the start-up and 
shutdown transients. The results of simulations of the HPDE 
runs with codes from the Q3D and TRANSIENT code families are 
compared to the experimental results. The results of the 
simulations are in good agreement with the experimental data. 

Additional critical phenomena analyzed in the 
AEDC/HPDE are the optimal load schedules, parametric variations 
around the simulations of HPDE Run 006-014, the parametric 
dependence of the electrode voltage drops, the boundary layer 
behavior, near electrode phenomena with finite electrode 
segmentation, and current distribution in the end regions. 

The U.S. U-25 experiment is characterized by 
computational simulations of the nominal operating conditions. 
The steady-state performance for the nominal design of the U.S. 
U-25 experiment is analyzed, as is the dependence of 
performance on the mass flow rate. 
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A NASA-specified 500 MW(th) MHD flow train is 
characterized for computer simulation and the electrical, 
transport, and thermodynamic properties at the inlet plane are 
analyzed. 


Issues for the scale-up of MHD power trains are 
discussed. The AEDC/HPDE performance to date is analyzed to 
compare these experimental results to scale-up rules. The 
optimum Mach number distribution is analyzed with emphasis upon 
its effect on part- load and transient behavior. The effects of 
alternate cross-sectional shapes on channel performance are 
also evaluated. 
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ANALYTICAL INVESTIGATION OF CRITICAL PHENOMENA 
IN HHD POWER GENI^ATORS 


1.0 INTRODUCTION 


U.S. MHO technology for commercial power generation 
has made a significant transition in 1980. Test data is being 
obtained for MHD power trains of sufficient scale, interaction, 
and durability to bridge the gap between prior laboratory scale 
cperiments and the requirements of coamiercial systems. In 

addition, detailed design studies are being initiated for 
advanced power train concepts which will ultimately lead to 
pilot-scale tests of a coal-fired MHD power plant. The 

significant power train projects in this transitional phase of 
MHO technology include the AEDC High Performance Demonstration 
Experiment (HPDE), the U.S. U-25 generator/diffuser; U.S. 
participation in the U-25B experiments, the UTSI Coal Fired 
Flow Facility (CFFF), the Advanced Power Train (APT) for the 
Componert Development and Integration Facility (COIF), and the 
Engineering Test Facility (ETF). 

These projects each represent an important advance 
beyond previous laboratory experience. In contrast to past 
laboratory experiments, these projects test MHD power trains in 
which the interaction between the electrical forces and working 
fluid will be strong enough to cause important departures from 
classical internal duct gasdynamics. There are at present only 
limited experimental data to guide the prediction of the 
performance and behavior of MHD channels in such moderate-to- 
high interaction situations. Likewise, there is a lack of 
experimental data against which the test results of these 
programs might be compared. Under NASA Contract DEN3-179, STD 
Research Corporatior has continued the application of its 



comprehensive analytical capabilities to predict and interpret 
the behavior of high interaction MHD power trains. This Final 
Report for Contract DEN3-179 describes the results of analyses 
of high interaction experiments and designs carried out by STD 
during this transitional phase of MHD technology. 

The planning, detailed analysis, and interpretation of 
the ABDC/HPDB and the U-25 experiments have been the subject of 
concern at STD Research Corporation from their initial 
conception [1-1]. During many years of contributions to the 
initiation and evaluation of these experiments, STD Research 
Corporation has analyzed many of the interim designs for the 
experimental hardware. Under U.S. Department of the Interior, 
and U.S. Department of Energy contracts since 1971, STD 
Research has made numerous specific recommendations for 
improved test parameters and procedures for the HPDE and U.S. 
U-25 tests. Many of the STD recommendations have resulted in 
changes in test hardware and/or procedures. 

Under DOE contract AC-01-79ET15501 and NASA contract 
DEN3-179, STD Research has analyzed the "as-built" 
specifications of the HPDE and U-25 experiments under nominal 
operating conditions. These analyses have quantitatively 
predicted the effects of critical, nonideal processes which 
will control the performance of the HPDE and U-25 tests. 
Moreover, the work undei these contracts has mapped the 
predicted nominal operating characteristics of those 
experiments over the likely ranges of mass flow rate, working 
fluid composition, diffuser recovery coefficient, wall 
temperature, and other operating conditions. 

Under Contract DEN3-179, STD Research has carried out 
analysis of preliminary data from the first tests of the HPDE. 
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These experimental data confirm earlier analytical predictions 
by STD Research Corporation, such as the presence of large 
axial voltages during the start-up transient and the presence 
of much- larger- than design electrode voltage drops. This 
experimental confirmation of effects predicted by the STD 
analytical tools further validates the accuracy of these tools. 

Studies of the nominal operation of the HPDE have 
yielded estimates of the relative importance of critical 
phenomena expected in the flow train tests. Unfortunately, the 
experiments can never run at precisely the nominal conditions. 
A significant test of an analytical tool is its capability to 
utilize Inputs from instantaneous measurements taken during a 
test, and then to predict other quantities measured at the same 
instant during the same test. Such data are now available from 
the HPDE, and comparisons between such data and the STD/liHD 
code results are presented in this report. 

Previous test data analyses at STD Research 
Corporation [1-2] have amply demonstrated the benefits of 
careful test data analysis for providing data consistency and 
understanding of the relative importance of t>^e operative 
physical processes. As data have been received fr<»n the HPDE, 
STD Research Corporation has undertaken a detailed analysis of 
appropriate, self-consistent data sets which have led to (1) 
better understanding of the significance of each test .nd (2) 
recommendations for Improved test procedures. 

The survey of the nominal performance of the U-25 
generator designed and constructed in the United States was 
carried to the point of evaluati.jn of the "as-built" 
performance of the ebanne) . The channel has not been operated 
to date, and therefore the analytical predictions ol the 
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studies of lUMBinal performance have not been checked against 
experimental data. Insights gained by the comparison between 
the analytical simulation of AEDC/HPDE experimental data were 
applied to the U-25 nominal operating conditions, however. 
These analyses focused on the implications of operating the 
U-25 channel with cold walls without modification of the 
original design loft. The results indicated that the electrode 
boundary layer voltage drops would assume a major role in the 
achievement of the 10 liW design goal with the “as-built” U-25 
nominal operating conditions. 

Perhaps the most important objective of analytical 
modeling of any physical process is the development of the 
ability to extend the results of experiments under one set of 
conditions to the conditions of other experiments. To do this, 
one must have the ability to exercise sufficiently rigorous and 
complete physical models which do not depend upon the 
application of adjustable, ^pirical factors in order to obtain 
agreement with experimental data. The use of such rigorous 
analytical tools, as are employed at STD Research Corporation, 
place a substantial burden on the experimentalist to provide 
precise, complete measurements for all of the physical data 
required as input to the calculation. Experience at STD 
Research Corporation [1-4] has shown that the combination of 
well diagnosed experimental data and rigorous analyses can pay 
dividends in understanding complex physical mechanisms and the 
experimental data itself. After the initial “shake down" 
phases of HPDE testing, the HPDE is beginning to produce such 
data. The simulations performed under Contract Di5N3~179 have 
been fruitful as a consequence. 

It is of interest to extrapolate the results of the 
AEDC/HPDE experimental data analyses to-date to conditions of 
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higher interaction. These conditions Include operation of 
pilot scale and commercial scale MHD generators. Under 
Contract 0EN3-179 and its predecessors, STD Kese;«rch 
Corporation has characterized the behavior of MHO power trains 
from laboratory scale to commercial scale utilizing a self- 
consistent set of analytical tools. It was shown in [1-3] that 
the performance and fluid behavior of MHD power trains 
correlates well with the appropriate interaction pare''-^tc . 
The performance and interaction parameters of the AEDC/HPDE 
relative to generators of other scaies are survey'*'* in the 
present report . 

Utilizing the results of calculations carried out 
under Contract AC-01-79ET15501 , the scale dependence of the 
performance of generators of alternate cross-sections 
(elliptical and hexagonal) also has been analyzed. The 
benefits of alternate cross-sectional geometries are strongly 
dependent on the boundary layer characteristics present at a 
particular station in a given channel. The powei generation 
potential of each of '^hree cross-sections in each of four 
channels ranging in scale from 20 MW thermal input to 2000 MW 
thermal input is compared. The comparison has been carried out 
between the power actually generated and (1) the power 
generation potential predicted by simplified electrical 
calculations, (2) power generation with rectangular geometry, 
and (3) ideal power generating capability based on center-line 
conditions at each axial station. While the analyses are 
indicative of the potential benefits of alternate cross- 
scctional configurations in MHD generators at various scales, 
the results do not represent optimum configurations. A number 
of additional computational degrees of freedom which were not 
considered in the previous study may also serve to improve the 
performance of generators of alternate cross-sections beyond 
the levels demonstrated. 


5 



Finally, an analytical study of the scale dependence 
of part load and transient behavior of MHD power generators is 
summarized in the present report. It is shown that a critical 
parameter for the successful part load operation of MHD power 
generation is the velocity or Mach number distribution in the 
MHD generator. The conclusion of this study has been that the 
key to successful off-design or part load operation of MHD 
generators is the maintenance of the Mach number distribution 
at or near the optimum value for the stagnation conditions 
existing at each station in the MHD generator. 

This Final Report is organized according to the major 
topics addressed during the course of Contract DEN3-179. 
Section 2 presents a brief description of the codes utilized 
durii g the analytical Investigations, and points out documents 
and literature in which fuller descriptions of the codes may be 
found. Section 3 focuses on the analysis of the "as-built" 
AEDC/HPDE, including the analysis of actual test data. Section 
4 is devoted to special -Studies based upon the simulations of 
AEDC/HPDE experimental data described in Section 3. Section 5 
describes the characteristics of the nominal, "as-built", 
operating conditions and specifications of the U.S. U-25 flow 
train. Section 6 describes the results of thermodynamic 
calculations of coal combustion products, including estimates 
of the electrical conductivity at the inlet of a NASA specified 
500 MWth MHD generator. The purpose of Section 7 is to put the 
results of the present contract and Contract AC-01-79ET15501 
into the context of the overall development program for MHD 
power generation. Sections 8 and 9 summarize the findings of 
the report and enumerate the major recommendations resulting 
from the studies presented. 
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2.0 DESCRIPTION OF THE STD/MHD CODES 

The STD/MHD codes are constructed from a large 
collection of modules or subprograms which address various 
aspects of the MHD problem. Taken together, these codes define 
and solve the MHD problem, including, where applicable: 


1. Viscosity and wall roughness 

2. Volume viscous losses 

3. Wall roughness effects on skin friction 

4. Other skin friction effects 

5. Compressibility 

6. Turbulence, including MHD effects 

7. Turbulence structure parameters and transport coefficients 

8. Convective heat transfer 

9. Radiative heat transfer 

10. Anisotropic, nonuniform, fluctuating electrical transport 

properties 

11. Nonuniform, fluctuating plasma/fluid properties 

12. Nonequilibrium effects 

13. Finite reaction rates, including ionization/recombination 

14. Electron energy relaxation 

15. Electron radiation losses 

16. Electron thermal and concentration diffusion 

17. Nonlinearity of the plasma as a conducting medium 

18. Electrode surface and sheath phenomena 

19. Electrochemical effects at the electrodes and slag layer 

20. Finite electrode effects 

21. Flow separation 

22. Current leakage due to imperfect segmentation or slag 

layer leakages 

23. Axial current leakage in the plasma and interelectrode 

arcing 

24. Trans-boundary layer arcing 

25. Shock waves and shock wave/ boundary layer interactions 

26. MHD instabilities Cincluding magnetoaerothermal) 

27. Ionization instability modes 

28. Current streamer formation 

29. Start-up transients, including sequencing of fuel, 

oxidizer, and seed valves, secondary Injection of fuel, 
seed and oxidizer, and load circuit switching 

30. Shut-down transients 

31. Perturbations from steady state operation, including steam 

plant upsets, load circuit faults, part-load operation 
and load following transients 

32. Generator/power conditioning system interactions 

33. Generator end region losses 

34. Slag layer dynamic phenomena 
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35. Channel geometry effects. Including nonrectangular 

cross-sectional shapes 

36. Transverse pressure gradients and nonuniformities 

37. Secondary flows induced by MUD or viscous effects 

38. Fundamental combustion processes, including particle 

burnout, particle trajectories, liquid layer formation 
on combustor walls, seed mixing and uniformity, pressure 
loss, etc. 

39. Aliv^n waves 

40. Hartmann flow 

41. Nucleatlon effects 

42. Electron capture by particles 

43. Radiative ionization and recombination 

44. Nonequilibrium radiation and coherent radiation 

45. Effects of nonideal mixtures in condensed phase chemistry 

46. "Faraday Catastrophe" load faults 

47. Three-dimensional and time-dependent Induced magnetic 

field effects 

48. Ion slip effects 

49. Nonideal plasma effects 

50. Combustor-generated inlet swirl 

The modules or subprograms describing these and other 
processes can be coupled in a variety of ways, depending upon 
the desired focus of a given calculation. 


Detailed descriptions of the formulations, solutions, 
and applications of the appropriate STD/HHD codes are contained 
in Ref. (2-1] through 1 2-5). A capsule description of each of 
the TRANSIENT, Q3D, PROPERTIES, FIN, INLET, BONDLAY, and 
GRAPHICS code families follows. 


Code Family : TRANSIENT 

Time-dependent magnetohydrodynamics, plasmadynamics, and 
compressible gasdynamics for internal flows with MHD 
power extraction and body forces. Comprehensive code 
family accounts for two-dimensional effects, wall effects 
(heat transfer and friction) and kinetics. Applicable to 
MHD power systems, gas-dynamic lasers, and other high 
temperature flow systems in which start-up, shutdown or 
unsteady operation are important. TRANDEE is one of the 
specialized codes within this family. The User's Guide to 
this code is available, Ref. (2-3). 


8 



Code Family : Q3P 


Quasi-three-dimensional , steady HHD, electrical, and fluid 
mechanical effects in equilibrium (combustion-driven) or 
nonequilibrium (noble gas) IIHO flow trains (combustor, 
nozzle, channel, and diffuser) of arbitrary cross section. 
This code family includes various cross sectional MHD 
generator performance calculations which may be coupled to 
predict three dimensional performance and behavior. While 
finite segmentation options exist, they were not utilized 
in the performance of the Q3D calculations described in 
this report. QUE30EE is one of the specialized codes 
within this family. A User's Guide is available for this 
code. Ref. (2-4J. 

Code Family : PROPERTIES 

Generalized transport and thermodynamic properties of 
equilibrium combustion flames. These codes consider 
reactions between more than 500 gas, liquid, or solid 
state species. Flexible inputs include initial fuel, 
oxidizer, and additive (seed, ash, etc.) compositions and 
states. Outputs include electrical and fluid transport 
coefficients and all thermodynamic state variables. 

Code Family : FIN 

Two-dimension?* 1 MHD electrode solutions with electron 
energy relaxation, electron ionization/recombination 
kinetics, effects of slag coatings, and finite electrode/ 
insulator width. Codes in this family apply to the region 
of the MHO generator channel which may be considered 
"periodic"; i.e., regions in which changes from electrode 
pair to adjacent electrode pair are essentially 
negligible. 

Code Family : BONPLAY 

Integral boundary layer models with MHD body forces. Joule 
dissipation, wall roughness effects. 



Code Family : GRAPHICS 

General and customized graphics packages for 
post-processing raw data files generated by other code 
families. Includes codes for contour generation and three 
dimensional projection of bivariate data. 

Code Family : INLET 

Codes in this family solve the two- dimensional electrical 
part of the MHD problem in the part of the MHD generator 
which cannot be considered periodic. These codes take 
given gasdynamic variables and compute the two dimensional 
distributions of electron temperature, Ohm's law coeffi- 
cients, current density, and electric field in the end 
regions of an MHD channel. 
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3.0 


HPDE DATA ANALYSIS AND INTERPRETATION 


STD Research Corporation has performed calculations to 
survey the power producing characteristics of the AEDC/HPDE 
over the nominal operating range as well as calculations to 
simulate actual experiments which have occurred to date. The 
results of these studies are discussed in this section. 

3 . 1 Channel Characterization 


Subsection 3.1.1 describes the physical parameters 
used for simulations of the nominal operation of the AEDC/HPDE. 
Descriptions of operating parameters which differed from these 
nominal conditions during actual tests are provided in 
subsequent subsections. 


3.1.1 AEDC/HPDE - Nominal Conditions 

The flow train of the Arnold Engineering Development 
Center High Performance Demonstration Experiment (AEDC/HPDE) 
consists of a combustor, transition nozzle, electrode test 
section, generator, generator/diffuser adapter section and 
diffuser. Fig. 3-1 is a schematic of the final, "as-built" 
internal dimensions of the device (axial locations are 
referenced from the combustor back plate). 

3. 1.1.1 Geometry 

The combustor for the HPDE is based on renovation and 
modification of what was formerly the LORHO burner and is 
described in 13-1] through [3-4]. For the purpose of 
performance analysis, the combustor may be char9cterized as a 
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0.610 m diameter cylinder of 1.195 m length. The transition 
nozzle is a three-dimensional contraction from the 0.605 m 
diameter burner to a 0.264 m by 0.489 m rectangular electrode 
test section. The electrode test section is of no consequence 
during the MHD experiments and may be regarded simply as an 
unloaded portion of the generator. The MHD generator channel 
is constructed in 5 sections, the dimensions of which are shown 
in Fig. 3-1. In the STD computer codes the x-axis is located 
at the center of the chancel cross-section. The y-axis is in 
the electrode-to-electrode direction with the positive 
direction being from cathode toward anode. The z-axis is in 
the sldewall-to-sidewall direction. The magnetic field vector 
is in the positive z direction. The channel height, width, 
area, and aspect ratio are plotted as a function of the axial 
coordinate in Fig. 3-2. 

There are 485 pairs of electrically accessible 
electrodes, of which only 417 pairs spanning 7.15 m in the 
Interior of the channel are loaded under the nominal operating 
load Schedule. 

The electrodes are electrically isolated from each 
otutr with refractory, and each electroae spans the transverse 
distance across the channel. The gap between electrodes is 
1.52 m throughout the channel [3-5]. Graphite caps are 
attached to each electrode pair to allow operation at a high 
surface temperature. 

The "pegwall” insulator walls are composed of 0.019 
square conducting pegs which are spaced 1.58 mm apart and are 
insulated with refractory material [3-4). 

The generator is fitted (via an adapter section) with 
a plane-walled diffuser with parallel sidewalls (extensions of 
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the electrode walls) and diverging top and bottom walls. The 
mechanical design 13-1) of the diffuser allows for setting the 
upper and lo*«er walls at any angle between 0 and 3®. Experi- 
ments to date have fixed the divergence angle at 1° on both 
walls. 

3. 1.1. 2 Magnetic Field 

The magnet for the MHD High Performance Demonstration 
Experiment is designed to provide (at 27 MW and 16 kA) a peak 
field of 6 T when ^/ulsed from 77 K (IJN2) and 3.5 T when 
operated continuously water-cooled. The magnet bore is 0.89 ra 
wide by 0.71 m high at the entrance, 1.4 m wide by i.l7 m high 
at the exit, ?nd the poles are 7.1 m long. The total length of 
the magnet is a 9.16 m. 13-6) 

Fig. 3-3 shows the agreement between the predicted and 
measured magnetic field distributions. STD simulav\ons of the 
AEDC/HPDE experiments to date have used, as inp.. t for the field 
distribution, the MEA design values (3j-i'jj. These values are 
scaled by the measurement of Hall probe H9 for each particular 
run under analysis. 

3. 1.1. 3 Loading 

The nominal loading scheme for the HPDE generator is 
in a Faraday configuration. For purposes of performance survey 
calculations, the load distribution is given as follows (3-5). 

x(m) RL (ohm-m) 

2.735 < X < 8.875 0.8 

8.875 < X ± 9.885 1.6 

where x is measured from the burner backplate. The electrode 
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spacing is given in [3-7] and is included in Table 3-1 for 
convenience. 


3. 1.1. 4 Working Fluid 

The working fluid for the channel is composed of the 
products of the combustion of toluene and oxygen with nitrogen 
Qiluent. The seed is potassium hydroxide dissolved in 
methanol. The fuel/seed irixture is injected into the com' .stor 
as a liquid at ambient temperature. The nominal mass flow 
rates for the facility are 45 kg/s for operation between 2 and 
3 T, 50 kg/s for opi ration between 3 and 4 T, and 53 to 54 kg/s 
for operation at 4 T and greater (3-5]. The nominal value for 
both the N 2 /O 2 mole ratio and the stoichiometry is 1.0. Table 
3-2 contains a list of the flow rates which describe the 
working fluid under nominal conditions (2 T, 45 kg/s). 

3. 1.1. 5 Wall Conditions 

The nominal values for the surface temperatures on 
both the insulating and electrode walls are shown in Fig. 3-4. 
These temperatures were determined by a one-dimensional 
ransient heating analysis performed by AEOC personnel and 
described in (3-7J. For all calculations the duration of 
heating was taken to be 15 s, after which there was a cool-down 
cycle during w» ich atmospheric air is forced through the 
channel. The temperature history for each wall is calculated 
at the channel inlet and exit, and the distribution of the wall 
temperatures along the channel were obtained by assiunlng a 
linear variation between these locations (for values corre- 
sponding to t « 15 s). The combustor operates with a wall 
temperature of 450 K, and the rorghness of all walls is 
charactr>rized by an equivalent sand roughness of 3 mm in the 
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electrode test section and channel, and as "smooth" in the 
combustor, nozzle, and diffuser [3-5]. 

3. 1.1. 6 Channel Inlet Conditions 

The character of the flow as it enters the UHD 
generator is not completely defined by the experimental data. 
Ideally, one would prefer to know with confidence the complete 
distribution of mass, momentum, and energy at the inlet plane. 
Available data are sparse due to the experimental difficulty of 
obtaining these measurements. The heat loss to the water- 
cooled components of the burner is on the order of 3% (3-4). 
The total heat loss to all components upstream of the first 
loaded electrode was calculated, as described in Section 
3. 1.2.5, by STD Research; the value of 5.7% was used for the 
nominal case. The original design calculations based on the 
nominal operating conditions indicated that the boundary layer 
thickness is approximately 12 mm at the inlet plane. 

3. 1.1. 7 Other Conditions and Assumptions 

The generator exit condition is determined by the fact 
thu.t the HPDE flow train exhausts to the atmosphere. Hence, 
the static pressure in the exit plane of the diffuser is 
assumed to have the nominal value of one atmosphere. 

3.1.2 AEDC/HPDE - Experimental Conditions 

The HPDE experimental conditions are chosen according 
to experience learned from each successive test and a test plan 
which includes variation of the key test parameters. Conse- 
quently, the HPDE runs have not, to date, operated at exactly 
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the "nomiaal" conditions pr .ented in tbe previous section. 
Post-experimental data analysis is the only way to determine 
the operating point of a particular experiment. It is the 
purpose of this section to set foi ch as completely as possible, 
using tbe data at hand, the experimental conditions under which 
the AEDC/HPDE has actually operated to this date. 

3.1.2. 1 Magnetic Field 

During power producing runs the magnetic fxeld 
strength is monitored via Hall probes, tbe locations of which 
are shown in Fig. 3-5. All probes are located on the floor of 
the magnet bore. A measure of the variation of the magnetic 
field in the interelectrode direction is provided by probes H2A 
and U3A which are displaced above and below the mid line (y ^ 0) 
by 0.014 m. The location of these probes with respect to the 
design magnetic field distribution is shown in Fig. 3-3. Note 
that they are near the peak magnetic field region. For the 
purpose of computer simulation of the experiments, the measured 
values of probes H2A and H3A are averaged over the time 
interval of interest. The resulting value of the magnetic 
field is then used to scale the design curve to yield the 
complete magnetic field distribution ior tbe experiment under 
consideration. Table 3-3 contains the values of the magnetic 
field so obtained for the AEDC experiments to date. 

3. 1.2. 2 Loading 

The load bank is a series of liquid rheostats 
consisting of polypropylene buckets, each containing two copper 
plates. The plate spacing can be varied to achieve a 
resistance range from 10 to 60 ohms. Table 3-4 contains the 
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resistance of each electrooe pair as set by AEUC personnel 
13-8|, (3-9|. These values have been normalized by a water 
conductivity of 150 pS/cm. Note that the reference resistances 
for electrode pairs numbered greater than 357 have been changed 
for runs subsequent to 006-007 and also that electrode pairs 
401 and above are open circuited. Table 3-5 contains the water 
conductivity during the runs for all runs to date. To deter- 
mine the resistance schedule for any particular run, the values 
from Tables 3-4 and 3-5 must be used in the following equation: 

“ruB " '•rat 

The product of t.he above load resistance and the electrode 
pitch (width from Table 3-1 plus the insulator thickness of 
1.52 mm) then uniquely specifies the loading of the channel in 
the Faraday configuration. It should be noted here that 
comparisons of the load resistances determined in the above 
manner and those determined from the vol tag^»-curren t 
characteristics of the experimental data are not very good for 
Run 006-008, but were better for subsequent runs. Further 
discussion on this point Is contained in Sectit.is 3.6 and 3.7. 

3. 1.2. 3 Working Fluid 

Table 3-6 contains a list of the flow rates of the 
various constituents of the working fluid for two typical 
experimental runs (006-008 and 006-014). These flow rate data 
are input to codes from the STD THERMODYNAMICS family to 
generate a thermodynamic data base and other data. Input flow 
rates for oxidizer, seed, and fuel were obtained by averaging 
the experimental data over the time period 4.54 £ Tg < 5.21 s 
for Run 006-014 and 9.5 ^ ^2 S. s for Run 006-008. 
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3. 1.2.4 Wall Conditions 


Fig. 3-6 illustrates the results of a transient 

heating analysis performed by AEDC personnel for the electrode 

surface including the 1/2" graphite caps (3-11 J. Experimental 

data have indicated that the heat transfer 

rates on the electrode surfaces at the inlet flange and at a 

2 

station 5.59 m from the inlet flange are 450 Etu/ft -sec and 
2 

100 Btu/ft -sec, respectively. Several of the power producing 
runs to date have yielded their best data about 5 seconds into 
the run. Reading the curves of Fig. 3-6 (at t = 5 s) for the 
heat transfer rates mentioned above yields electrode surface 
temperatures of 1130 K and 540 K at the inlet flange and a 
point 5.59 m downstream, respectively. For the purpose of 
performance calculations, the complete temperature distribution 
is assumed to be linear along the channel between these two 
points. These data were the best available at the time these 
simulations were performed. 

All simulations beginning with the simulations of Run 
006-014 utilized the temperature profile described in the above 
paragraph. Simulations of runs prior to 006-014, which were 
carried out before this information became available, utilized 
the nominal temperature distribution. Fig. 3-4. 

As the HPDE runs were analyzed, better wall 
temperature data became available. A further revision to the 
electrode surface temperatures has been published in (3-12] and 
is reproduced in Fig. 3.7. This estimate is based on a two- 
dimensional, fully viscous calculation including the effects of 
pressure gradient. 
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3. 1.2. 5 CbanDel Inlet Conditions 


To determine the boundary layer thickness at the first 
loaded electrode, a series of computations have been made with 
a code from the STD BONDLAY Family. Several runs were made 
assuming various values for the equivalent sand roughness in 
the components upstream of the channel w? ch the following 
results at the first loaded electrode: 

Sand Roughness Height 5( electrode wall) 6( sidewall) 


(mm) (mm) (mm) 


0.0 (smooth) 17.6 1.9.3 

0.1 24.0 

1.0 34.9 38.4 

2.0 40.0 44.3 

3.0 43.6 48.4 


These calculations also predict that the total 
enthalpy loss up to the first loaded electrode ranges between 
approximately 5.7% for the smooth wall calculation to 
approximately 7.2% for the 3.0 mm roughness calculation. The 
latter value compares well with the datum of 7.5% measured by 
AEDC (3-12]. 
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density, integrated power output, and interelectrode voltages 
are determined from tne primary measurements and other 
geometrical data. 

3.2.3 STD Interface with HPOE Data Systems 

The HPDE data acquisition system data are processed by 
codes from the STD GRAPHICS family and plotted as a function of 
both space and time. These codes allow the data and 
computational results to be displayed with any degree ot time 
and space resolution necessary, contingent of course, on the 
ultimate density of the data, to examine events of importance. 
Fig. 3-8 is illustrative of the manner in which critical events 
can be followed throughout their lifetime using this type of 
display. 

A plot of the data acquisition system data for a typical 
AEDC run (006-008) is given in Fig. 3-9. On this figure the 
time coordinate increases from lef t-to-right and starts at T 2 = 
0 (Tg is the time in seconds from combustor ignition). Also, 
the axial distance along the flow train is plotted from 
top-to-bottom. On this figure, the value of the data is 
plotted vertically above the corresponding point in the 
time-distance coordinate plane. It is possible to define and 
plot subsets or "windows" in the data as small as necessary to 
resolve an event of interest if the data are available. Such 
options have as yet gone unused for the HPDE data due to the 
relative coarseness of the experimental data in space and time. 
Much higher data rates will be required t"> resolve the space- 
time variation of critical phenomena such as HHD instability 
modes, electrical transients, etc. 
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3 . 3 Steady State Pe rfor mance at the Nom i na l Ope rati ng 
Conditions 

Prior to receiving the experimental data from the 
initial powered runs, a study was conducted in order to 
determine the probable performance range of the AEDC generator 
during the first phase of testing. The nominal operating 
conditions were defined by AEDC personnel and are documented in 
Section 3.1.1. During these computations the peak magnetic 
field was varied between 2.0 and 4.0 T and the mean flow rate 
was varied between 45 and 54 kg/s. Due to the unavailability 
of sophisticated FIN or ARRAY calculations in the early stages 
of te HPDE analysis, a variety of wall condition models were 
employed in Q3D for the initial calulations of the HPDE 
performance under the nominal operating conditions. 

With the arc mode current transfer option, the current 
transverse through the boundary layer flows with negligible 
loss when the local electric field exceedes a critical value - 
Based on previous work, the critical value used for these 
computations is -12 kV/m. For the diffuse current mode option, 
the current transverse through the boundary layer is 
computationally forced to be uniform across the width and 
length of the electrodes. 

Figs. 3-10 and 3-11 present the summary of the predicted 
performance of the HPDE at the nominal operating conditions for 
2 T, 3 T, and 4 T operation. Performance predictions were 
carried out as a function of mass flow rate with the STD 
TRANSIENT family of codes in the quasi-steady mode of 
operation. The predicted values fall within the bounds of 
uncertainty defined by a series of more sophisticated, 
quasi-three-diraensional calculations with the Q3D family of 
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codes. Depending upon the relative importance of seed 
condensation and current transport mechanisms, the predicted 
power output was found to vary between 4.28 to 15.3 MW at 2 T 
and between 25.96 to 39.76 at 4 T. Pig. 3-11 indicates that 
the experiment should not be expected to obtain 15% enthalpy 
extraction at 4 T and the nominal operating conditions. 

Figs. 3-12 to 3-24 are results of a calculation, 

COBQPXRIHU, with a code from the Q3D family for the 2 T nominal 

condition (Point A (2 T) in Pig. 3-10). These are axial 

profiles of certain key variables which are output from the 

calculation . The origin for the axial coordinate in these 

figures is located at the inlet flange (end of electrode test 

section in Fig. 3-1). Note how the sectioning in the channel 

construction (see Fig. 3-1) is directly reflected in the Mach 

number and velocity distribution of Pig. 3-13. The load factor 

decreases from 0.8 to 0.64 over the active channel (Fig. 3-15). 

Fig. 3-16 shows that approximately 15.3 MW are produced at the 

nominal 2 T conditions at an enthalpy extraction ratio of 

0.059. Maximum normal current densities of approximately 0.5 
2 

A/ cm and maximum Hall fields cf approximately 1.8 kV/m ar^ 
present during operation as shown by Pig. 3-17. Maximum 
voltage drops* of approximately 750 V occur in the rear end of 
the generator (Pig. 3-21). 


♦Except as otherwise noted, voltage drops in this report are 
defined as the difference between the Faraday voltage and the 
voltage difference across the channel if the centerline 
transverse voltage gradient were extrapolated to the 
electrodes. 


24 



3.4 


Hall Voltage Overshoots 


The mechanism for the occurrence of Hall field 
overshoots during start-up and shutdown was described in Vol. 
I, pp. 5-12 to 5-16 of the STD Research Corporation report 
FE2243-17 under contract EX-76-C-01-2243 (May 1978). During 
periods in which the channel is operating at design mass flow 
rates, but is unseeded or slightly seeded, high velocities are 
experienced in the channel. In the presence of the magnetic 
field, B, these higher-than-design velocities can lead to much 
higher-than-desi«;n Hall fields: 

= PUB(1-K)/G 

2 

Taking note that p ^ B/p, we find that (1) ^ B , and (2) 
every factor in the above expression (except the nonuniformity 
factor G) changes in the direction to increase the magnitude of 
Ey as the velocity increases. Such unseeded or slightly seeded 
conditions are present during the start-up and shutdown 
sequences because the seed valve is opened last and closed 
first. 


In the calculations described in the forementioned 

report. Hall fields exceeding 20 kV/m locally were observed 

at the nominal 6 T condition during the start-up and shutdown 

events, and Hall fields exceeding 4 kV/m extended over more 

than half the channel. Scaling these fields by the square of 

the ratio of magnetic flux density in Run 006-008 to the value 

2 

used in the previous study (2.35 T/6 T) , we might expect Hall 
fields of the order of 3 kV/m over the entire channel due to 
the increased extent of the supersonic flow region in AEDC/HPDE 
Run 006-008 (see Section 3.6 for a complete analysis of this 
run). Thus, instantaneous total Hall voltages up to 23 kV 
might be expected in the 006-008 start-up and shutdown events. 
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Fig. 3-25 is a plot of the anode- to- ‘ground voltage for 
electrode pair number 336 during AEDC/HPDE Run 006-008 in which 
the diffuser became grounded during the run. Figs. 3-26 and 
3-27 show the anode and cathode voltages at all instrumented 
electrodes (for the Data Acquisition System) over the entire 
duration of the run. It is clear that the overvoltages 
occurred throughout the channel. The reduced data acquisition 
system data, which are presented in all plots of data 
acquisition systems in this report, are averages over five 
acquisitions in order to filter the noise. Therefore, while 
voltages as high as 25 kV are not displayed in these figures, 
the high standard deviations associated with the peaks in the 
Fig. 3-25 data imply higher instantaneous voltages than the 
averaged values depicted in Fig. 3-25. 

The sequence of events that took place during Run 
006-008 are reasonably well understood and are depicted 
schematically in Fig. 3-28. This figure presents the actual 
measured diffuser-to-ground voltage along with a curve which 
represents the ideal generator voltage which would occur . f 
there had not been an external arc from the diffuser to ground. 
Due to an approximate 2 s delay between fuel injection and seed 
injection into the combustor flow, sufficient time is available 
to establish steady, unseeded combustion flow in the HHD 
channel during start-up. If no electrical faults develop, the 
Hall voltage would have reached 25 kV as discussed above. This 
is indicated by the solid trace in Fig. 3-28. Instead, an 
electrical breakdown from diffuser to ground developed as the 
Hall voltage increased in Run 006-008. The exact breakdown 
voltage is not known from the experimental data due to the long 
averaging times of the data acquisition system and panel meter 
data. If sufficient information about the geometrical 
configuration of the region near the electrical breakdown were 
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known, the air breakdown voltage could be calculated. This is 
not a critical parameter because the breakdown problem at this 
location has since been corrected. 

During the period between the electrical breakdown and 
the onset of seed flow to the combustor, the external arc is 
sustained by passage of a modest amount of current through the 
plasma froL> the combustor to the diffuser breakdown location. 
Independent assessments of the electrical conductivity in the 
plasma during the unseeded period yield an estimate of 0.03 S/m 
+ 0.02 S/m. These estimates were arrived at by (1) examination 
of the STD THERUODYNAMICS code family predictions for the 
electrical conductivity of the unseeded combustion flame, and 
(2) an analysis of the Faraday current and voltage measurements 
during the unseeded portion of Run 006-008 using calculated 
values of the gas velocity. 

By assuming that the electrical conductivity is rela- 
tively constant in the generator during the unseeded portion of 
the run, by assuming that the conductivity is relatively 
uniform over the channel cross-section, and by taking note that 
the measured voltage to ground of anode 1 and anode 10 differ 
by 1000 volts (the corresponding difference on the cathode side 
is 1200 volts), it 3 possible to estimate a total current flow 
during the unseeded period of approximately 10 + 5 A. Such a 
"trickle" of current would be sufficient to maintain the 
external arc from the diffuser to ground. 

As the seed first enters the combustor, the bulk conduc- 
tivity of the plasma increases, and the current available to 
the external arc increases substantially. The steady-state 
current in the external arc, corresponding to the time when the 
seed flow has reached its steady-state value, is estimated from 
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Q30 computations to be in excess of 850 A, and possibly as high 
as 1200 A (see Section 3.6). The bur .ing voltage of the 
external arc decreases as the current increases, according to 
the standard V-I characteristic of electrical arcs. During the 
seed-on transient, the voltages to ground of the loaded elec- 
trodes in the generator undergo a voltage transient which is 
indicated in Figs. 3-25 to 3-27. Unlike these floating 
electrode circuits, the diffuser voltage is fixed by the ar 
burning characteristics; and, consequently, the diffuser does 
not experience this "seed-on" transient. 

During the HPDE shutdown sequence, the seed valve is 
shut prior to fuel cut-off. As the seed is purged from the 
channel , the conductivity required to sustain the high current 
discharge gradually diminishes, and the diffuser voltage 
follows the V-I characteristics of the external discharge. 
Were there are no electrical faults present, the Hall voltage 
would again rise to approximately 25 kV. However, the burning 
voltage of the arc rises to some lower value, and then the 
external diffuser-to-grouud discharge is extinguished. When 
the discharge is extinguished, the diffuser-to-ground voltage 
rises but is Intercepted by the voltage trace which would be 
present without external electrical discharges. Tuis trace 
^alls as the fuel supply to the combustor diminishes and 
combustion ceases. 

The practical lesson from all this: over-voltage 

protection for many tens of thousands of volts rather than 
hundreds or thousands of volts should be provided between power 
train and ground. 
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3.5 


Shock Location 


The range of possible mass flow rates and magnetic field 
strengths obtained in the AEDC/HPDE facility allow for 
operation in both the supersonic and subsonic flow regimes. 
The nominal operating conditions for the initial testing phase 
(described in Section 3.1.1) were designed to insure operation 
in the supersonic flow regime. Ideally, the sh< ck-down to 
subsonic flow would occur downstream of the active channel. 
While surveying the nominal operating conditions with a code 
from the TRANSIENT family, particular attention was given to 
following the position of the shock as it ' aried with mass flow 
rate and magnetic field strength. 

Figs. 3-29 and 3-30 are representative of the results of 
this study at the 2 T nominal condition. The origin of the 
axial coordinate for these plots is the burner backplate 
(channel exit is at 9.8E ra). Fig. 3-31 summarizes the 
predicted location of the normal shock during operation at 2 T, 
3 T, and 4 T with the nominal operating conditions and various 
mass flow rates. It is seen from these figures that the shock 
is predicted to lie outside the active portion of the 
generator, except during 4 T operation. The 4 T shock can be 
pushed out of the generator by a small increase of the mass 
flow rate within the facility limits. Also s* vn are the shock 
locations predicted for unseeded operation for the nominal N/0 
ratio and for a combustion mixture, designated N/0 = 1.25, 
which simulates the off-design conditions for the first hot 
flow tesc of the generator. The data from this test show \ a 
pressure disturbance indicative of a shock system at 
approximately 9.1 m. The STD calculation based on the best 
available interpretation of experiment?. 1 conditions predi'- cs a 
normal shock to exist 1.6 m downstream of the pressure anomaly. 
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This is ^'onsistent with the empirically determined observation 
in supersonic flows that the duct stalls 2-3 diameters upstream 
of a normal shock computed by one-dimensional gas dynamics. 

3.0 Analysis of Run 006-008 


During start-up of Run 006-008 an electrical breakdown 
from diffuser to ground developed as the Hall voltage increased. 
Calcul ons were initiated with codes from the Q3D family to 
determine the leakage current from the diffuser to ground for 
this run and the associated power loss due to operating with the 
grounded diffuser. An input data set, which contains magnetic 
field and load resiotance distributions and thermodynamic data 
from the AEDC Data Package [3-9], was constructed for use in 
these computations. These inputs are discussed more fully in 
Section 3.1.2. 

It should be remarked here that this simulation was per- 
formed prior to receiving the revised wall temperature schedules 
of Pigs. 3-6 or 3-7 and as such the calculation was done using 
the nominal electrode wall temperature schedule of Fig. 3-4. In 
addition, the assumption of smooth walls upstream of the genera- 
tor was made according to the nominal conditions. Subsequent 
simulations of the Bur. 006-014 data suggest that the upstream 
conditions are more accurately characterized by rough walls. 

As a first approximation to simulate the effects of the 
grounded diffuser, a series of Q3D computations was performed in 
which a constant value of is assumed to exist at each cross- 
section of the generator. The actual leakage current would be 
determined as that current which yields distributions of normal 
current density, power, Faraday voltage, and Hall potential most 
closely resembling those measured during the experiment. 
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At the time these calculations were performed, the data 
for the electrode potentials during the run appeared to contain 
a number of uncertainties. The data from the meter panel and 
data acquisition systems were each incomplete. Depending on how 
the data were combined, an estimate fcr the Hall potential 
difference across the active channel between 1100-1600 V might 
be obtained. In addition, calculations to check the load 

resistance at several electrode pairs indicated a few 
discrepancies between the data from the two systems. Many of 
the load resistances inferred from the measured currents and 
voltages differ substantially, particularly at the rear end of 
the generator, from the values used in the simulations. These 
input values were provided to STD on the basis of separate 
measurements of load resistance prior to the run and are 
documented in Section 3.1.2. 


According to AEDC persornel (3-13) there were two causes 
of data inconsistencies in Run 006-008: 

(1) The current measurements by the panel 
meters are subject to a large uncertainty 
because currents of the order of 24 A are 
t--ing measured on a 200 A (full-scale) meter, 
resulting in approximately a + 5 A uncertainty 
in the reading. 

(2) The current transductors used in the data 
acquisition system were probably measuring 
such low currents as to be in a nonlinear 
range of tfc ? probe characteristics. 


It is important to note that at the time of the STD Run 
006-008 simulations, the HPDE might be considered to be in a 
■’shake-down" mode in which the data acquisition and reduction 
procedures were beicg perfected. Some of the data utilized in 
this simulation have been confirmed by AEDC to be inconsistent. 
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These inconsistencies have been corrected by ABDC since the time 
of the STD simulations. When the data acquisition system data 
and the meter panel data for Run 006-008 were refined by AEDC, 
better agreement with the load resistances was obtaiaed. 

Two Q3D runs were made in an attempt to simulate this 
experiment: CMSIPOFIKW with I, * 0 and CMSIQBEIKW with » 850 
A. The axial variation of the gasdynamic, plasma, and 
electrical parameters from CMSIPUPIKW are given in Pigs. 3-32 
through 3-43. The axial variation of these same parameters 
from computation CMSIQBEIKW are given in Pigs. 3-44 through 
3-55. 


Comparison of the axial variations of the gasdynamic 
parameters presented in Pigs. 3-32 and 3-44 indicates that the 
values are nearly the same with a slightly higher velocity, 
Mach number, stagnation temperature, and stagnation pressure 
from CMSIQBEIKW than for CSiSlPUPIKW. Comparison of the axial 
variation of the plasma properties in Pigs. 3-33 and 3-45 
indicates that the Hall parameter is larger at the exit for 
CMSIQBEIKW. As seen by comparison of Pigs. 3-35 and 3-47, the 
Paraday load factor and electrical conversion efficiency are 
not identical when there is a net 1^. Por CMSIQBEIKW, the 
Paraday load factor decreases to about 0.4 at the back end of 
the channel. 

Due to the large value of 1^^ and power dissipation in 
the arc from diffuser to ground, the values of the power 
extraction parameters are much smaller for CMSIQBEIKW, Pig. 
3-48, than for CMSIPUPIKW, Pig. 3-36. The axial variations of 
the electrical field variables in Pigs. 3-37 and 3-49 indicate 
that the net 1^^ causes the Hall field to be negative at the 
front and back ends of the channel. In addition, the Paraday 
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voltage and normal current density are lower with a net 
Comparison of the axial variation of the electrode potential 
from Figs. 3-38 and 3-SO shows the effects of a ne- I_: the 

A 

potentials decrease at the front end of the channel (due to 
negative E ); potentials are smaller throughout the channel; 

sL 

and the differences from cathode to anode (Faraday voltage) are 
smaller throughout the channel. 

Comparison of the electrode voltage drops in Figs. 3-39 
and 3-51 indicates that the voltage drops are the same at the 
entrance, but peak at about a 20% lower value for QIIQBEIKW. 
Comparison of the electrode boundary layer parameters on Figs. 
3-41 and 3-53 indicate that the net causes an increase in 
the anode shape factor and a decrease in the cathode shape 
factor. Comparison of the sidewall boundary layer parameters 
on Figs. 3-42 and 3-54 indicates that the sidewall boundary, 
displacement, momentum, and enthalpy thicknesses are all 
greater with a net 1^^. 

Computation CMSIPUFIKW simulated the conditions which 
might have been obtained in the HPDE if the external short from 
diffuser to ground had not occurred. Computation CMSIQBEIKW is 
a first approximation to simulate the conditions of Run 006-008 
at the time corresponding to T 2 of approximately 9.0 s. This 
computation is only a first approximation because the actual 
caused by the shorted diffuser is not known. Therefore, 
precise correspondence between calculated and measured values 
cannot be expected. 

Fig. 3-56 illustrates the static pressure distribution, 
both measured and predicted, in the active generator. The 
excellent agreement of theory and experiment indicates that 
there were no large-scale anomalous fluid mechanical phcncmena 
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present in the channel during this portion of the run, and 
suggest that the differences between the experimental data and 
the computations described below may be mainly electrical in 
nature. 


The measured and predicted potential of the cathode are 
depicted in Fig. 3-57. The experimental distribution from the 
AEDC meter panel data system was obtained from an AEDC plot, 
which is given in Fig. 3-58. It should be remarked here that 
the meter data for the potential is unreliable over the first 
10 instrumented electrodes, since negative potentials could not 
be accurately measured at the time of the run. Also plotted in 
F:g. 3-57 is the cathode potential recorded by the data 
acquisition system. Both experimental distributions have been 
plotted with an offset from the recorded values such that the 
potential at electrode no. 168 is zero. 

The Q3D calculations as originally planned would 
determine the constant value of in the generator such that 
the Hall potential difference computed from entrance to exit 
compared well with the measured values. The unreliability of 
the meter data at the inlet and the sparsity of the data from 
the acquisition system at the exit prevents a precise 
determination of the potential difference. Simultaneous 
consideration of both sets of data indicate a Hall potential 
difference of about 1200 V as a reasonable target (see Fig. 
3-57). The Q3D calculation CMSIQBEIKW was made with I, * 850 
A, nominal wall temperatures, and diffuse current transfer. 
The predicted Hail potential shows good qualitative agreement 
with the data. Further computations at an increased would 
probably lower the potential difference to th._ larget value. 



Figs. 3-59, 3-60 and 3-61 are plots of the experimental 
and the calculated distributions of normal current density, 
Faraday voltage, and power respectively. The meter data 
distributions were obtained from an AEDC figure, which is 
reproduced as Fig. 3-62, while those for the data acquisition 
system were taken from (3-9]. Again, there is good qualitative 
agreement between theory and experiment for these quantities. 
Clearly, the theory and experiment lack quantitative agreement, 
especially in the rear-end region of the generator. Possible 
causes of this lack of agreement are (1) failure to match I^, 
(2) unreliable electrical data, or (3) the presence of 
phenomena not revealed by the present computation. 

Certainly, some of the discrepancy is expected since the 
Hall potential comparison of Pig. 3-57 indicates that a higher 
should be used in the computation. in addition, there are 
discrepancies in the experimental data recorded by the two 
measuring systems. Perhaps this is best illustrated in the 
current and voltage graphs of Figs. 3-59 and 3-60. While the 
voltage measurements at ?2 = 8 s agree fairly well, the current 
distributions vary differently along the channel. In addition 
to the fact that the data is in some instances perhaps not yet 
well understood, the occurrence of additional phenomena such as 
internal shorting and current loops (resulting in nonconstant 
I^) in the generator cannot be excluded. STD computation 
CMSIQBEIKW utilized a leakage current of I_ = 850 A to predict 

A 

a Hall potential difference of approximately 1400 V. An 
additional estimate for the magnitude of can be obta* ed 
from the equation 


I 


X 


Vv 

A 


L 
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where L is the length over which the short occurs and A is the 

cross sectional area occupied by the current. At the back end 

2 

of the channel a typical cross sectional area is 0.8 m . The 
measurements from the data acquisition system at - 8 s 
indicate an anode to ground potential (V) of approximately 1400 
V at 6.22 m from the first loaded electrode. The axial 
distance from this station to the diffuser (grounded) is 2.67 
m. STD computations yield an electrical conductivity ( ) in 
this region of approximately 3 S/m. Assuming the leakage 
current fills the entire cross-sectional area of the channel, 
application of Eq. (1) yields - 1300 A. Similarly, at the 
front end of the generator typical values of the cross- 
sectional area and gas conductivity are 0.14 m and 11 S/m, 
respectively. The experimental data at Tg * 8 s indicates an 
anode to ground potential of approximately 400 V at 1.82 m from 
the channel inlet. Using these values, application of Eq. (1) 
predicts a leakage current of 340 A. 

To obtain a valid value of the leakage current and its 
effects, it is necessary to take into account the structure and 
location of the leakage current with sophisticated models. 
Although the bulk conductivity *'veraging presented in the 
previous paragraph is crude, the order of magnitude of the 
leakage current may be estimated. 

3.7 Steady State Q3D Simulation of Run 006-014 

AEDC Run 006-014 was a powered run of 8 seconds duration 
which arced in the breech at approximately 5.3 s after the seed 
entered the channel. The experimental data show two distinct 
performance regimes (see Figs. 3-63 to 3-72) 

(1) A period between 4.5 £ Tg £ 5.21 s in which 
the channel is ungrounded and the data is 
relatively steady in time. 


36 



(2) A period between 5.5 ^ T.^ ^ 6.5 s in which 
the channel is operating in a grounded condition. 


The Q3D simulations performed under this contract have 
pertained only to the channel performance when in the 
ungrounded condition. 

The input data for the Q3D simulations are fully 
described in Section 3.1.2. These data were obtained by 
averaging the experimental data over the time period 4.54 ^ Tg 
< 5.21 s. The effect of the bleed resistor was accounted for 
in the calculation with a constant axial current of 23.54 A. 
This value was obtained from the average diffuser-to-ground 
voltage (9416 V) and the known resistance of the #12 wire used 
as the bleed resistor (400 ohms). 

Heat transfer measurements by AEDC personnel have 

resulted in a modified electrode surface temperature as shown 

in Fig. 3-6. The sparseness of the heat transfer data 

prohibits exact specification of the electrode surface 

temperature. The Q3D simulations have, at the suggestion of 

AEDC personnel, used temperatures corresponding to the 450 
2 2 

Btu/ft -sec and 50 Btu/ft -sec curves (at T = 5 s) at the inlet 
and exit flanges respectively. A linear variation of 
temperature with axial location was assumed between these 
po ats. The nominal temperature distribution was retained for 
the insulating walls (Fig. 3-4). 

A list of the nominal load resistances, supplied by 
AEDC, is contained in Table 3-7 for runs subsequent to Run 
006-006. These values of resistances have been normalized a 
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nominal water conductivity factor of 150 Ks/cm. Also listed in 
Table 3-7 ai e the average load resistances obtained from the 
experimental data at these locations. To conform as closely as 
possible to the experimentally observed values of the load 
resistance, and yet retain the detailed electrode-by-electrode 
distribution provided by the nominal schedule, STD adjusted the 
nominal schedule with a conductivity factor of 159 ps/cm for 
these initial calculations. This agrees with the conductivity 
values derived from the data acquisition system values for 
Faraday voltage divided by Faraday current over the first 168 
electrode pairs. The load schedule used in these simulations 
appears to be at variance with the experimentally indicated 
values in the downstream half of the generator, and the 
performance predictions may be expected to be less reliable 
toward the aft end of the generator. A plot of input load 
resistance schedule and the resistance as determined by the 
experimental data at each instrumented electrode pair is given 
in Fig. 3-86. 

Three Q3D calculations have been made during the 
theoretical simulation of Run 006-014: 

(1) STD Computation BTAYUYDJBE: assumes all walls 
prior to inlet flange can be characterized as 
smooth walls. Inlet boundary layer thickness 
is 18 mm. Enthalpy loss upstream of the 
generator inlet flange is 5.7%. 

(2) STD Computation CHPQUUUJBO: assumes roughness 

of all walls prior to generator inlet flange 
can be characterized by an equivalent sand 
roughness of 3 mm. Enthalpy loss upstream of 
generator inlet flange increases to 7.2% and 
inlet boundary layer thickness increases to 44 
mm . ( See Sec tion 3. 1.2.5^. 
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(3) STD computation CHPQQFLJDV; Assumes roughness 
of all walls prior to generator inlet flange 
can be characterized as smooth walls. Inlet 
boundary layer thickness is 18 mm. Enthalpy 
loss upstream of the generator is increased to 
7.2% to account for a 3% combustion 
inefficiency due to incomplete droplet-burnout, 
as discussed in [3-14]. The chemistry was 
maintained the same, however. 

Figures 3-74 to 3-89 contain the results of a typical 

Run 006-014 simulation calculation (CHPQUUUJBO) in the form of 

axial profiles of the key gasdynamic and electrical variables. 

Fig. 3-75 shows that the Faraday load factor, K, decreases 

along the channel from approximately O.I to 0.6. Figs. 3-76 

and 3-77 show that the interaction parameter based on velocity 

and the interaction parameter based on pressure achieve values 

at the end of the channel of approximately 2.3 and 0.95 

respectively. This interaction is sufficiently high to cause 

noticeable electrode boundary layer asymmetries (e.g., the 

electrode shape factors in Fig. 3-79) and nearly zero blockage 

at the channel exit (Fig. 3-76). Fig. 3-86 shows that the 

generator power output under these conditions is approximately 

22 MW. Maximum normal current densities and Hall fields of 

2 

approximately 0.7 A/cm and 2.0 kV/m respectively exist at 
steady-state conditions (Fig. 3-87). Fig. 3-89 indicates that 
boundary layer voltage drops in excess of 2000 V may exist at 
the end of the generator for this relatively cold wall 
operation. 

The boundary layer voltage drop measurements are perhaps 
not sufficiently refined at the present time to discriminate 
between various models of current transport at the walls. The 
diffuse discharge model selected for the Q3D simulation of Run 
006-014, has yielded voltage drops within the range of 
uncertainty of the data and has been approximately confirmed 
with the more sophisticated FIN and ARRAY calculations 
described in Section 4.5. 
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3.8 Steady-State TRANSIENT Simulation of Run 006-014 

Under conditions of higher interaction than were 
studied with Q3D in the neighborhood of the Run 006-014 
operating conditions, the flow will become transonic in the MHD 
generator. The Q3D family of codes requires additional 
computational resources to iterate on the downstream pressure 
condition in transonic or subsonic flows. To economize, the 
high interaction cases for the Run 006-014 parametric study 
were carried out with the TRANSIENT family of codes, with the 
quasi-steady option. The starting point for the TRANSIENT 

parametric variations is a simulation of Run 006-014 utilizing 
data generated by the Q3D simulationn CHPQUYDJBE to account for 
multidimensional phenomena. Figures 3-90 through 3-98 are 
plots of the results of a computation, BRCYTRCJEB, with a code 
from the TRANSIENT family to simulate AEDC Run 006-014. The 
origin of the axial coordinate in these graphs is the burner 
backplate. The first loaded electrode is at x = 2.7 m and the 
last is at X = 9.4 m (see Fig. 3-90). These calculations were 
made for conditions similar to those of CHPQUYDJBE (smooth 
walls upstream of generator and 100% efficient combustion), 
which were discussed in the previous section. 

Most of the key parameters resulting from TRANSIENT 

simulations shown in Figures 3-90 through 3-98 are also quite 

comparable to the Q3D simulation CHPQUUUJBO, which was 

described in the previous section. In particular. Figure 3-96 

2 

shows that a maximum current density of 0.685 A/cm and a 

maximum Hall field of 2.005 kV/m are predicted with the 

TRANSIENT code. This compares well with the J of 0.7 

2 y,max 

k/cfsT and the E^^ of 2.007 kV/m predicted with Q3D (see Fig. 

3-87). Similarly, the exit Mach numbers predicted by TRANSIENT 

and Q3D are 1.73 and 1.74, respectively (Figs. 3-78 and 3-91). 
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These figures also show the shock is located 10.74 m from the 
burner backplate. The total heat loss to the wall predicted by 
TRANSIENT is 31. 4t MW, while that predicted by Q3D is 30.97 MW. 
The TRANSIENT calculation yields 23.6 MW for the integrated 
power output while the Q3D calculation CHPQUUUJBO predicts 21.3 
MW or a difference of approximately 10%. Most of this 
discrepancy is due to the pre-generator wall roughness 
characterization which was for smooth walls in the TRANSIENT 
calculations and for rough walls in CHPQUUUJBO. The power 
output from the TRANSIENT calculation is more comparable to the 
result of the Q3D simulation CHPQUYDJBE, 24.6 MW. 

The simulation of the generator operation with the 
TRANSIENT code is useful because it gives a broader view, if 
with slightly less detail, of the entire flow train. It should 
be noted that the accuracy of TRANSIENT and the agreement 
between results from computations with the TRANSIENT and Q3D 
code families is due to inputs to TRANSIENT from Q3D 
calculations to account for fundamentally multidimensional 
phenomena such as wall losses and plasma nonuniformity factors. 
This is necessary to ensure the accuracy of any quasi-one- 
dimensional calculation. 

3.9 Comparisor of Run 006-014 Measurements with Q3D 
Simulations 


The results from Q3D computations discussed in Section 
3.7 were analyzed and compared to the experimental data from 
HPDE Run 006-014. In this section the results of this 
comparison are presented with emphasis placed on two Q3D 
calculations, CHPQUUUJBO and CHPQQFLJDV. The conditions and 
underlying assumptions for these separate computations (see 
Section 3.7) are identical except for the following: 
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(1) The Inlet boundary layer thickness and enthalpy 
flux are determined by a smooth wall boundary layer 
calculation upstream of the generator for 
computation CHPQQFLJDV (rather than rough with 3 mm 
equivalent sand roughness assumed in CHPQUUUJBO) . 

(2) A 3% combustion inefficiency due ^o incomplete 
droplet-burnout is assumed for computation 
CUP^FLJOV (rather than a 100% efficient combustor 
assumed in CHPQUUUJBO) . 

The experimental data for the static pressure, a primary 
measurement of the data acquisition data, during the time after 
ignition, 4.1 jC Tg £ 5.2 s is given in Fig. 3-99. The static 

pressure is a primary measurement of the data acquisition data. 
The axial coordinate in Fig. 3-99 is referenced to the 
generator inlet flange. Also given in Fig. 3-99 are the axial 
distributions of static pressure from computations CHPQ.<UUJBO 
and COBQQFLJDV. There is excellent agreement between the 
experimental and the computed values throughout the channel. 

The experimental HPDE meter panel data from Run 00 014 

at Tg * 4.8 s are compared to the values obtained from 

computations CHPQUUUJBO and COBQQFLJDV in Figs. 3-100 through 
... xll. The axial coordinate in these figures is referenced to 
the inlet flange. The Faraday vc'.tage (Fig. 3-100) and the 
current per electrode (Fig. 3-101) are primary measurements of 
the meter panel data. The normal current density (Fig. 3-102) 
and the integrated power output (Fig. 3-105) are determined 
from the primary measurements and other geometrical data. The 
computed values from computations CHPQUUUJBO and COBQQFLJDV 
fall within the measurement tolerances of the experimental 
values. Some local discrepancies of up to 5% are evident, but 
it is our Judgment that the accuracy of the input data, 
especially for load resistance, as well as the + 2% fluctua- 
tions in current and voltage during the period 4.5 s < Tg < 
5.21 s could explain such differences. As shown in Fig. 3-100, 
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the computations tend :o underpredict the Faraday voltage in 
the front half of the generator. This can bf explained by 
consideration of the uncertainty of the channel loading. 

The difference between the voltage-current characteri- 
stics measured by the two data-ga^hering systems and the 
nominal values is discussed in Section 3.7. Tb,^ calculatec 
values of the load resistances from the meter panel data yield 
values that are typically 10% higher than those computed from 
the data acquisition system (see Fig. 3-103). This is 
especially true in the latter third of the generator. In 
addition, the value of a load resistance computed at a 
particular electrode pair can vary 3-5% with time. 

A rough calculation, under the assumption that the fluid 
properties change only slightly at this rather low interaction 
level , indicates that if the downstream load resistance were 
matched more accurately, a lower current and a higher Faraday 
voltage could be expected, the differences being approximately 
3 '5%. The integrated power output, shown in Fig. 3-105 would 
be relatively unaffected (the difference is estimated to be 
less than 1%) by this adjustment to Input load resistance. 

The integrated power distribution for the meter panel 
data measurements can be regarded as approximate, but better 
than the distribution obtained from the data acquisition system 
measurements. Ihese points were calculated by numerically 
integrating the power through the loads for the instrumented 
electrode pairs. To obtain a val"e for thy integrated pow..r, 
the power was assumed to vary linearly between each set of 
instrumented electrodes (10 electrode pairs apart). The 
scatter In the power per electrode measured by the meter panel 
data, and plotted in Pig. 3-10', indicates that such a regular 
behavior is probably a weak assumption. 
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Computation CUPQUUUJBO simulated a bleed resistor fr(»n 
the diffuser-to-ground with a net axial current of 24 A. The 
power dissipated ia this bleed resistor was calculated to be 
0.2 MW by CHPQU0UJBO. This ic less than 1% of the power 
produced in the channel. 

It should be reemphasized that these preliminary 
computations have assumed diffuse-mode current transfer. At 
t’’e relatively cold electrode surface temperatures indicated by 
the experidtental heat transfer measurements to this point, arc 
mode current transfer is a distinct possibility, particularly 
on the anodes, and should be considered. This will be 
discussed further in Section 4.5 This will be discussed 
further in Section 4.5. The results described in Section 4.5 
further suggest the possibility of the existence of very small 
arcs, particularly on the anodes. 

The axial distribution of total electrode voltage drop 
from computation CHPQUUUJBO is given in Fig. 3-106. 
Experimental values of total voltage drop, obtained from the 
HPDE program manager, at four axial locations are also given in 
this figure. The electrode voltage drops computed by STD 
Computation CHPQUUUJBO appear to agree with the measured 
voltage drops within the experimental tolerances. 

The axial distributions of centerline, cathode, and 
anode electrical potential are given in Figs. 3-107, 3-108, and 
3-109, respectively. The potentials have been adjusted to set 
ground (0.0 V) at the centerline of the first loaded electrode 
pair. There is excellent agreement between the meter panel and 
the results of computation COBQQFLJDV. Also given on Fig. 
3-109 are the data obtained by the data acquisition system for 
the time period 4.1 £ Tg _< 5.2 s. 
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TABLE 3-1 


Chanaei Electrode Width 


Electrode Reference 
No. ♦ 


Electrode Width** No. of Electrodes 

(m) 


-30 

to 

26 

0.0254 

57 

27 

to 

58 

0.01905 

32 

59 

to 

162 

0.01588 

104 

163 

to 

394 

0.0127 

232 

395 

to 

444 

0.0254 

50 

445 

to 

454 

0.02223 

10 


*Electrodes are numbered consecutively upstream and downstream 
from electrode 0 which is the first electrode connected to a 
load. 

** The insulator thickness is 1.52 nmi between adjacent 
electrodes 


0-4355 



TABLE 3-2 


Nominal Flow Rates of Combustion Gases 


Component 


Flow Rate 


Fuel 

Oxidizer 

Seed 

N2/O2 

Pilot (CH4) 

Seed Concentration 


5.641 kg/s 
37.247 kg/s 
2.069 kg/s 
1.0 

0.034 kg/s 

0.454 kg KOK/kg CH3OH 


Total enthalpy of combustion gases (no losses) 271.7 MW 

Specific enthalpy of combustion gases 6.04 MJ/kg (adiabatic 
conditions) 

■Aflame * ^°21 K 


0>4356 
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TABLE 3-3 


Nominal Magnetic Field Used for AEOC Simulations 


Run no. 

2 

Time Interval 

Ref. 

006-006 

1.537 

4.5 

< Tg < 

7.5 

13-81 

006-008 

2.327 

9.5 

< Tg < 

10.2 

13-91 

006-010 

2.697 

4.0 

< ^2 < 

4.5 

[3-lOj 

006-013 

2.810 

4.5 

i Tg < 

5.5 

13-111 

006-014 

3.231 

4.5 

< Tg < 

5.3 

13-111 


♦Due to nonlinear bebavior of the instrumentation, the numbers 
presented in this table are slightly in error. Refined data. 
13-12 j. were not available at the time these simul:;tions were 
performed. 


0-4357 
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TABLE 3-4. 


A. HPDE Load Resistance Settings for Runs Prior to 006-008 
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Electrode 

No. 

Reference 
Res is. fi 
(K = 150) 

Date 

1979 


Electrode 

No. 

Reference 
Resis. R 
(K = 150) 

Date 

1979 

58 

o8 

Oct 

.11 


102 

53 

Oct 

.15 

59 

47 




103 

51 



60 






104 





61 






105 


f 



62 






IOC 

53 



63 






107 

56 



64 

50 




108 

58 



65 

47 




109 

53 



66 

50 




110 

51 

1 

f 

67 

47 




111 

50 

Oct 

.16 

68 






112 

50 



69 






113 

52 



70 






114 

50 



71 

50 




115 





72 

47 




116 





73 






117 

J 




74 






118 

1 

' 



75 






119 

52 



76 






120 

50 



77 






121 





78 






122 





79 






123 





80 

50 




124 





81 

51 

Oct 

.15 


125 





82 






126 





83 






127 

i 




84 






128 

52 



85 






129 

50 



86 






130 





87 

53 




131 





88 

51 




132 





89 






133 

52 



90 






134 

52 



91 




I 


135 

50 



92 

53 




136 





93 

51 




137 





94 






138 





95 






139 





96 






140 





97 






141 





98 






142 

48 


1 

99 






143 

52 



100 






144 

50 



101 




L 


145 

52 

1 

f 


0 - 4330*1 


51 









Reference 
Resis. R 
(K * 150) 


Date 

1979 


Reference 
Resis. R 
(K = 150) 


Date 

1979 


Electrode 1 
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TABLE 3-4. (cont. 


B. HPDE Load Resistance Settings for Runs 
Subsequent to 006-007 


Electrode 

No. 


Reference 
Resis. Q 
(K » 150) 


Open (®) 




















TABLE 3-5 


Water Conductivity in Load Buckets for AEDC Experiments 


Run No. 

Krun 

Ref. 

006-006 

i78.0 

13-8] 

006-008 

175.4 

13-9] 

OOG-010 

178.6 

13-11] 

006-011 

approx. 183.0 

13-11] 

006-012 

172.4 

13-11] 

C06-013 

150.0 

13-15] 

006-014 

160.0 

13-16] 


0-4358 
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TABLE 3-6 


Plow Rates of Combustion Gases and 
Runs 006-008 and 006-014 


Component Run 006-008 

Fuel (kg/s) 6.60 

Oxidizer (kg/s) 40.99 

Seed ^kg/s) 1.84 

Pilot (CH^) (kg/s) 0.034 

N 2 /O 2 ratio 1.0 

Seed Concentration 0.3 

(kg KCH/kg CH 3 OH) 

Equivalence Ratio 
(fuel lean) 1.5 

Total enthalpy of combustion 
ga:>es (no losses) (MW) 299.3 

Specific enthalpy of 
combustion gases (HJ/kg) 6.05 

Tflame <^) 3048 


0-4359 


006-014 

6.84 

39.82 

2.54 

0.034 

0.929 

0.3 

1 . 11 

303.8 

6.16 

3063 
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TABLE 3-7 


Load Resistances for AEDC Run 006-014 


Data 

Acquisition 

Station 

Number 

Electrode 

Pair 

Nunber 

Nominal 
Resistance 
» 150) 
^^ohms) 

Averaged 
Experimental 
Res j stance 
4.54 1 T2 £ 5.21 
(ohms) 

Indicated 

Conductivity 

Factor 

(pS/cm) 

1 

42 

39.1 

36.64 

160.0 

2 

84 

50.9 

48.19 

158.4 

3 

126 

50.0 

47.97 

156.3 

4 

168 

69.1 

65.63 

157.8 

5 

210 

61.9 

61.25 

151.6 

6 

252 

67.8 

66.02 

154.0 

7 

294 

65.5 

65.82 

149.3 

8 

336 

65.5 

65.70 

149.5 

9 

376 

141.6 

153.71 

138.2 


0-4360 
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J'iK. 3-1. AKDO/HPDE internal geometry with designated cross sections 






9-3187 


Fig. 3-2. Axial variation of the channel geometry in the 
AEDC/HPDE at the nominal 2 T operating condi- 
tions. 
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Fig. 3-3. Magnetic field strength distribution along magnet bore 
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VM>OR BARRIER 



Location of the Hall probes on the floor of the magnet bore 
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Computed time history or pressure distribution throughout the 
flow train of typical 300: MWth iMHD generator during start-up. 
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Electrical power output as a function of mass flowrate and 
magnetic field for the nominal operating conditions of the 
AEDC/HPDE. Code families: TRANSIENT and Q3D. Codes: :QS 

(solid symbols) and Q3DXY > Q3DXZ /pen symbols), hll 
computations completed before 18 September 1979* 

STD RtSfARCH corporation* 
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Fig. 3-12. Axial variation of the channel geometry in the AEDC/ 
HPDE at the nominal 2 T operating conditions; Q3D 
calculation . 
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Fig. 3-13. Axial variation of the centerline gasdynamic va-iables 
in the AEDC/HPDE at the nominal 2 T operating condi- 
tions; Q3D calculation. 


73 











o| o 
id! • 


N 

O 

Sg ^ 

UJO ^ • 


a: o 

tij « 


lyg 

i: >0. 

CEW >-•• 
€C K 


-e >-8 

U V 
CCM 39 
X o 


lU o 
z u 
1-^ o 
-I*? Ul®. 
Oo 

Ul oc 
»- QC 
Z liJ 
lU > 
Oo <lo 
o o 


CO CM 


ol o 
ol o 



9-3192 


X (M) 


Fig. ?-14. Axial variation of the plasma parameters in the AEDC/ 
HPDE at the nominal 2 T operating conditions: Q3D 

calculation . 
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Fig. 3-16. Axial variation of the power extraction parameters in 
the AEDC/HPDE at the nominal 2 T operating conditions: 
Q3D calculation 
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Fig. 3-17. Axial variation of the electrical field variables in 

the AEDC/HPDE at the nominal 2 T operating conditions: 
Q3D calculations. 
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Fig. 3-18. Axial variation of the short circut fields and the 
local internal impedance in the AEDC/HPDE at the 
nominal 2 T operating conditions: Q3D calculation. 
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Fig. 3-19. Axial variation of the electrical stresses on the 
electrodes in the AEDC/HPDE at the nominal 2 T 
operating conditions: Q3D calculation. 
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Fig. 3-20. Axial variation of the electric potential in the AEDC/ 
HPDE at the nominal 2 T operating conditions: Q3D cal- 

culation . 
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Fig. 3-21. Axial variation of the electrode voltage drops in the 

AEDC/HPDE at the nominal 2 T operating conditions: Q3D 

calculation . 
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Fig. 3-22. 


Axial variation of the electrode boundary layer para- 
meters in the AEDC/HPDE at the nominal 2 T operating 
conditions; Q3D calculation. 


82 





S. w 


ntS£*«CM CORPOfUtlON 


O O O O O M 

»C tft, Nj 


flEOC-NPOE NOrUNIL 2T CONDiriSNS / COPOrKlllHIi 


2j 8j J V 


gji gj7 gJ 


. 8 




5 £ 


O. 

W 

au 

r ^ 

2^10 

»-iC3L =OL 

t 2C 
Co. x:^ 

X* 5* 


>C(D Ho 

OnOL r£ML 

<3C • 

:5 2 

i 2 


S S 

, X 

i#> Oi0 

rjou c*^ 

s :i 

2 i 

«.5o 

o. 


sj-sj 


•* Ui 

z 

O Zu 

« . Jr . 

Z 

o “ 


z ^ 

safe 

i* 

(X n 


— 1» -!« 
o*> 

- 2 
= 2 
5o 

o«^ 


(0 \tn 


oj oj oj 


0| o 
o| D 


OI O 
Of O 



-3193 


Fig. 3-23. Axial variation of the sidewall boundary layer para- 
meters in the AEDC/HPDF. at the nominal 2 T operating 
conditions: Q3D calculation. 



Axial variation of the heat transfer parameters «in the 
AEDC/HPDE at the nominal 2 T operatine conJit'ons: Q3D 

calculation . 
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Fig. 3-25, Anode to ground voltage of electrode 336 for AEDC/HPDE 
run 006-008 
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Fig. 3-28. Total Hall voltage for AEDC/HPDE Run 006-008. (Dashed 
lines represent vo tage history to ground.) 
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Fig. 3-32. Axial variation of the centerline gasdynamic variables 
from computation CMSIPUFIKW for the AEDC/HPDC. 
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Fig. 3-34. Axial variation of the short circuit variables and 
internal impedance from computation CMSIPUFIKW for 
the AEDC/HPDE. 
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Fig. 3-35. Axial variation of the electrical loading 

from computation CliSIPUFIKW for the AEDC/HPDE. 
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from computation CMSIPUFIKW for the AEDC/HPDE. 
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fig. 3-37. Axial variation of the electrical field variables from 
computation CMSIPUFIKW for the AEDC/HPDE. 
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Fig. 3-40. Axial variation of the heat transfer parameters from 
computation CMSIPUFIKW for the AEDC/HPDE. 
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Axial variation of the plasma properties from 
computation CMSIQBEIKff for the AEDC/HPDE, 
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Fig. 3-48. Axial variation of the power extraction parameters 
from computation CMSIQBEIKW for the AEDC/HPDE. 
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Fig. 3-51. Axial variation of the electrode voltage drops from 
computation CMSIQBEIKW for the AEDC/HPDE. 
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Fig. 3-52. Axial variation of the heat transfer parameters from 
computation CMSIQBEIKW for the AEDC/HPDE. 
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Fig. 3-53. Axial variation of the electrode boundary layer 
parameters from computation CMSIQBEIKW for the 
AEDC/HPDE. 
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Fig. 3-54. Axial variation of the sidewall boundary layer 
parameters from computation CMSIQBEIKW for the 
AEDC/HPDE. 
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Fig. 3-55. Axial variation of the channel geometry from 
computation CMSIQBEIKW for the AEDC/HPDE. 
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Fig. 3-»5S • Anode and cathode potential Run 006*008 
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Fig* 3—63 • Pow^r HPDE Bun 006—014 J5TD R£$eA«CM coapohation 
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Fig. 3>69. Voltage anode to ground HPDE Run 006-014 enm 
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Fig. 3-73. Load resistance for HPDE Run 006-014 
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Fig. 3-74. Channel geometry HPDE Run 006-014 


134 


ORIGINAC PACE f« 

OF poo:i c’’ • ' rv 


8 

# 



0-4136 



ELECTRICAL LOROINO 

HPOE RUN OOS-014 



X (N) 


Fig. 3-75. Electrical loading HPDE Run 006-014 
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Fig. 3-77. Interaction parameter based on pressure HPDE Run 
006-014 
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Fig. 3-78. Centerline gasdynaniic variables HPDE Run 006-014 
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Fig. 3-79. Electrode boundary layer parameters HPDE Run 
006-014 
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Fig. 3-81. transfer paraicters HPDE Run 006-014 
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Fig. 3-82. Plasma properties HPDE Run 006-014 
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Fig. 3-83. Nonuniformity factors HPDE Run 006-014 
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*ig. 3-84. Short circuit variables and internal impedance 
HPDE Run 006-014 
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Fig. 3-85. Electrical stresses on the electrodes HPDE Rur 
006-014 
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Fig. 3-88. Electric potential HPDE Run 006-014 
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Fig. 3-89 . Electrode voltage drops HPDE Bun 006-014 
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Fig. 3-90 . Flowtrain geometry for AEDC Run 006-014 







Fig. 3-91. Axial distribution of average gasdynamic variables for 
Run 006-014 
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Fig. 3-93. Axial distribution of boundary layer heat transfer paraiueters 
for .4EDC Run 006-014 
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Fig. 3-94. Axial distribution of plasma properties for AEDC Run 006-014 
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-95. Axial distribution of electrical loading for AEDC Run 006-014 
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Fig. 3-99. Comparison of the axial distributions of static pres 
sure from the data acquisition system and STD Simula 
tions of HPDE Run 006-014. 
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Fig. 3-100. Comparison of the axial distributions of Faraday volt- 
age from the experimental data and from STD simulations 
of HPDE Run 006-014. 
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Fig. 3-101 • Comparison of the axial distributions of current per 
electrode from the experimental data and STD simula- 
tions of HPDE Run 006-014. 
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Pig. 3-102. Comparison of tMf. axial distributions of normal current 
density from tl i experimental data and STD simulations 
of HPDE Hun OOt-014. 
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Fig. 3-103. Comparison of the axial distributions of load resist- 
ance from the experimental data and STD simulations of 
HPDE Run 006-014. 
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Fig. 3-104 


Comparison of the axial distributions of power per 
electrode from the experimental data and STD simula- 
tions of HPDE Run 006-014. 


X + 



ACSCAPCN CORPORATION 


Q POWER EXTRACTION PARAMETERS 

. HPDE RUK P06-OI4 



Fig. 3-105. Comparison of the axial distributions of integrated 

power output from the experimental data and STD simu- 
lations of HPDE Run 006-014. 
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Fig. 3-106. Comparison of the axial distributions of 2-D electrode 
voltage drop VD(V) for HPDE run 006-014. 
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Fig. 3-107. Comparison of the axial distributions of centenine 
potential from the experimental data and STD simula- 
tions of HPDE Run 006-014. 
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Fig, 3-108. 


'i.iparison of the axial distributions of ca,hode (STD 
. ion) potential from t’.ie experir.iental data and 
Old .^mutations of HPDE Pun 006-014. 
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Fig. 3-109. Comparison of the axial distributions of anode (STD 
notation) potential from the experimental data and 
STD simulations of KPDE Run 006-014. 
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Fig. 3-110. Comparison of the axial distributions of anode inter- 
electrode voltage from the experimental data and STD 
simulations of HPDE Run 006-014. Experimental data 
represent averages over 10 electrodes. 
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Fig. 3-111. Comparison of the axial distributions of the Anode 
Hall field from the experimental data and STD simu- 
lations of HPDE Run 006-014. Expermiental data re- 
presents averages over 10 electrodes. 


4.0 CRITICAL PHENOMENA IN THE AEDC/HPDE 

The experimental results of the AEDC/HPDE to date 
reveal the importance of a number of performance-controlling 
mechanisms in large-scale, high-interaction HHD power 
generators. Among these are boundary layer voltage drops and 
plasma nonuniformity, near-electrode phenomena including finite 
electrode effects, end effects, and the effects of various 
operating parameters on the optimization of the HPDE 
performance. This section is devoted primarily to the 
elucidation of these effects through comparisons with and 
parametric variations from the experimental operating 
conditions of HPDE Run 006-014. 


4.1 PARAMETRIC VARIATIONS AROUND THE Q3D SIMULATION OF RON 

006-014 


The successful completion of the Run 006-014 
simulations gives an opportunity to exploit the computer models 
to predict performance in other operating regimes. Using the 
results of the detailed Q3D simulation of HPDE Run 006-014 
described in Section 3.7 as a starting point, single-parameter 
variations in wall temperature, load factor, N/0 ratio, mass 
flow rate, and magnetic field have been carried out. These 
computations demonstrate the sensitivity of HPDE performance to 
most of the operating conditions under the control of the 
experimentalists. 

Fig. 4-1 presents the results of the parametric 
computations. All of the computations represented in Fig. 4-1 
are variations from computation CHPQUUUJBO, which is the 
nominal, simulation of HPDE Run 006-014. The Faraday load 
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factor distribution was altered to a constant value of 0.55 
(computation COBQFRQJEA), 0.60 (computation COBQAVWJDY), and 
0.70 (computation COBQBGDJDY . The electrode wall surface 
temperatures were increased to a constant value of 1250 K 
(computation COPQRPHJBJ), to a constant value of 1500 K 
(computation COBQJQKJDY), and to a constant value of 1750 K 
(computation COBQDTBJDl ) . Also shown is computation 
COBQAIYJEl, in which the electrode wall surface temperatures 
were assumed to rise to ”20-sec" values, which were determined 
by the same method applied to obtain the "5-sec" values of the 
nomical simulations. The peak magnetic field strength was 
increased to 4.02 T (computation COBQLRNJDZ) and to 4.28 T 
(computation COBQIJOJEA). The N/0 ratio was decreased to 0.80 
(computation COBQMLZJDZ) and increased to 1.00 (computation 
COBQIFBJEA). In contrast to the other parametric computati ~>ns, 
the N/0 variations were centered about computation COBQAVWJDY 
with a constant load factor of 0.60. 

These new computations have shown that a constant 
Faraday load factor of 0.60 yields slightly better performance 
than the nominal simulation, which has a varying load schedule 
with an average load factor of approximately 0.61. A constant 
Faraday load factor of 0.70 yields worse performance. A load 
factor of 0.55 increases enthalpy extraction at the expense of 
isentropic efficiency. 

As shown in Fig. 4-1, an increase in wall temperature 
to 1500 K results in higher values for both enthalpy extraction 
and isentropic efficiency than for the nominal simulation or 
the simulation with 20 second wall temperatures, but lower 
values than for 1750 K walls. 
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An increase in mass flow rate to 55 kg/s results in 
lower isentropic efficiency and lower enthalpy extraction 
ratio. While a computation with decreased mass flow rate was 
not carried out, it is clear that some improvement in the 
performance would probably be obtained by dropping the mass 
flow rate under the Run 006-104 conditions. The reason is that 
a reduction in the mass flow rate would drive the Mach number 
distribution toward the optimum, which in the case of the 
AEDC/HPDE appears to be transonic (see Sec. 4.4). 

An increase in magnetic field strength (COBQLRNJDZ) 
results in higher values for both generator isentropic 
efficiency and enthalpy extraction ratio. 

A variation of N/0 ratio to 0.80 leads to an increase 
in enthalpy extraction ratio, but a decrease in isentropic 
efficiency. The opposite is observed with an increase in N/0 
ratio to 1.00. 

It is important to distinguish the genf'rator 
isentropic efficiency displayed in Fig. 4-1, which expresses 
how efficiently the generator made use of the total pressure 
drop across the generator, with the power train isentropic 
efficiency, which expresses how efficiently the generator made 
use of the pressure drop from the combustor to atmosphere. The 
latter quantity varies almost directly with the enthalpy 
extraction rate and is perhaps less informative than the 
generator efficiency. The results in Pig. 4-1 are easily 
interpreted if this distinction is kept in mind. 

The conclusion from this survey is that, at least in 
the neighborhood of the Run 006-014 operating point, the 
pay-off from providing hot walls in this experiment is 
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significant and Justifies a change in operating procedures or 
hardware to achieve elevated electrode surface temperature. It 
Is apparent that increases of the magnetic field and electrode 
surface temperatures are the primary means of attaining both 
the enthalpy extraction and Isentropic efficiency goals of the 
HPDE. 


There are methods for mitigating the voltage drops, 
even in cold-wall designs (e.g., load schedule optimization). 
The sensitivity of the electrode voltage drops to the 
parameters studied thus far are described in the following 
section. 
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4.2 Parametric Dependence of Electrode Voltage Drops and 
Plasma Uniformity 

The Q3D parametric survey of performance around the 
Run 006-014 conditions provides opportunity to investigpte the 
sensitivity of the electrode voltage drops and the plasma 
nonuniformity to the important parameters of the experiment. 
In this section, the wall temperature, magnetic field, and load 
factor dependence of the total electrode voltage drop and two 
parameters, which characterize the effects of the plasma 
uniformity on generator performance, are examined. Before 
investigating the sensitivity of the total electrode voltage 
drop and the two nonuniformity factors to these three 
variables, it is useful to review the importance of the voltage 
drops and the nonuniformity factors in characterizing the 
performance of the MHD generator. 

The primary importance of the electrode voltage drop 
is that it is a measurable quantity. The electrode voltage 
drop is measured by subtracting the Faraday voltage from the 
total transverse voltage one would obtain by extrapolating the 
core electric potential gradient all the way to the walls. 
Clearly, high electrode voltage drops correspond to reduced 
power production. However, there is not a simple way of 
relating the magnitude of the electrode voltage drops to the 
reduction in power output. Nevertheless, because the electrode 
voltage drops are measurable and give an indication of the 
power penalty for given operating conditions, it is a useful 
parameter to compute and to compare with the experimental data. 
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The cootrlbuting factors to the electrode voltage drop 
may be easily seen In the simplified MHD electrical equations. 
We can express the electrode voltage drop in terms of the 
formula for the transverse electric field Eyi 


E. 


1+1 


2 


Jy + UB - pE^ 


(4-1) 


y <r 

It is seen from this formula that the voltage drops are equal 
to the difference between the centerline value of By and the 
cross-sectional average value of By, since the Faraday voltage 
is simply the average transverse electric field, 
multiplied by the transverse dimension of the channel, D. Thus 
the voltage drop is expressed as follows: 


''d • <V> 


(4-2) 


It is seen from this formula that the electrode voltage drop 
can be considered to be composed of three components. These 
components are: the current density times the effective 
resistivity difference of the cross-section. 



(4-3) 


the induced field deficit in the boundary layers 


[(U)^-(U>]b (4-4) 

and the variation of the Hall parameter through the electrode 
boundary layers due to pressure and temperature gradients in 
the boundary layers, 
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(4-5) 


[»)t-<P>]E, 

It Is perhaps Interesting to tabulate these components 

of the electrode voltage drops under various conditions. 

However, the variation of these expressions does not yield a 

great deal of Insight or quantitative Information about the 

effects of experimental parameters upon the overall performance 

of the generator. Again, this Is because the electrode voltage 

drop Itself Is not simply related to the performance of the MHD 

generator. Also, because J and E are factors In these 

y * 

expressions which are themselves affected by the voltage drops, 
these expressions can hardly be called Intrinsic 
characteristics of the channel or useful Indications of the 
performance penalties associated with the voltage drops. 

A much more useful quantity for expressing the effects 
of experimental parameters on the generator performance Is the 
plasma nonunlformlty factor. The plasma nonunlformlty factor 
Is defined as 

G » <<^>( , (■*-6) 

where the brackets Indicate cross-sectional averages of the 
quantities enclosed within. The concept of the plasma 
nonunlformlty factor Is a generalization of the original 
expression by Rosa [4-1]. 

In the case where variations of the Hall parameter In 
the cross-section are weak (usually a good assumption), the 
expression for the plasma nonunlformlty factor can be 
simplified to the following expression: 

G or g + (g-l)P^ (4-7) 
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where 


g = <o->(l/«r> 


is the conductivity nonuniformity factor. Thus, the plasma 
nonuniformity factor, G, may be separated into contributions 
from the nonuniformity of the conductivity, g, and the strength 
of the Hall parameter, which is directly proportional to the 
magnetic field and approximately inversely proportional to the 
static pressure. 

The plasma nonuniformity factor, G, is a useful 
parameter because it directly divides the power density at any 
given station in a Faraaay generator =0) 

, <T>(u/ B^l-K)K (4-9) 

The conductivity uniformity factor, g, is an intrinsic 
parameter which depends only on the distribution of 
conductivity in the channel cross-section, and not upon 
specific knowledge of the MHD interaction. 

The scale-dependence of g has been investigated by STD 
Research Corporation in [4-2]. There it was shown that a 
number of loss mechanisms, including cold boundary layers, 
axial current leakage, surface voltage drops, and other effects 
could be summarized in this simple parameter, g. Also, as the 
scale of the generator increases (with fixed wall conditions), 
the v?lue of g was shown to generally decline as the Inverse 
square root of the thermal input to the generator. In most 
applications, the dependence of the Hall parameter in the 
channel cross-sectlou is weak enough to make the parameter, g. 
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a very useful characterization of the potential performance of 
the BIHD generator duct. 

The following sections examine the sensitivity of V^, 
G, and g to variations in the electrode wall temperature, the 
magnetic field strength, and the load factor. It is to be 
understood that the effects of Hall parameter variation, 
induced field deficit, and effective resistivity variations are 
lumped into the total electrode voltage drop, <n the 
following sections. 

4.2.1 Dependence on Electrode Wall Temperature 

This section compares the total electrode voltage 
drops and the plasma and conductivity nonuniformity factors 
resulting from three parametric variations in electrode wall 
temperature around the nominal Run 006-014 simulation. The 
nominal calculation has been discussed in Section 3.7. The 
three computations carried out with electrode wall temperature 
parametrically altered to constant values of 1250 K, 1500 K, 
and 1750 K have been described in Section 4.1. The conclusions 
regarding the sensitivity of electrode voltage drops and plasma 
uniformity presented in this Section are based upon those four 
computations, and apply primarily to conditions in the 
neighborhood of the operating conditions of Run 006-014. 

To display the sensitivity of the electode voltage 
drops and the plasma uniformity, we have chosen to tabulate 
their values at four stations in the HPDE. Three of these four 
stations were locations at which the HPDE instrumentation 
provided transverse profile measurements. The fourth station, 
which varies about 0.1 m from computation to computation, was 
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the location at which the maximum electrode voltage drop 
occurred in each STD Q3D computation. 

The location of the three volta<?e profile measurement 
stations and the approximate location the maximum voltage 
drop computed by STD Research Corporation is shown m Pig. 4-2, 
which has been reproduced from (4-3 J. Fig. 4-2 also presents a 
useful definition of the various coordinate systems used by STD 
Research Corporation and AEDC in citing positions in the 
channel. The lower axis in Fig. 4-2, which has been marked 
"STD Q3D coordinates" uses as a reference point the inlet 
flange of the MHD generator. The STD coordinate system also 
begins at the inlet flange of the MHD generator, but is 
expressed in metres rather than in inches. A second frame of 
reference has as its origin the first loaded electrode. This 
frame of reference in shown in the upper axis which contains 
the location and stations numbers, in inches, of the electrode 
voltage drop measurements made at the HPDE The three stations 
at which the electrode voltage drops have been tabulated are, 
in STD coordinates, x = 1.43 m, x = 2.78 m, and x = 4.62 m. 
These locations correspond to STA 24.7 in., STA 77.4 in., and 
STA 215 in. 

Fig. 4-3 il lustr-.tes the calculated variation of the 
total electrode voltage drop at the four stations in the MHD 
generator as a function of electrode surface temperature. It is 
noted that the maximum electrode voltage drop declines from 
approximately 2100 volts to approximately 1200 volts as the 
electrode surface temperature increases from the "5 sec" value 
of approximately 600 K to the assumed temperature 1750 K. At 
the front of the generator, x * 1.43 m, the electrode voltage 
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drop declines from approximately 500 volts to approximately 3C0 
volts as the surface temperature decreases from then's sec"value 
of approximately 1000 K to the assumed wall temperature of 1750 
K- The decrease in electrode voltage dropr is approximately 
linear with temperature at each station. For comparison the 
experimentally measured electrode voltage drops at stations 
1.43, 2.78, and 4.62 m are displayed at the same wall 
temperature that was assumed for the Q3D simulation of Run 
006-014. The agreement between the measured and computed 
voltage drops appears to be within the experimental uncertainty 
of the measurements. 

As noted above, the impact of the electrode wall 
surface temperature on the performance of the generator can be 
suimnarized by the value of the plasma nonuniformity factor G, 
which results from the nonuniformity of the conductivity field 
due to cold boundary layers. Fig. 4-4 presents the average 
value of the plasma nonuniformity factor throughout the 
generator as a function of the average electrode wall 
temperature. The error bars surrounding each average point 
correspond to the standard deviation of the value of wall 
temperature or plasma aonuniformity factor through the 
generator. For the "5 sec" walls conditions of the Run 006-014 
simulation, the average conductivity nonuniformity factor is 
approximately 2. This means that any given station in the HPDE 
is operating with half the power density that would otherwise 
be obtained if the HPOE could provide perfect conductivity 
uniformity. Of course, doubling the local inte**action at every 
station in the UHO generator would profoundly alter the 
flowfield in the generator, so chat one cannot state 
categorically that the cold walls cause a 50% reduction in 
powei input. 
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As the wall temperature increases, the average plasma 
nonuniformity factor declines to a value of 1.60 at 1250 K, 
1.38 at 1500 K. and 1.22 at 1750 K. These reductions in the 
plasma nonuniformity factor correspond to potential increases 
in the MHO interaction in the generator of 27%, 47%, and 66%, 
respectively, over the "5 sec" wall temperature distribution. 

The variation of the conductivity nonuniformity 
factor, which is an intrinsic characteristic of the channel 
flow, is shown in Pig. 4-5. As the wall temperature increases 
from the "5 sec" values to 1250, 1500, and 1750 K, the average 
conductivity nor.uniformity factor, g, declines from 1.26 to 
1.17, to 1.11, and to 1.06, respectively. As was shown in 
(4-2], values of this parameter implicit in the original design 
of the HPDE were typically on the order of 1.03, corresponding 
to wall temperatures between 1900 K and 2100 K. The "5 sec" 
wall temperature distribution is much more characteristic of 
small scale generators such as the Avco Mark VI and the UTSI 
channel . 

4.2.2 Dependence on Magnetic Field 

The dependence of the total electrode voltage drop upon 
the magnetic field is explained by two principal competing 
mechanisms. First, because of the increase in Hall parameter, 
the effects of conductivity nonuniformity are magnified by the 
contribution of the Hall parameter term in the plasma 
nonuniformity factor. Compensating this effect is t.ie fact 
that an increased magnetic field creates increased current 
density in the electrode boundary layers. This tends to heat 
the boundary layers through Joule dissipation and improve 
somewhat the overall uniformity of the conductivity in the 
cross-section . 
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In the neighborhood of the HPDE Run 006-014 operating 
conditions, these two competing effects tend to slightly 
increase the electrode voltage drops as the magnetic field is 
increased to 4 T. Pig. 4-6 illustrates the almost linear 
dependence of the electrode voltage drop at four stations in 
the MHD generator as the peak magnetic field is increased from 
the Run 006-014 value of 3.27 T to 4.0 T and to 4.25 T. 

Relative to the total induced voltage, however, the 
electrode voltage drop does not increase very dramatically as 
the magnetic field increases. Fig. 4-7 illustrates the near 
insensitivity of the relative voltage drop to the peak magnetic 
field value between 3.27 T and 4.25 T. This result suggests 
that one-dimensional models which rely upon a reduction in the 
induced field in order to represent the electrode voltage drops 
may be adequate in the neighborhood of the HPDE Run 006-014 
conditions as long as the axial variation of the reduction 
factor is accounted for. 

The fact that the relative voltage drop is weakly 
dependent on peak magnetic field is related to the fact that 
the plasma nonuniformity factor is essentially constant with 
magnetic field, as is Illustrated in Fig. 4-8. Again, the 
competition between the declining conductivity nonuniforroity 
factor due to Joule dissipation (shown in Fig. 4-9) and the 
increasing Hall parameter due to the increase in the peak 
magnetic field gives a nearly constant overall plasma 
nonuniformity factor variation. 

The speculation that at high magnetic fields the voltage 
drops would be reduced by increased Joule dissipation in the 
boundary layers has not been borne out, at least in the 
neighborhood of the Run 006-014 operating conditions. The Hall 
parameter increase is an equal or even dominant effect over the 
g factor reduction due to Joule dissipation. 
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4.2.3 Dependence on Load Factor 


The variation of the electrode voltage drops and the 
plasma uniformity with electrical load factor provides a vivid 
illustration of the influence of Joule dissipation in the 
boundary layers upon the performance of IIHD generators. As the 
load factor increases and the generator goes more toward open 
circuit, the Joule dissipation in the boundary layer decreases. 
As the boundary layt.^s become cooler, the plasma nonuniformity 
factor, G, and conductivity nonuniformity factor, g, increase 
substantially. In addition, with the reduction of current 
density associated with increasing load factor, less electron 
nonequilibrium can be sustained at the electrodes, so that the 
near electrode region becomes substantially more resistive. 

The effect of increasing the Faraday load factor on 
the conductivity nonuniformity factor is shown in Fig. 4-10. 
At a load factor K * 0.55, the average conductivity 
nonuniformity factor g = 1.22. As the load factor increases to 
0.60, g » 1.24. At the Run 006-014 simulation conditions, in 
which the average load factor is approximated 0.61, the 
average value of g increases to 1.26. When the . ..ad factor is 
held constant at 0.70, the magnitude of g increases to 1.34. 
This is an important effect for cold wall generator designs 
which depend on Joule dissipation to alleviate some of the 
thermal gradients in the electrode wall boundary layers. 

The effect of increasing load factor on the overall 
plasma nonuniformity factor, G, is Illustrated in Fig. 4-11. 
The increase in G is even more dramatic than the increase in g 
because, as the channel goes toward open circuit, the flow 
velocity increases and the static pressure in the generator 
generally decreases. Hence, as the channel is unloaded, not 
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only is there less Joule dissipation in the boundary layers, 
but the Hall parameter also tends to increase. This results in 
a steeply increasing plasma nonuniformity factor. The 
sensitivity of G to load factor reemphasizes the need for 
calculational procedures which can account for the Joule 
dissipation in the boundary layers. A calculation which 
assumed that the value of g and/or G was independent of B or K 
would clearly be incorrect for this high interaction device. 

The behavior of the total electrode voltage drops as a 
function of load factor reveals a very interesting phenomenon 
associated with the accumulation of Joule heat through the HHD 
generator. Fig. 4-12 illustrates the values of the electrode 
voltage drop at four stations in the AEDC/HPDE. Four 
calculations are presented. The Run 006-014 simulation, 
CHPQUUOJBO, and three parametric variations with constant load 
factor throughout the generator are shown. Also shown are the 
data from Run 006-014, assuming that the load factor at each 
measurement station is as it is computed by the simulation. 

The computed dependence of the electrode voltage drop 
on load factor is rather weak for the cases in which the load 
factor is held constant. (Of course, a constant voltage drop 
with increasing current density at lower load factors implies 
increased power losses.) 

Ou the other hand, the simulation of Run 006-014, 
which incorporated a variable load factor (approximately 0.69 
at le inlet to approximately 0.55 at the exit) deviates 
i..creasingly from the constant load factor trend lines toward 
the back end of the generator. The maximum voltage drop in the 
Run 006-014 simulation is approximately 2200 volts, whereas the 
maximum electrode voltage drop with constant load factors 
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occurs at a load factor of approsimately 0.60, and has a value 
of approximately 1800 volts. The implication of this result is 
that the declining load factor profile throughout the generator 
in the Run 006-014 conditions cause less Joule dissipation 
upstream in the generator than would occur if the load factor 
were held constant at the avt:rage Run 006-014 values. This 
leads to generally colder electrode boundary layers in the 
downstream regions. This result suggests that the optimum load 
factor profile from the point of view of minimizing voltage 
drops would be one which begins initially lower than the 
average value and increases throughout the generator. While 
this is considered to be a second order effect on the overall 
performance of the device, it certainly bears further 
investigation for the future loading of the AEDC/HPDE. 

4 . 3 Parametric Variations Around the Transient Simulation 
of Run 006-014 


The response of the HPDE to parametric variations in the 
neighborhood of the Run 006-014 operating parameters was 
analyzed with the Q3D code family in the previous section. 
That study identified increases in the wall temperature and 
magnetic field as the principal means for achieving increased 
enthalpy extraction and isentropic efficiency of the generator. 
It is desirable to extend these results to the ultimate 
capabilities of the facility, i.e. , 6 T operation. This 

section presents such results. 

In the neighborhood of the Run 006-014 simulation 
studied with Q3D, the flow in the HHD generator is uniformly 
supersonic. Par from this neighborhood (i.e., at much higher 
magnetic fields or wall temperatures) the flow may become 
transonic or subsonic. Because of the additional expense and 
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effort required for Q3D computations in the transonic or 
subsonic modes, which require iterative matching of downstream 
pressure conditions, it was decided to extend the parametric 
studies far from the neighborhood of the Run 006-014 conditions 
by using the quasi-steady option of the TRANSIENT family of 
codes . 


To make such a parametric survey meaningful , it was 
necessary to incorporate as much of the information from the 
quasi-three-dimensional calculations using the Q30 family of 
codes into the TRANSIENT calculations. In order to include the 
effects of variable wall temperature on the electrical solution 
in the quasi-steady TRANSIENT calculation, values of the 
conductivity nonuniformity factor, g, obtained from Q3D 
simulations, were used as input. For the parametric survey of 
wall temperature, the g values from Q3D Run CHPQUUUJBO were 
used for the "cold wall" or "5-sec" wall electrode wall 
conditions. Q3D calculation COBQFAVJDY was used to obtain a 
"hot wall" or constant 1500 K electrode wall temperature 
distribution of g values. 

The starting point for all of the TRANSIENT parametric 
studies was Run BRCYTRCJEB described in Section 3.8. Utilizing 
the g factors corresponding to cold walls, the magnetic field 
distribution was increased progressively in steps from the 
nominal value of 3.27 T to 4 T, 5 T and then to 6 T. At each 
step in the magnetic field, sufficient time was provided to 
obtain a steady-state solution before the next increase in 
magnetic field was applied. Similarly, the hot wall simulation 
was carried out by first instantaneously modifying the g factor 
distribution to correspond to 1500 K walls of Q3D Run 
COBQFAVJDY, and then increasing the peak magnetic field 
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strength in steps from 3.27 T to 4 T, 5 T and then to 6 T. The 
results of this parametric variation in wall temperature and 
magnetic field are shown in Fig. 4-13. Throughout all of these 
parametric variations, the load schedule was fixed according to 
the Run 006-014 distributions, and the g distributions were 
held fixed at either the "hot wall" or "cold-wall" values. In 
view of the dependence of g on magnetic flux density discussed 
in the above sections, the consequence of holding g fixed is 
probably a slightly conservative estimate of performance. 

As is clearly evident in Fig. 4-13, the reduction of the 
plasma nonuniformity factor due to the increase in wall 
temperature from the "5 sec" values to a constant 1500 K 
distribution makes the difference between achieving the 
enthalpy extraction goal of the HPDE and not achieving it. The 
6 T operation of the HPDE with 1500 K walls is predicted to 
achieve nearly 15.5% enthalpy extraction with the Run 006-014 
nominal load schedule. With "5 sec" walls the HPDE achieves 
less than 13% enthalpy extraction. A detailed list of the 
isentropic efficiency, enthalpy extraction ratio, and the exit 
Mach numbers are shown in Table 4-1. 

It is seen from Table 4-1 that between 4 and 5 T with 
hot walls, there is a transition from supersonic to transonic 
flow with a shock in the generator. While the performance is 
predicted to increase despite the entry of the shock into the 
generator, the ramifications of operating a high-interaction 
MHD generator with a normal shock in the active portion of the 
generator remain to be investigated. 
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4.4 


Optimal Load Schedules 


The goal of the HPDE is performance, not duration. 
Therefore, the optimization of the loading of the HPDE should 
be aimed at optimizing performance without the constraints of 
"rule of thumb" or conventional limitations on current density. 
Hall field, or other parameters. 

For given conditions of velocity, conductivity, and 
magnetic field, the most important parameter in a Faraday 
generator is the electrical loading factor, K » <Ey>/<UB>, 
which defines the operating point on the load line of a local 
section of the generator (K = 1 is open circuit, K = 0 is short 
circuit). For given conditions of conductivity, velocity, and 
magnetic field, the maximum enthalpy extraction per unit length 
is achieved at a load factor of K * 0.5. For maximum 
efficiency, the load factor is somewhat hig ler than 0.5. 

In general, for generators of moderate and high 
interaction, the load-line is not linear. Changing the loading 
of the generator profoundly alters the distributions of 
conductivity and velocity in the generator, and 
multidimensional phenomena enhance the nonlinearity of the 
electrical performance. Therefore, in order to survey the 
performance of the device as a function of loading, the model 
employed must be capable of accurately predicting the 
interaction between the MHO forces and the fluid. For the 
present studies of the HPDE loading, a code from the TRANSIENT 
code family is used in order to facilitate rapid surveys of the 
Impact of the loading of the generator. The nonlinear coupling 
between the loading and the fluid behavior is described in this 
section. 
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While sophisticated, constrained optimization procedures 
could be utilized for optimizing the loading schedule of the 
HPOE. it is felt that the variation of the load factor in the 
axial direction produce second order effects on the 
performance. The most important parameter is the average load 
factor. This conclusion is borne out by comparison with 
constant load factor calculations and variable load factor 
calculations, as was discussed in Section 4.2. Therefore, all 
of the load survey studies that have been conducted to date 
have been carried out with the specification of a constant load 
factor throughout the generator. 

Many results have been presented in separate reports for 
the HPDE "nominal" operating conditions. Perhaps more 
pertinent are the results obta.,ned when parametric variations 
of the load factor around the simulation of an actual high- 
powered run of the HPDE are calculated. Section 4.1 
illustrates the load factor influence around the operating 
conditions simulated by Q3D for the Run 006-014. It was shown 
in that section that the reduction of the load factor in the 
Run 006-014 experiment could have improved the enthalpy 
extraction somewhat, but at a penalty in isentropic efficiency. 
In Section 4.3, the performance at various magnetic fields up 
to 6 T was compared with the TRANSIENT code with hot and cold 
wall assumptions. These studies utilized the load resistance 
schedule of Run 006-014. 

Relative to the achievement of the performance goals of 
the HPOE, it was suggested that the variation in the load 
factor, particularly at the "hot wall", high magnetic field 
points in Section 4.3, might allow both the enthalpy extraction 
and isentropic efficiency goals to be achieved simultaneously. 
Consequently, the electrical load factor was varied around each 
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of the 4 T, 5 T and 6 T points with hot walls described in 
Section 4.3. The results of these load factor variations are 
shown in Figs. 4-14, 4-15, and 4-16 for 4 T, 5 T, and 6 T, 
respectively. Also shown on Figs. 4-14 through 4-16 are the 
points computed with the actual Run 006-014 load schedule. In 
each case the Run 006-014 load schedule yields slightly better 
performance than the maximum obtained with constant load 
factors. However, it is felt that there is not enough of an 
improvement to warrant a global optimization study with the 
quasi-one-dlmenslonal TRANSIENT code. It would be more 
appropriate to carry out an optimization with a 
multidimensional code which could independently account for the 
boundary layer dissipation effects described in Section 4.3.3 

The importance of operating the UHD generator near the 
optimiim Mach number profile down the generator is emphasized by 
Figs. 4-14, 4-15 and 4-16. The flow regimes of the generator 
at the various load factors are indicated by the density of the 
trend curves. It is apparent that the enthalpy extraction of 
the generator optimizes in each case for 4, 5, and 6 T at 
transonic flow conditions. The optimum Mach number for maximum 
enthalpy extraction rate appears to be that load factor which 
provides transonic flow. At 4 T, the optimum load factor is 
approximately 0.65; at 5 T, it is approximately 0.72; and at 6 
T, it is approximately 0.80. The meaning of this is that, as 
the magnetic field increases, the load factor must also 
increase (i.e., vhe generator must be less and less heavily 
loaded) in order to keep the flow field distribution close to 
the optimum design values. This concept of maintaining the 
flow field at or near the optimum Mach number distribution is 
also central to the understanding of the part load and off- 
design operation of generators of all scales, as is discussed 
in Section 7.2. 
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4.5 Near Electrode Phenomena with Finite Electrode 
Segmentation 

In order to establish a basis for the Q3D current 
tra»iCport models, near electrode phenomena were studied under 
the operating conditions of the AEDC/HPDE Run 006-014. Two 
principal code familes were utilized for the study of near 
electrode effects with finite electrode segmentation. The FIN 
family of codes computes two-dimensional distributions of 
electrical and plasma transport variables in the region between 
a single pair of electrodes. The assumption implicit in this 
code family is that the variation in behavior from electrode to 
electrode is sufficiently weak that periodic boundary 
conditions can be imposed on the upstream and downstream ends 
of the region being computed. On the other hand, the ARRAY 
family of codes considers a string of electrodes without the 
requirement of periodicity. The upstream and downstream 
boundary conditions are obtained from Q3D or from the 
imposition of Neumann boundary conditions specifying local ..y 
invariant electrical properties. 

Two FIN calculations were performed at the x *» 2.77 m 
station in the AEDC/HPDE. Input data for gasdynamic profiles 
and initial conditions were obtained from the Q30 Run 
CHPQUUUJBO, the nominal simulation run for the HPOE Run 
006-014. Two computations were carried out. Computation 
BRCYICLJEA utilized the option with finite reaction rates. 
Computation BRCYLARJD4 utilized the option for instantaneous 
reaction rates for electron reactions. Table 4-2 shows a 
comparison between the Faraday voltage, total voltage drops. 
Hall field, and local power density, computed by the two FIN 
calculations and the Q3D simulation. We note that the 
inclusion of finite electrode effects in the FIN calulations 
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results in the reduction of the Faraday voltage by 68 volts. 
This Far’^lay voltage reduction is associated with an increase 
in the total electrode voltage drop of 128 volts. The Hall 
voli in the FIN calculations is approximately 50 V/m higher 
in mi ^ de, and the overall power density at this station is 
slightly decreased (4.5%) by the inclusion of the finite 
electrode effects. The relative small differences between the 
Q3D simulation CHPQUUUJBO and the FIN calculations indicates 
that the diffuse discharge model is suitable for uiations 
under the Run 006-014 conditions. 

On the other hand it is possible that the inclusion of 
finite electrode effects and near electrode effects may explain 
differences between the Run 006-014 measurements and the Q3D 
computation CHPQUUUJPO which were discussed in Section 3.9. It 
is anticipated that, as the boundary layers grow and the bulk 
nonuniformity of the boundary layers begins to dominate the 
near-electrode effects on the voltage drops, the importance of 
near-electrode effects in the back end of the generator will be 
less. Insufficient resources were available to investigate the 
near-electrode behavior throughout the entire MHD generator 
under the present effort. It would be of interest to complete 
thi analysis throughout the entire length of the generator 
under the Run 006-014 simulation conditions. 

Two ARRAY calculations were carried out to investigate 
near-electrode phenomena without the requirement of 
periodicity. Computation BRCYYMVJEC also considered the 
nominal conditions of operation in the neighborhood of x - 2.77 
m in the AEDC/H>»PE Run 006-014. A region of 2.5 channel 
diameters around x « 2.77 m was chosen for simulation. An 
expanded view of the current pattern in the near-electrode 
regions computed by the ARRAY family of codes at x * 2.77 m is 
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shown in Figs. 4-17 and 4-18 for the cathode region and anode 
region. Instabilities in the near-anode region is apparent in 
Fig. 4-17. The operating characteristics of the electrodes at 
X = 2.77 m is quite comparable to the FIN’ results discussed 
above and with the measurements of Run 006-014. 

Fig. 4-19 illi’strates the current pattern in the entire 
computational region of computPtion BRCYYMVJEC. In this 
figure, the details of the near-electrode behavior is 
suppressed by the coarseness of the plotting grid. 

Fig. 4-20 shows the corresponding equipotentia ■ s in tnis 
ARRAY computation. It is noted that the equipotentials are 
rather uniform. The steep inclination of the equipotential in 
the core in due to the relatively low Hall parameter (1.2) at 
this station. The curvature of the equipotentials near the 
walls reflects the presence of voltage drops due to the imposed 
gasdynamic boundary layer profiles. 

The importance of relaring the periodicity assumption is 
emphasized in Figs. 4-21 and 4-22. These Figures illustrate 
the result of calculation BRCYYQPJEC, in which the exact 
conditions of th previous ARRAY calculation were varied only 
by assuming that a single electrode pair went to short-circuit. 
The imposition of a short-circuit on this region causes a gross 
distortion of the equipotential pattern of Fig. 4-20, as shown 
in Fig. 4-21. Walls which normally operate at a potential 
difference of 1700 V or more are suddenly forced, in a very 
small region, to operate at the same potential. Very consider- 
able electric fields developed due to this condition are likely 
to cause severe electrical breakdowns in the neighborhood of 
the Faraday short. 
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In addition, all of the power being produced in a region 
approximately one or two diameters upstream and downstream of 
the MHD generator is available to be coupled into the Faraday 
short. As a result, the surge current through the Faraday 
short increases by two orders of magnitude over the normal 
operating current of the electrode pair. It is evident that 
the increase in the current to the external load circuit would 
be enough to cause either mechanical or thermal damage to the 
load circuit connections. This and similar faults must be 
protected against in tuture HPDE operations. 

The result of this "Faraday catastrophe" vividly 
illustrate^ that fine electrode segmentation is both a blessing 
and a curse. Fine electrode segmentation increases the overall 
current uniformity and decreases Hall effect reductions in the 
power output of the generator. On the other hand, if Faraday 
faults do occur, they will occur with more damaging 
consequences with finer electrode 

4-6 Current Distribution in the End Begions 

The HPOE inlet eddy currents were investigated, 
utilizing a two-dimensional, finite electrode, finite reaction 
rate code from the INLET family of codes. The operating 
conditions were specified from the outputs of the Run 006-014 
simulation by Q3D (CHPQUUUJBO) . The presence of open circuit 
electrodes in the magnetic field gradient of the end regions of 
the HPDE leads to internal circulating currents in the plasma. 
These are depicted in Fig. 4-23 for the computation BRCYZGGJEC. 
Because the passage of current in the downstream end of the 
current loops is through a higher magnetic field than the 
return current through the upstream portion of the loops, there 
is a net reduction in the total pressure of the power train. 
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This loss of pressure does not produce useful work. Under the 
Run 006-014 operating conditions, it is estimated, based on 
this computation, that this pressure drop is of the order of 
0.07 atmospheres. 

Fig. 4-24 shows the equi potentials associated with the 
current pattern in Fig. 4-23. It is interesting to note that 
the lack of current through the upstream boundary layers make 
these boundary layers very good insulators. The cold upstream 
boundary layers effectively shield the core from the 
open-circuit electrodes, which serve as insulators, and promote 
the circulation of current inside the plasma. It is noted that 
the voltage in the core at the first loaded electrode persists 
for a considerable distance upstream of the generator. 
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TABLE 4-1 


Results of TRANSIENT Computations Performed for Task I 
Nominal Loading from HPDE Run 006-014 
and Combustor Thermal Input of 303.2 HW 


Wall 


Run 

max 

^T) 

M . ^ 
exit 

Temp 

<K) 

BRCYTRCJEB 

3.274 

1.794 

"5 s 


4.0 

1.648 

•*5 s 


5.0 

0.651 

"5 s 


6.0 

0.654 

"5 s 

BRCYUXWJEB 

4.0 

1.124 

1500 


5.0 

0.650 

1500 


6.0 

0.699 

1500 


0-4361 


Power 

(MW) 

Enthalpy 

Extraction 

(%) 

Efficiency 
(1) . 

23.60 

8.24 

37.46 

28.95 

10.11 

40.72 

34.17 

11.94 

42.25 

36.39 

12.71 

45.11 

36.80 

12.83 

43.88 

41.75 

14.56 

51.56 

44.39 

15.48 

53.81 


199 



TABLE 4-2 


Electrical Performance Parameters at x * 2.77 m 
in the AEDC/HPDE Run 006-014: Comparison Between 
FIN Calculation BRC7ICLJEA and Q3D Calculation CHPQUDUJBC 


Quantity 

FIN 


Run No. 

BRCYICLJEA 

CHPQUUUJBO 

Total Electrode 

Voltage Drop, V 

830 

702 

Faraday Voltage, V 

1676 

1743 

Centerline Hall 

Field, V/m 

984 

959 

Local Power , 

Density, WH/ur 

16.05 

16.78 


0-4362 
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RcSEARCH COMPOWATtON 



5 6 7 8 9 10 11 

0-4125 ENTHALPY EXTRACTION RATIO, % 

Fig. 4-1. Generator isentropic efficiency vs. enthalpy extrac- 
tion ratio for simulation of HPDE Run 006-014 and 
simulations with variations in the conditions from 
those of Computation CHPQUUUJBO. Three-letter code 
by each point is the abbreviated STD Computation num- 
ber. 201 
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4-2. Coordinate systems and stations at which total voltage drops are measured 
and computed voltage drops are tabulated. 




TOTAL ELECTRODE VOLTAGE DROP. 








AVERAGE PLASMA NONUNIFORHITY FACTOR 


RESEARCH corporation 


3.00 


|e» 


2.00 h 


1.00 *- 
500 

0-4170 

Fig. 


AEDC/HPDE 


X RUN 006-014 SinULATION 
O PARAHETRIC VARIATION OF 
ELECTRODE SURFACE TEMPERATURE 



^ i I 

1000 1500 _ 2000 

AVERAGE ELECTRODE WALL TEMPERATURE, (K) 

4-4. Variation of the plasma nonuni forraity factor G with 
electrode wall temperature: Parametric variations 
around the AEDC/HPDE Run 006-014 simulation conditions. 
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AVERAGE CONDUCTIVITY NONUNIFORMITY FACTOR 


RCSEARCM CORPORATION 


AEDC/HPDE 


X RUN 006-014 SIMULATION 
O PARAMETRIC VARIATION OF 
ELECTRODE WALL TEMPERATURE 




1.10 



1.00 


1000 


1500 


0-4171 AVERAGE ELECTRODE WALL TEMPERATURE. T^ (K) 

Fig. 4-5. Variation of the conductivity nonuniformity factor 
with electrode wall temperature: Parametric varia- 
tions around the AEDC/HPDE Run 006-014 simulation 
conditions. 




RESEARCH CORPORaTsOM 


AEDC/HPDE 



0-4120 PEAK MAGNETIC FIELD, BIT) 


Fig. 4-6. Variation of the total electrode voltage drop 

with peak magnetic field at various stations in 
the AEDC/HPDE; parametric variations around the 
Q3D simulation of the HPDE Run 006-014 
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RELATIVE TOTAL ELECTRODE VOLTAGE DROP. Vh/UBD 


fit St ARCH CORPORA! iON 


AEDC/HPDE 



3.0 A.O 5.0 6.0 


0-4123 PtAK MAGNETIC FIELD. B (T) 


Fig. 4-7. Variation of the relative electrode voltage drops 
with peak magnetic field at various stations in 
the AEDC/HPDE: parametric variations around the 
Q3D simulation of the HPDE Run 006-014 
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AVERAGE PLASMA NONUNIFORMITY FACTOR 6 » g + (g-1) /3 


ReSCARCM corporation 


AEDC/HPDE 



3.0 ^.0 5.0 6.0 

0-4122 PEAK MAGNETIC FIELD. B (T) 


Fig. 4-8. Variation of the plasma nonuniformity factor with 
peak magnetic field in the AEDC/HPDE: parametric 
variations around the Q3D simulation of the HPDE 
Run 006-014 
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CONDUCTIVITY NONUNIFORMITY FACTOR, g =<or><\/o> 


RESEARCH COP'»ORaT‘ON 


AEDC/HPDE 



3.0 ^,0 5.0 6.0 

0-4121 PEAK MAGNETIC FIELD, B(T) 


Fig. 4-9. Variation of the conductivity nonuniformity 

factor with peak magnetic field in the AEDC/HPDE: 
parametric variations around the Q3D simulation 
of the HPDE Run 006-014 


2)3 



CONDUCTIVITY NONUNIFORMITY FACTOR, g ^<o><\/(r> 


«ESt*»CN ON 


AEDC/HPDE 



Fig. 4- 10. Variation of the conductivity nonuni formity 
factor with Faraday load factor; parametric 
variations around the Q3D simulation of the 
AEDC/HPDE Run 006-014. 
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PLASMA NONUNIFORMITY FACTOR. G = g + (g-l)/J 


RCSEAACH COBPOPATiON 


AEDC/HPDE 



Fig. 4-11. Variation of the plasma nonuniformity factor with 
Faraday load factor: parametric variations around 

the Q3D simulation of the AEDC/HPDE Run 006-014. 
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RESEAPCH COPPOPAIiON 


AEDC/HPDE 



X = 4.26 M 
(STA 215) 


•/ = 2.78 M 
(STA 77.4) 

X = 1.43 M 
(STA 24.7) 


0-4119 FARADAY LOAD FACTOR, K = <Ey> /<UB> 


Fig. 4-12. Variation of total electrode voltage drop at various 
stations in the AEDC/HPDE with Faraday load factor: 
parametric variations around the Q3D simulation of 
Run 006-014. 


212 



Ml'jCAnCM COnPOHATiCN 



Fig, 4-13. Generator isentropic efficiency vs, enthalpy extraction ratio 
for simulation of HPDE Run 006-014 and simulations with 
variations in the magnetic field and wall temperature distribution 
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BtSCABCH ConPOOA’iOM 


PEAK MAGNETIC FIELD; 4 TESLA 



Fig. 4-14. Generator isentropic efficiency vs. enthalpy extraction 
ratio for simulation of AEDC/HPDE Run 006-014 with 
variations in load factor and a magnetic field of 4 tesla 
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BfS€»oC«* COBPOB*’ iiN 


PEAK MAGNETIC FIELD; 5 TESLA 



ENTHALPY EXTRACTION RATIO (%) 


Fig. 4-15. Generator isentropic efficiency vs. enthalpy extraction 
ratio for simulation of .4EDC/HPDE Run 006-014 with 
variations in load factor and a magnetic field of 5 tesla 
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ncSCAOCH CORPORATION 


PEAK MAGNETIC FIELD; 6 TESLA 



ENTHALPY EXTRACTION RATIO (%) 


Fig. 4-16 . Generator isentropic efficiency vs. enthalpy extraction 
ratio for simulation of AEDC/HPDE Run 006-014 with 
variations in load factor and a magnetic field of 6 tesla 
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CONTOUR INTERVAL: 50 V 



5.0 


CRITICAL PHENOMENA FOR THE U.S. U-25 EXPERIMENT 


An MHD generator for installation in the Soviet U-25 
facility was designed and fabricated in the United States as 
part of the Joint U.S./U.S.S.R. cooperative MHD Program. STD 
has simulated the "as built" performance of the U.S. U-25 
channel to determine the nominal operating performance and to 
analyze the critical phenomena of this experiment. 

5 . 1 Channel Characterization: U-25 Nominal Operating 
Conditions 


The values of the physical and nominal operating 
parameters used for simulations of the U.S. U-25 channel are 
discussed in this section. 

5.1.1 Geometry 

The major components of the flow train sections of the 
U.S. U-25 experiment are the KS-3 combustor (Soviet design and 
fabrication), the U.S. U-25 nozzle and channel assembly, and 
the U.S. U-25 diffuser assembly. The KS-3 combustor is a 
cylindrical unit with an inside diameter of 1.7 m [5-1). It is 
connected to the main oxidant feed line, which is of 1.5 m 
inside diameter, by a short conical diffuser section. The 
combustion zone lies entirely within the conical divider, which 
is perforated to meter the oxidant flow to the upstream and 
downstream parts of the combustion zone. The combustion zone 
and the first exit transition piece is of circular section (1.2 
m in diameter) at the inlet and of rectangular section (640 mm 
by 760 mm) at the exit. The nozzle of the U.S. U-25 flow train 
assembly is connected directly to the flange at the exit of the 
first transition piece. The nozzle inlet plan section is 624 
mm by 744 mm (5-2). The nozzle contours are defined by the 
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753.7 mm tangent radius on the sidewall contour and a straight 
taper from 744 mm at inlet to 730 mm at nozzle exit on the 
electrode wall contour. The exit of the nozzle mates direcilv 
with the inlet of the channel. 

The channel frame number, width, height and 
diagonallzation angle are given as functions of axial distance 
in Table 5-1 15-3J. The origin of the axial distance 
coordinate is the 0-0 station, which corresponds to the front 
vertical edge of the magnet pole. This coordinate runs through 
the center of the channel, midway between both the electrode 
walls and the sidewalls. The U.S. U-25 generator is 
rectangular along its entire length. 

5.1.2 Magnetic Field 

The axial distribution of the magnetic flux density is 
given in Table 5-2 15-4). 

5.1.3 Loading 

The U.S. U-25 experiment is planned to be run with a 
single load at a nominal load current of 2840 A (5-5). The 
loading is connected with distributed current collection 
through diodes (International Rectifier 501VFR200) on frames 1 
through 67 on the front end and diodes (International Rectifier 
501VFR200) on frames 366 through 432 in the back end (5-3j. 

5.1.4 Working Fluid 

The working fluid for the U.S. U-25 experiment is 
defined in this section. The mass flow rate for the nominal 
conditions is 50 kg/s [5-31. 
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The composition of the fuel is given by volume as 
[5-41, [5-6): 

CH^ 92.4 - 94.5% 

CgHg 3.33 - 2.11% 

C’3Hg 1.38 - 0.87% 

0.36 - 0.14% 

COg 0.70 - 0.70% 

Ng 1.5 - 1.5% 

O.^ 0.3 - 0.2% 

The fuel has a lower heating value of 45.98 to 46.51 MJ/kg and 
is supplied to the combustor in the gaseous phase at ambient 
temperature [ 5-4 j . 

The oxidizer is composed of air enriched with oxygen 
so that the total oxygen in the mixture is 40% by weight of the 
combined weight of the air and added oxygen [5-5j. The oxygen 
for enrichment is assumed to be pure O 2 [5-5 J. The air is 
taken from the atmosphere and has a nominal moisture content of 
0.89% by volume (equivalent to 50% relative humidity at 10®C) 
[5-5J. The total oxidizer flow is preheated to 1200 ^C before 
it enters the combustor [5-3], [5-5|. 

The seed is potassium carbonate in purified water 
solution [5-5J. The concentration of K 2^®3 water is 1:1 

by weight [5-5). The flow rate of seed is such that K 2 CO 3 is 
2 % of the total mass flow rate [5-5). The seed solution is 
injected at ambient temperature [5-5). 

For the nominal conditions the combustion is 
stoichiometric [5-3), [5-5). 
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5.1.5 


Wall Conditions 


The combustor wall temperature is 1773 K for all 
surfaces {5-5]. The nozzle, channel and diffuser wall 
temperatures are 644 {5-5J, 644 to 672 (5-5|, and 514 K (5-3J, 
respectively, for each wall and for the entire length of each 
section. The wall roughness is characterized by an equivalent 
sand roughness height of 3 nun in the combustor and channel and 
smooth walls in the nozzle and diffuser. 

5.1.6 Channel Inlet Conditions 

The conditions at the inlet of the U.S. U-25 are 
determined by the total mass flow rate, the thermodynamics of 
the working fluid, the inlet dimensions, and the enthalpy flux. 
The enthalpy flux was determined by assuming complete 
combustion and applying a 4.0% heat loss [5-1] for the 
combustor . 

5.1.7 Other Conditions and Assumptions 

The design value of the diffuser recovery coefficient 
is 0.45. The pressure at the end of the diffuser is 2 to 3 kPa 
below atmospheric [5-5j. 
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5.2 Steady-State Performance for Noainal Design 


The U-25 steady-state perforraance for the nominal design 
was simulated with the TRANSIENT code family in the 
quasi-steady mode of operation. The initial computations were 
performed under the optimistic condition of an assumed conduc- 
tivity nonuniformity factor of 1.06. This value was taken from 
preliminary D-25 calculations with the Q3D family of codes 
performed previously to this contract and corresponds to 
arc-mode current transport with a critical field of 12 kV/m> 
Definitions of arc-mode current transport and critical eiectric 
field are given in Section 3.3. For the cold walls of the 
nominal operating conditions, the arc-mode current transport 
assumption is likely to yield the highest estimate of power 
output. 


Under the nominal conditions, transonic flow is 
predicted for the U.S. U-25 experiment by computation 
BRCYKDXJGJ. The axial distributions of the gasdynamic and 
plasma property variables for the U-25 under the nominal 
conditions are given in Pigs. 5-1 through 5-10. The variables 
displayed in these figures are grouped into the following 
classifications: (1) open and short circuit variables; (2) 
electrical loading, (3) plasma properties; (4) electrical field 
variables; (5) average gasdynamic variables; (6) electric 
potential; (7) channel geometry; (8) heat transfer parameters; 
(9) power extraction parameters; and (10) optimization 
variables. The origin of the axial coordinate in these figur<-‘s 
is the upstream end of the KS-3 combustor. For this axial 
coordinate system, station 0-0 is located at an axial distance 
of 2.789 m. the upstream end of the nozzle is located at an 
axial distance of 2.100 m, the end of the channel is located at 
an axial distance of 9.851 m, and the end of the diffuser is 
located ac an axial distance of 12.827 m. 
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Experience with simulations of the HPDE indicate a 
higher conductivity nonuniformity factor for cold wall chancels 
than was utilized in computation BRCYKOXJGJ. Computation 
BRCYEHZJGK used the n<XDinal conditions and, in addition, had an 
assumed conductivity nonuniformity factor of 1.24, which was 
established from simulations of the HPDE under cold wall 
conditions. The results of computation BRCYEHZJGK are 
presented in Figs. 5-11 through 5-20. The variables presented 
in Pigs. 5-11 through 5-20 correspond to the same variables 
presented in Figs. 5-1 through 5-10. The flow remains 
transonic and the power output decreases to 7.74 MW 
(BRCYEHZJGK) from 10.97 MW ( BRCYKDX JG J ) . 

By decreasing the load current to obtain maximum power 
output, some of the power difference can be recovered. 
Computation BRCYJCGJGK assumed a conductivity nonuniformity 
factor of 1.24 and also assumed a load current of 2128 A. The 
results of computation BRCYJCGJGK are presented in Figs. 5-21 
through 5-30. The flow remains transonic, but the power 
increases to 9.44 MW, which is only slightly below the goal of 
10 MW. 


The isentropic efficiency is shown in Fig. 5-31 as a 
function of the gross channel power for these three 
simulations. This figure illustrates that of these three 
simulations of the nominal operating conditions only 
computation BRCYKDXJGJ predicts an output power above the goal 
of the experiment, but that simulation BRCYJCGJGK, with higher 
conductivity nonuniformity factor and lower load current, 
predicts that the output power will be very close to the goal. 
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5.3 Quasi-Steady Performance as a Function of Mass Flow Rate 


A series of five preliminary performance calculations 
with the nominal conditions, except for mass flow rate, were 
carried out with the TRANSIENT code family in the quasi-steady 
mode of operation. The range of mass flow rates was taken to 
be from 35 to 55 kg/s. The performance of the U-25 over this 
mass flow rate range was evaluated with an assumed conductivity 
nonuniformity factor of 1.06. This value was taken from 
previous preliminary U-25 calculations with the Q3D family of 
codes and is typical of values obtained with the assumption of 
arc-mode current transport with a 12 kV/m breakdown criterion 
in the U-25. The arc-mode current transport assumption is 
likely to be the most favorable current transport mechanism for 
cold wall channels and is likely to yield the highest estimate 
of power output. 

The five computations in this series are: A5MYISTIKM, 50 
kg/s; A5MYNEPIKN, 45 kg/s; A5MYNGCIKN, 55 kg/s; A5MY0NTIKN, 40 
kg/s; and A5MYRAW1KN, 35 kg/s. 

The transonic flow predicted by these calculations for 
all mass flow rates above 40 kg/sec is illustrated by Fig. 
5-32. Fig. 5-32 shows that the flow chokes from 3.2 m to 6.6 m 
downstream of the combustor backplate (the origin of axial 
coc 'dinate for these STD computations) and recovers to 
atmospheric pressure through a normal shock in the diffuser. 
Because the power extraction and MHD interaction in these 
computations are near maximum, it is expected that more 
conservative assumptions for the wall current transport 
mechanism will lead to even more supersonic flow conditions. 

The power output computed in these preliminary 
calculations is presented in Fig. 5-33. It is significant that 
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these prelim, nary calculations predict that power output in 
excess of 10 MW might be achieved even at a mass flow rate of 
40 kg/s. It must be reemphasized, however, that these 
calculations are based on the assumption of the most favorable 
current transport mechanism for this cold wall channel. 

The isentropic efficiency is given as a function of mass 
flow rate in Fig. 5-34. The maximum value of isentropic 
efficiency calculated was nearly 40% from A5MYNEPIKN. As is 
shown in Fig. 5-32, the shock enters the channel for the 
nominal conditions and for mass flow rates less than 44 kg/s. 
The decrease in isentropic efficiency caused by the shock 
entering the active channel is seen in Fig. 5-34. 
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TABLE 5-1. Design Data for the U.S. U-25 MHD Generator 


FRANC 

AXIAL 

DISTANCE 

NIOTH 

HEIGHT 

OIAGONALIZATIOli 

ANGLE 

NUNBCR 

(•) 

(■) 

U) 

(radians) 

1 


A.SSSiiT 

0»999t0 

#*90190 

1 

••.•aaTs 

AelSSlS 

0»?9000 

•»10#|2 

% 

•f.ATSIo 

A.S95«« 

o^liooo 

0.A0AA9 

A 



t^79AA0 

O.AlltO 

% 

•«.0*77« 

OeSSm 

0I79AAO 

0.AI7AO 

A 

•t.0S40t 

t.SSTOT 

0,19AOO 

•*A9A11 

? 

•t.OIoAS 

AelSTM 

0«79000 

A.A9090 

• 


A.J5T95 

0«79A00 

#*49oOo 

f 

•.0A74S 

AeSSAlO 

OeTlOOO 

•.##900 

|0 

4.42134 

0.19AA* 

0«790AO 

0*AA070 

ll 

4.4SSS0 

0«S99S9 

0e79000 

0.A9017 


4.44424 

0el59A| 

0»I9000 

0.AO2S2 


4.04332 

A.SA610 

0^79000 

O.AOaOA 

1: 

4.47734 

•elftOAl 

0«79000 

0,47516 

4.44152 

OeSAtSS 

0»79000 

0.40145 

lA 

4^ 14344 

OeSAlBS 

OeTlvOO 

0.40772 

• T 

4.HO40 

AeSA^S9 

0eT9O0O 

0.40907 

1? 

4.I34U 

A.SA294 

0»T9000 

0.90020 

4.14447 

AeSASSO 

O.79O00 

0.90641 


4.14243 

••Se4vA 

Oe79AOO 

0*91290 

fl 

4.17724 


0*79000 

0.91076 

fi 

4.t4|7o 

»«S#S2S 

0*79000 

0*92400 


4.24421 

0»9ASA6 

0.79000 

0.99102 

r 

4.22077 

0«9A4»A7 

0*79000 

0.99712 


4.23334 

0.9AT11 

0*79000 

0.94910 

|A 

0.23043 

0a|A77S 

0*79000 

0.94025 

>T 

0.24474 

0eS«A90 

0*79000 

0.95520 


0.27453 

0*SA9<»6 

0*79000 

0,5612V 


4.24436 

A*S6977 

0,79‘^Om 

0.56720 


0.30028 

o*57oai 

0*79000 

0.57924 


0.32-I22 

0.9711? 

0*79000 

0.57910 

Y? 

4.33422 

OeS7lA$ 

0*79000 

0.5051U 


4.33424 


0*79000 

0.9O0OV 


4.34041 

Oe9792A 

0*79000 

0,50606 


4.34454 

0*S7AaS 

0*79000 

0.602/1 


0.30083 

A«|7A79 

0*79000 

0.60053 


4.41313 

0.S7557 

0*79000 

0,61493 


0.43o3o 

A*S7A5S 

0*79000 

0.62OII 

tA 

4.44303 

AeS77t9 

0*79000 

0.62506 

«o 

4.44142 

A.IT797 

0*79000 

0,69190 

• 1 

4.47448 

0«97AA9 

0*79000 

0.69720 

«2 

4,44240 

oesms 

0*79000 

0.64200 

•s 

4.34424 

0*9»090 

0*79000 

0.64069 

•A 

4.32444 

AeSAt37 

0*79000 

0.65027 

«s 

4.33044 

0.9A224 

0*79000 

0.65007 

•A 

4.33570 

0.SA91I 

0*79000 

0,66945 

•T 

4.57172 

0.9S996 

0*79000 

0.67101 

•• 

4.54773 

0»SAAA9 

0*79000 

0.67695 

#A 

4.44340 

A.9A574 

0*79000 

0.60206 


4.42403 

AelAAA? 

0*79000 

0.60794 


4,43427 

0»SA79o 

0*79000 

O.6O9OI 


4.45244 

o«seas9 

0*79000 

0.60044 

fi 

4.44404 

AelAAS* 

0*79000 

6.70906 

fA 

4.44540 

0.S90I7 

0*79000 

0.70025 


4.70244 

O.lAfO* 

0*79000 

0.71461 

fA 

4.71443 

0.9919$ 

0*79000 

0.71004 


4.73335 

0.9920$ 

0*79000 

0.72925 

?• 

4.73211 

0.99975 

0*79000 

0.79094 


4.744O0 

0. 9^905 

0*79000 

0.7950O 

*0 

4.7450c 

0.19$$5 

0*79000 

0.74104 


9-3336 
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TABLE 5-1. Design Data for the U.S. U-25 MHD Generator (cont) 

axial diagonal ization 


FRAME 

DISTANCE 

WIDTH 

HEIGHT 

ANGLE 

NUMBER 

(•) 

It) 

(■') 

(radians) 

it 


OeStnnS 

0.73000 

0.74629 

a 

O.NNpOI 

• eSmt 

0.73000 

9.75144 


•;ns7iv 

•eS^tat 

0.73000 

0.75699 

I* 

o^NSaa? 

OeWlT 

0.73009 

0.76173 


»«e7lS2 

0»4000t 

0.73000 

0.76695 


e,ea« 2 & 

0»40t<*0 

0.73000 

0.77193 


»,<0N77 

0.4019? 

0.73000 

9.77700 


»»«<4S9 

0.40204 

6.73000 

0.79203 


e.9«7os 

0. 403^0 

0.73049 

0,79497 

r 

O.wu 

0.4v3tt 

0.73192 

0.7994U 


0.975117 

0.40401 

0.73270 

0.79540 

72 

C.99JH 

0.40421 

0.73399 

0.7954U 

7? 

I.OOtttP 

0.4044? 

0.73906 

0.79940 

r 

t.OAOSj 

0.4046? 

0.73623 

0.79540 

7» 

t.oocoa 

9.404M 

0.73741 

9.79540 

7* 

1.0S720 

0.40SO) 

0.73959 

0.79540 

V 

1.07)59 

OoAOS?; 

0.73977 

0.79940 


1.00900 

0.40S45 

0.74094 

0.79540 


t.|092) 

0 .4036(2 

0.74212 

0.7054O 

no 

t. 12257 

0.40967 

0.74330 

0.7954V 

M 

l,^»•'92 

0.40606 

0.74449 

0. 79540 

9 I 

I.IS52P 

0.40629 

0.74565 

0.79540 


I.t7l*| 

0.40696 

0.74693 

0. 79540 

;• 

1.1079) 

0.40672 

0.74901 

0. 79540 


». 29929 

0.44694 

0.74919 

0.79540 


l.22>'6« 

0.40719 

0.79036 

0.79540 

hT 

1.2)990 

0.40736 

0.79194 

0.7954V 

ne 

I.25)J5 

»a.4m790 

0.19272 

0.79540 

nO 

l.29'*t,7 

O.OCTHO 

0,75390 

0.79540 

♦0 

l.28i)( 1 

'*.4oa»»? 

0.79909 

0.79540 

• i 

i.)0?)i> 

6.40024 

0.75625 

0.7954O 

f? 

l.)l«Tc 

0.4VO4O 

0.75743 

0.79540 


l.))505 

0. 40960 

0.79961 

0.79540 


I.)51S9 

4.40991 

0,75979 

9.79940 


I.)077j 

0.40913 

0.76096 

0.79540 


t.)09oe 

0.40935 

0.76214 

0.79540 


t.OOPOO 

0.40998 

0.76332 

o « 79540 

•• 

I.0U77 

0.40991 

0.76490 

0.79540 


1 .9JSJI 

0.41003 

0.76967 

0.79540 

too 

1.40946 

0.41026 

0.76695 

0.79540 

tol 

I.4O500 

0.41090 

0.76903 

0,79540 

lo2 

1.40214 

0.41073 

0.76921 

0. 79540 

to^ 

1.49049 

0.41097 

0.m39 

0.7954O 

IO» 

I.SI40) 

0.41120 

0.77156 

0.79540 

to* 

1.5)110 

0.41144 

0.77274 

0.79540 

lOO 

1.54752 

0.41168 

0,77392 

0.70540 

10^ 

l,5*}0o 

«.4I19? 

0,77509 

0.79540 

loO 

1,501.21 

0.41216 

0.77627 

0*79540 


1.5905) 

0.41240 

0.77749 

0. 79540 

IfO 

1.01290 

0.41264 

0.77463 

0.79540 

»p 

1.02924 

0.41299 

0.77990 

0.79540 

llO 

1.0455a 

0.41313 

0.79099 

0. 79540 

t\y 

1.0019) 

0.41337 

0.79216 

0.79540 


1.07027 

0.4136? 

0,79334 

9.79540 

»i» 

1(09402 

0.41396 

0.79491 

0.79540 

lit 

1.71090 

0.41411 

0.79969 

0. 79540 

itT 

1.727)1 

0.41436 

0.79697 

0.7954U 


1,74)0) 

0.41401 

0.78909 

0.7954O 

l|0 

1.75999 

0.41496 

0.79922 

0.79540 

too 

1.770)1 

0.41911 

0.79040 

0.7954O 
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TABLE 5-1. 

Design Data 

for the 

U.S. U-25 

MUD Generator (cont) 



AXIAl 



OIAGOHALIZATION 

^RANC 

DISTANCE 

HIOTH 

HEIGHT 

ANGLE 

NUMSEn 

(•) 

(a) 

(a) 

(radians) 


it 

i.mM 

a.41519 

6ft79t56 

6^70540 


u 

UtOfOJ 

a.4tt92 

0.79276 

9.7B540 


A’ 

l»AiAST 

a.4l5«7 

0.79S93 

t.76540 


M 

i.Min 

teBlBtS 

a.795ll 

a.?a»** 


(s 

I.ASOoa 

».4iasa 

0.79629 

9.7B540 


»• 

t.674«0 

6.4I666 

6.79747 

a.rasaa 


»T 

UB90TS 

4.41990 

0,79Ma 

a.mao 


r 

l,9pTa9 

0.41719 

a,79*aa 

B.7B540 



t««2SlS 

€•417^2 

0.60IOO 

a.Tas4« 


U 

i.*we 

9#4176B 

C.60216 

a.tasav 

l4l 

I.W612 

9e4l79S 

a.aoSis 

B.7B540 


i.m«T 

0.41621 

0.60453 

a.m«o 

1 

is 

i.90«ei 

0.41646 

0.60571 

a.rasqo 

t 

«• 

t.OOSta 

0.41674 

o.aaaaa 

0.7BS4O 

IV 

2«02lStf 

0.4I9UI 

O. 6 O 6 O 6 

0.7B540 

IV 

s.ojraq 

0.41926 

9.60924 

0.7B54O 

tv 

i.0S4|« 

0e4l95S 

s«aia42 

a.7as«« 

tv 

A.OTPSJ 

0.41962 

0.01160 

0.76540 

iV 

2«0B6B» 

0.42009 

0.61277 

6.76540 

tav 

2.1 0122 

0.420SO 

i..ati9S 

0.76S40 

lal 

2.t:*ia 

0.42064 

o.aiiil 

0.7B540 

fZ 

2. lisa J 

0.42091 

0.61631 

0.7B540 

Ul 

2.tS221 

0.42119 

0.61746 

0.7654U 

1«« 

t.l^Bao 

0.42146 

0.61666 

0.TBS4O 

us 

2.16494 

0.42174 

0.61964 

0.76540 

la* 

2.20I2« 

0.42202 

0.62102 

0.7654O 

t«7 

2.2IT«S 

0.42290 

0.62219 

0.76340 

tiA 

2.7119T 

0.42296 

0.62337 

0.76540 

IV 

2.25 ;• 12 

0.42267 

0.62453 

0.76540 

iV 

2 •26060 

0.42319 

0.62373 

0.76S4O 

lal 

2.2SSM 

0.42344 

0.62«»9 

0. 76340 

iV 

2.296V., 

0.42371 

0.62C06 

0. 763-0 

iV 

2.31bo« 

0.42394 

0.62920 

0,76540 

1< 

i« 

2.1S2«q 

0.42427 

0.63044 

0.76S40 

tV 

2.S40SC 

0*42436 

0.63161 

0.76540 

tv 

2.S6473 

0.42465 

0.63279 

0.76540 

IV 

2.1«»»7 

0.42513 

U. 63397 

0.76340 

iV 

2,mii 

0.42542 

0.63513 

0. 76540 

iV 


0.42571 

0.63632 

0,76540 

tap 

2.41f 10 

0.42600 

0.63750 

0.76540 

lal 

2.<i«6qi 

0.42630 

0.63666 

0.76540 


2.‘»*27« 

0.42659 

0.63986 

0.76540 

us 

2.«79ll 

0.42669 

0.64103 

0. 76540 

U« 

2.«9S*# 

0.42716 

0.64221 

0,76540 

iV 

2.5* l»2 

0.42746 

0.64339 

0. 76540 

iV 

2.52»17 

0.42776 

0.64457 

0. 76540 


2,5«flS| 

0.426tt8 

0.64574 

0. 76540 

UA 

2.s*»a» 

0.42636 

0.64692 

0,76540 


2.57720 

0.42669 

0.64610 

C, 76540 


2.59150 

0.42699 

0.64926 

0.7B54O 

tv 

2.699,9 

0.42929 

0.65045 

0.70540 

*z* 

2.92621 

0.42960 

0.65163 

0.7B540 

tv 

2«9«25a 

0.42991 

0.63261 

0.7B540 

tr 

2.95692 

0.43022 

0.65399 

0.TB54O 

I2« 

2.97526 

0.43033 

0.65316 

0. 76540 

tv 

2.99(61 

0.43064 

0.65634 

0.70540 

itT 

2.70795 

0.43115 

0.65752 

0.76540 

11 

• 

2.72010 

0.43147 

0.65670 

0. 76540 

ll 

rA 

2.70060 

0.43176 

0.65967 

0,76540 

lao 

2.75696 

0.43210 

0.66103 

0.76540 
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TABLE 3-1. Design Data for the U.S. U-25 MHD Generator (cont) 



Axm 



OiAGONAlIZATlO/. 

FRAHE 

OISTAKCC 

MIOTH 

HEIGHT 

ANGLE 

NUMBER 

(■) 

(o) 

(•) 

(radians) 

tut 

*;T7JM 

4.4S242 

^.96223 

0,705411 

tll2 


4.43274 

0»96)4i 

9.79540 



4.433*7 

0.964S9 

9,79540 

tM 


0.43330 

0.96576 

9.79540 

mS 

2.«SeT| 

0.43371 

0.96694 

9.79540 

IM 

2.«5S0% 

0.43403 

0.06312 

0.19540 

IHT 


0.43437 

0.96929 

0.79540 

UR 

t.aaTT« 

0.43470 

0.67047 

0.79540 

in* 

S.aaaoa 

0.435*1 

0.97165 

0.79540 

HO 

2.72043 

0.43534 

0.97293 

9.79540 

HI 

2.7J0T7 

0.43540 

0.97400 

0.79540 

H2 

2.7SSI1 

0.437*1 

0.97519 

0.79540 

Hi 

2.70740 

0.4*735 

0.976)0 

0.79540 

H* 

2.7450'* 

0«4S6^A 

0.977S4 

0. 7954*1 

hs 

3.07215 

0.43702 

0.97971 

0 . 795 * 

IH 

3.01040 

0.4I7S0 

0.97999 

0.79540 

hT 

3.03443 

0.43770 

U.99107 

0.79^40 

IfR 

3.05110 

0.430U4 

0.99225 

0,7t»40 

hO 

3.04T52 

0.43039 

0.99)4) 

0.T9540 


3.04307 

0.43973 

0,99460 

6,79540 

2o< 

3.10021 

0.43** >9 

0.99579 

• .79540 

Hi 

3.II0S« 

0.43942 

0^99696 

0.79540 

to9 

3.IS27.'* 

0.43977 

0.99914 

0,79540 


3.10724 

0.44012 

0.999)1 

0,79540 

2o9 

3.17550 

0.44047 

0.99049 

0,79540 

?06 

3.14103 

0.4409) 

0.99167 

0,79540 

2 oT 

3.17020 

0.44f 19 

0,90295 

0.79540 


S.2t4«>2 

0.44tS4 

^l,9«^402 

0. 7954 11 


S.23t9o 

t»41|90 

0,#*V52« 

0.79540 


5.24731 

i«.412?<3 

6,906)9 

0.79*ino 


3.27345 

O.4i2o? 

0,90756 

0.79540 


5.2040m 

0.44299 

C.9997) 

0.79540 

oh 

5.27434 

0.44)35 

0,99991 

0.79540 

0i« 

3.31240 

0.44)7t 

0,90109 

0.79540 


3.32043 

0.44409 

0.90227 

0. 79540 

0 ‘H 

3.34537 

0.44445 

0,90)44 

0.19540 


3.37172 

0.44492 

0,90462 

0.79540 

01 

5.57004 

0.445H 

0,9059o 

0.79540 

OlO 

5.37441 

0.44556 

0,90699 

0.79540 

o>o 

5.41C75 

0.44594 

U.9091S 

0,79540 

0>i 

3.427*0 

0.44631 

0,909)3 

0.79540 

o»o 

3. •14341 

0.44609 

0,91051 

0.79540 

OOS 

3.45070 

0.44707 

0,91169 

0.79540 

0>« 

3.47613 

0.44745 

0,91296 

0. 79540 


5.47247 

0.4479) 

0,91404 

0.79540 

OH 

3.50001 

0.44921 

0,91522 

0. 79540 

0>T 

3.52514 

0,44960 

0,91640 

0.74540 

00« 

3.541 5-i 

0.44999 

0,91757 

0,76540 

0»« 

3. 55705 

0.44939 

0.91975 

0.795^0 

OfO 

3.57410 

0.44979 

0, 91993 

0,79540 


3.57053 

0.45019 

0.92111 

0. 79540 

oh 

3.70700 

0.45050 

0,92229 

0,79540 

OtS 

3.72322 

0.45099 

0.92346 

0.79540 

Oi« 

3.73057 

0.45139 

0,92464 

0,7*5«*:) 


3.75501 

0.45179 

0.92592 

0,79540 


3.77225 

0.45219 

0,92699 

0,79540 


3,74040 

0.45269 

0,9291/ 

0.70540 

oh 

3.7*404 

0.45)'JO 

0,92915 

0,79540 

OtO 

3.72120 

0.45)42 

0,93053 

0,79540 

o«o 

3.737*3 

0.45)9) 

0,95170 

0.79540 
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TABLE 5-1. Design Data for the U.S. U-25 MHD Generator (cont) 


FRAME 

AX!«l 

OISUICE 

MIOTM 

HEIGHT 

OIAGOMAt IZATIC-II 
ANGLE 

MUN8ER 

(9) 

(•) 

(•) 

(*adians) 

»•! 

}.T9S«« 


••65266 


ui 

S.T7*« 


••61«pft 

•,7«S4* 


J.T9990 

f.mOM 

6.65564 

*.7*S40 



»e«$S99 

6.6H4I 

•.7«M« 




6.65759 

•,7«S4« 


l.tSSTo 


•.65677 

•.7*»4* 


>.M2C« 

•»4S4T^ 

6.65695 

•.7»940 


>•*•>»• 


6.661 16 

•.7*S*4 



•s4$T#| 

6.66216 

•.7»S49 


l.**l*7 


•.66546 

•,7»$44 

9%l 

».*|7*2 


••96666 

•,7«S4* 

Ai 

I.7JJ7* 

f«49H| 

•.66565 

•.7*544 




••66761 

*.7*54« 

P%^ 

!.•••*% 

P.^HTT 

6.666|9 

*.7*540 


I.M2T9 

P»4PP2| 

6.669J7 

*.7*54« 

>0 

I.W|« 


6.659S5 

*,7*5*« 


4(tlS97 

PePPtOp 

6.65172 

*.7*54« 



• •HtSB 

•.65266 

*,7*54« 


4.«<eiT 

OePPlPT 

0.656#6 

*,7*5*« 


4.tM$0 


6.65525 

•.7*54* 

tiil 

4.46S4* 


6.65645 

*.7*5*9 


4.m2* 

PePPSSi 

6.657ft| 

*.7*$44 


4.IISS4 

••4«S7p 

6.65679 

*.7*54* 


4.IW4 

t.4P42| 

•.95996 

*.7*549 


4.I4WJ 

e«4P46T 

9.96116 

*.7*54« 


4.1M57 

••44S12 

•.96252 

*.7*5*0 

II.T 

4.17442 


•.9655« 

*,7*54* 


4,|952* 


6.96667 

*.7*54« 


4.III41 


0.96565 

*.7*540 


4,2?74% 

0«PP4»V| 

U.9ft70$ 

*.T*5«v 


4,24*29 

e,467Pi 


4.7459b 

8 T^ 

«,29<b« 

v*^67Vi» 

g.96'»56 

*.T*b4l> 


«.27»9b 

e.r>PI7 

ft.97»Sft 

*.7*5*« 


4.29552 


•.97174 

*.7*54» 


4.549*7 

fePPPS 

6.97292 

*.7*5«« 

^2^ 

4.52941 

• •4H7v 

•.97469 

*.7*54* 

ftT 

4.5425b 

te«7P27 

6.«752f 

*.7*5** 


4.5547- 

t«A7t7$ 

•,•>645 

*.7*5*4 

f. 

4.575*4 

•tP7t2S 

•.97765 

*.7*54* 

»K« 

4.59159 

••47l7| 

•.976ao 

*.7*599 

HI 

4.44775 

0.4722P 

0.97996 

*.7*S4C 

t «2 

4,41*47 

Pe«726P 

6.96116 

•.7*5*» 

hs 

4.441*2 

••47SU 

6.96254 

*.7»59* 

H« 

*,*5*7b 

Pei7SP7 

9.96551 

*.7*54* 

hs 

4.47512 

••P7«IP 

6.96669 

t.7*54« 


*,*•995 

••P7«6% 

6.99567 

*.7*54» 

HV 

4.5»579 

••P751S 

6.96765 

*.7*5*9 

*..• 

4.92214 

P.i7SPS 

•.96622 

*.7*5*# 

Ht 

4,95444 

••<7615 

t,9694» 

*,7*545 

H« 

4.55*44 

•••766ft 

•.99056 

*,7*5*# 

HI 

4.57117 

•*•7717 

6.99176 

*.7*5*9 

H2 

4.9«75l 

••47766 

6.99291 

*.7*54c 

H) 

4.445*7 

•••7616 

6.99911 

*.7*5*« 

H« 

*<4f(>24 

•»47676 

6.99529 

*.7*54* 


4,45b5« 

C*4762l 

6.99647 

*.7*54* 


t.t«p49 

•••7971 

6.99764 

*.7*5*4 

*«2 

*,%»92J 

4,9*559 

•••6425 

6.99662 

4. 7*54# 

H« 

•••••77 

1 .••••• 

4. 7*54* 

H^ 

*.741*2 

^•••iSt 

:.66I1« 

*,7*S*# 


4.71424 

••<6t65 

1.66255 

*.7*5*« 
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TABLE 5-1. Design Data for the U.S. U-25 BIHD Generator (cont) 


9-3341 


FRANC 

AXIAL 

DISTANCE 

WIDTH 

HEIGHT 

diagonal I2AT ION 

ANGLE 

NUN6E(^ 

(a) 

(•) 

(a) 

(radians) 


4.7S*ioi 

U.08237 

1.09353 

0.79549 


4.75<M« 

u. 09293 

1.00071 

0.79599 


4%7a79l 

C. 09307 

t«00599 

9.70599 


4.7®SM 

o.o9oai 

1.00799 

9.70549 

^•9 


0.09055 

1.09929 

0.79549 



«.09St0 

1.09942 

0.79549 

9o7 

4.aSNM 

0.09999 

1.01899 

0.79549 

9®9 

4.®®®»1 

•i.o9621 

1.91177 

9.79599 

ftO 

4.a®597 

0. 09972 

1.01295 

0.79599 

4.aai7s 

0.09733 

l.v|4|S 

9.79549 

Sfl 


4.29799 

I.OI53I 

9.79599 

9|} 

*«• 

tt. 18949 

1.01649 

0.79590 

«i9 

4.®S*s*7€ 

0.19902 

1.91769 

0.79549 

9i« 

«.®47«® 

4.48990 

1.01994 

0.79549 

7t9 

4.49IS9 

0.09017 

1.02992 

0.70599 

4.®7®78 

i. 09075 

1.02129 

0.79549 

9i7 

4.99612 

•♦.19131 

1.02237 

9.79549 

9?9 

S.otAS® 

1.09192 

1.92355 

0.79540 

S9 

$•02891 

“»o925| 

1.02473 

0.79549 

5.e«5|S 

.09311 

1.42591 

0.79599 


5.49109 

c. 09371 

1,02799 

0.79549 


S.‘»7e55 

19033 

1.02931 

0.78193 

?>9 

5.^9629 

0.09902 

1,02959 

0./8OO8 


5.ll39« 

-.49999 

1.03097 

0.77953 


5.13157 

V. 19037 

I.U3213 

0.7T297 


5.10913 

*•09709 

1.03349 

0.76990 

5»7 


V. 19775 

1.03499 

0.76593 

5>* 

5.l9‘«0S 

•09905 

1.03591 

0,:.224 


5.^« »0| 

».19«#15 

I.V3717 

0.7596b 


*>..Mi»T3 

•.19V«S 

l.03«4l 

!>. 15505 



•*5v -So 

I.U3M5 

0.75145 



.'jt 129 

1 ,00o89 

0.7af93 


5.27-31 

-•5v2o0 

1 .v*i2l3 

0.74421 

51* 

5.2^739 

c.5'*27o 

1.00336 

0. 74057 

'i* 

5.3 Ml 

•U50 309 

t .00459 

0.73993 


5.32138 

4.50023 

i *o459| 

0.73329 

^<7 

5.33929 

tt.54097 

i.ooiol 

0,72993 

Si® 

5.355|a 

0.50573 

1.00920 

0.7259? 

Si9 

5.37190 

o.5‘>909 

1.04945 

0.72239 

s«o 

5.39409 

*^.54725 

i.oSeoS 

0.71991 

s«t 

5.U538 

'••b*>90| 

1.05199 

0.71495 

sw 

5.^2202 

1.5* 879 

1 .05309 

0.7112$ 

s«s 

5.13900 

'^•5*»959 

t .05425 

0.70 755 

s«« 

5.15513 

V .5 1035 

1.05544 

0.70394 

'•9 

5.07102 

.51114 

1.05663 

0.70013 

S<9 

5. >0800 

>.51193 

(.05791 

0.99»4I 

S«7 

5.50402 

7.51273 

1,05999 

0.09299 

s«® 

5.52i TO 

.51350 

l,06Of7 

0.9449$ 

S«9 

r 53772 

• .51035 

1.06135 

0.64521 

'i® 

5.55310 

j.5ISt9 

1.09251 

0.9914$ 

79* 

5.59902 

(.*•51599 

1.09399 

0.67779 

7M 

5.58555 

t .51977 

1.09494 

0.67199 

M3 

5.60163 

.51759 

1.09909 

0.67019 


5.01766 

-.51301 

1.0971$ 

0.96639 

'9» 

5.O3360 

1.5*925 

1.06930 

0.96291 

s?® 

5.04957 

M .52049 

1.0694$ 

0.95991 

S®7 

5.O650O 

•U52093 

1.07059 

0.65591 

S5® 

5.69l3t 

.52177 

1.07173 

0.95121 

sn® 

5.09710 

>*•52263 

1.07297 

0.64719 

s%0 

5.71285 

w .52397 

1.07013 

••643$f 
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TABLE 5-1. Design Data for the U.S. U-25 MHD Generator (cont) 


9>3342 


FRME 

AXIAL 

OISTAICE 

mOTH 

HEIGHT 

OIAGOiALIZATIOa 

AMGIC 

lUNBEft 

(a) 

(a) 

(■} 

(radians) 


%»wu 

6.52562 

i.tmt 


f.Taaaa 

6.52960 


6.99591 


5.n«as 

6.52696 

i.mw 

6.99207 

S4« 

S.TTSaa 

0.52979 

t.tms 

6.92629 


a.maa 

0.59119 

l.CWM 

6.62617 

M 

$«aaaa2 

0.59259 

I. Min 

6.92651 


s«aataa 

0.59661 

I.MMI 

6.91990 

ilB 

$e®STaa 

0.59562 

I.M«»T 

6.91277 


s.asaat 

6.59966 

t.M»M 

6.96666 

a.aaaaj 

0.59625 

l.*M>l 

6.90699 

w 

5.6as2a 

6.59965 

l.aM*! 

6 . 90:10 


S.aaaai 

6.56100 

I.OMM 

6.59720 


5.ai5aj 

6.56260 

I.OMM 

6.59220 


$.a?aT9 

0.56965 

I.0.1T2 

6.56916 


S.aaja* 

0.56525 


6.56590 


s.naa? 

0.56909 

i,0Mie 

6.56159 

>r 

Searaoa 

6.56602 

i.eMgi 

6.57100 


a.aaaei 

6.56966 


6.57190 


a.aas9« 

6.59077 

i.gMw 

6.59972 

fno 

a«aiaa7 

0.55215 

I.OMoT 

6.59577 

Ml 

a.aUtip 

0.55952 

I.IMM 

6.59161 

)•? 

a.aaaai 

6.556611 

I.MIM 

6.55765 

Sil^ 

a.eaMi 

0.59025 


6.55167 


a«o7eia 

0.99761 

• iMgg 

6.59990 


a.aaaai 

0.55690 

..loMi 

6.56592 

111* 

a. 19792 

0.59612 

I.IMM 

6.59199 


9.12226 

0.99197 

l.tpYS* 

0 .59799 

»«• 

a.lSaai 

0.999«*2 

I.IM.9 

6.51999 

Sii» 

a.tsssi 

0.59617 

1,1 #MT 

0.52999 


a.iaoo? 

0.59970 

i.iiiao 

0.52591 

1^1 


O.S97*»6 

1.1 12M 

0.52192 


a.iasip 

0.59696 

l.mR3 

0.5176V 


9.20950 

0.99971 

I.IIMI 

6.51160 

s#« 

9.22S9u 

0.97106 

I.IIST* 

6.50969 


9.2S97U 

0.572H 

I.IMM 

0.50569 

m 

9.29272 

6.97996 

I.IIB 12 

0.50179 

»#T 

9.29709 

0.57509 


0.69771 


9.26195 

0.5799? 

I.I2M* 

6.99100 


9.29501 

0.57769 

I.I2IM 

0.60990 


9.90969 

0.57696 

1.12278 

6.96559 

•fti 

9e9290a 

0.56025 

I.12JM 

0.66I67 

•at 

9.99622 

0.56190 

I, 1258* 

6.67719 

•aJ 

9.99297 

0.56260 

1, 12*28 

6.67191 

aa« 

a.99oaa 

0.56619 

i.i2n8 

6.60929 

•a’ 

a.96o59 

0.56560 

I.I28S8 

6.69515 

•a* 

9.99002 

0.56975 

1.I28M 

6.69105 


O.A06O9 

0.56605 

I.II888 

0.65O95 

•aj 

9,62299 

0.56999 

i.ilive 

0.65265 


0.«9902 

0.59092 

I.I13II 

0.66675 

aja 

9.99097 

0.59190 

I.1S825 

6.96909 


a.aaaaa 

0.59916 

I.IS»M 

6.66a52 

a{f 

9.a?09o 

0.59660 

i.lMSi 

6.69660 


9.69222 

0.59576 

I.1S788 

6.61226 

•r 

9^50900 

0.59701 

l.lMTT 

6.62615 

•P 

9.91907 

0.99626 

1,I1M8 

0.62602 


9.59909 

0.59955 

1. 14102 

6.61960 

•1^ 

0.96761 

0.6006? 

1. 18218 

0.61579 

•i! 

9.90110 

0.60206 

I.I8S28 

0.91 too 


9.57609 

0.90S95 

I.I88JI 

0.66799 

ata 

9.96659 

0.60601 

1.1 8M8 

6.90929 
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TABLE 5-1. Design Data for the U.S. U-25 MHD Generator (cont> 


FRAME 

AXIAL 

OISIAICC 

MIDTM 

HEIGHT 

01 AGONAL IZATIOR 
ANGLE 

NUMBER 

(•) 

(•) 

(0) 

(radians) 

•ts 


0.AO9»» 

I.IO00O 

0.10010 



A»A»7ta 

U 10772 

0*10000 

*.*<«*} 

•aAOMlA 

1.10001 

0*10001 


*.••502 

9«A0R6a 

1.10991 

0*10905 


*.*s«9a 

AaAtm 

1.11100 

0*10200 


•.*X#12 

••AtRia 

1.11210 

0*17010 

•P 

a.aas6« 

teA.SST 

1.11120 

0*17012 

a.awia 

••AtAAl 

l.llOlO 

0.10000 


a.Tiaai 

••AI9I*S 

1.15100 

0*10179 


a.T 2 « 0 * 

0«A17v9 

1*11050 

0*10117 


a.TjTaa 

6.AUS2 

1.11701 

0*15710 


a.mai 


1.11072 

0*11119 

•tJ 

a,7a<2a 

6«A29t9 

1,11901 

0*10099 


a.TTTW 

9ea^^•^ 

1.10090 

0*10079 


a.Tatai 

0eA2»f 

1*10100 

0*10o59 


a.aaasj 

a»*2AAe 

1.10107 

0*ll9l9 

•V 

a.*iTM 

d*A2S7a 

1.10010 

0*11210 

•1® 

a.a»c *2 

AaARAAf 

1.10121 

0*12797 

•tt 

a.aaaia 

0*420 IS 

l*10«»ll 

0*12175 

••0 

a.ASTaa 

0*02917 

1*10711 

0*11910 

ut 

a.araaa 

O*0lOS9 

1*10000 
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TABLE 5-2 


0 - 4^'*3 


The U.S. U-25 1-D Magnetic Field Distribution 
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.Axial variation of open and short circuit variables in the U-25 at 50 kg/s 
operating conditions. STD computation BRCYKDXJGJ 
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Fig. 5-3. Axial variation of the plasma properties in the U-25 at 50 kg/s operating 
conditions. STD computation BRCYKDXJGJ 
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Fig. 5-5. Axial variation of the average gasdynamic variables in the U-25 at 50 kg/s 
operating conditions, STD computation BRCYKDXJGJ 



RE8CARCH CPRPORATlON 



2'3 



RtSiAACH COA^OAATIOM 




249 


Fig. 5-7. Axial variation of the channel geometry in the U-25 at 50 kg/s operating 
conditions. STD computation BRCYKDXJGJ 
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Fig. 5-8. Axial variation of the heat tr»nrfer parameters in the U-25 at 50 kg/s 
operating conditions. STD computation BRCYKDXJGJ 
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Pig. 5-9. Ax ar ation of the power extraction parameters in the U-25 at 50 kg/s 

operating conditions. STD computation BRCYKDXJGJ 
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Fig. 5-11. Axial variation of open and short circuit variables in the U-25 at 50 kg/a 
operating conditions. STD computation BRCYEHZJGK 
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Fig. 5-13. Axial variation of the plasma properties in the U-25 at 50 kg/s operating 
conditions. STD computation BRCYEHZJGK 
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Axial variation of electricai field variables in the U-25 at 50 kg/s 
operating conditions. STD computation BRCYEHZJGK 
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Fig. 5-15, Axial variation of the average gasdynamic variables in the U-25 at 50 kg/s 
operating conditions. STD computation BRCYEHZJGK 
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Fig. 5-17. Axial ’'ariation of the channel geometry in the U-25 at 50 kg/s operating 
conditions. STD computation BRCYEHZJGK 
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Fig. 5-18. Axial variation of the heat transfer parameters in the U-25 at 50 kg/s 
operating conditions. STD computation BRCYEHZJGK 
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Fig. 5-20. Axial variation of the optimization variables in the 
operating conditions. STH computation BRCYEHZJGK 
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Fig. 5-21. Axial variation of open and short circuit variables in the U-25 at 50 kg/s 
operating conditions. STD computation BRCYJCGJGK 
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conditions. STD comru at ion BRCVJCGJGK 
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23. Axial variation of the plasma properties in the U-25 at 50 iig/s operating 
conditions. STD oomputaticn BRCYJCGJGK 
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STD computation BRCYJCGJGK 
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Fig. 5-25. Axial variation of the average gasdynumic variables in the 
operating conditions. STD computation BRCYJCGJGK 
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Fig. 5-27. Axial variation of the channel geometry in the U-2f» at 50 kg/s operating 
conditions. STD computation BRCYJCGJGK 


nrSIAflCH CORPORATION 







RESEARCH CORPORATION 



271 


RCbf ARCH CORPORATION 
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Fig. 5-30. Axial variation of the optimization variables in the U-25 at 50 kg/s 
operatin.- conditions. STD computation BRCYJCGJGK 
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Fig. 5-34. Enthalpy extraction as a function of mass flow 
rate for the U-?5 experiment. 
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6.0 


PERFORMANCE OP NASA-SPECIFIED 500 MW(th) MHD GENERATOR 


6.1 Channel Characterization 


The full load operating characteristics of the NASA- 
specified 500 MW(th) MHD channel are given in Table 6-1. The 
generator is defined to be square, i.e. rectangular with aspect 
ratio of 1.0, from inlet to exit. The axial distribution of 
magnetic flux density increases from 4 T at the inlet to 6 T at 
2.5 m from the inlet, then decreases linearly to 4.6 T at the 
exit of the generator. A list of the thermodynamic data to 
define the working fluid for this NASA-specified 500 MWvth) MHD 
channel are given in Table 6-2. The specifications contained 
in Tables 6-1 and 6-2 may not be consistent, i.e., it is 
possible that not all of these specifications can be attained 
simultaneously. 

6.2 Generator Inlet Plane Electrical Conductivity 
Calculations 


The specif Ictlons of the working fluid and inlet 
thermodynamic and fluid mechanical conditions, given In Tables 
6-1 and 6-2, were reduced to species inlet mole ratio, listed 
in Table 6-3, and Inlet thermodynamic conditions. The thermo- 
dynamic conditions are inlet static temperature of 2710 K and 
Inler static pressure of 4.35 atm. 

The inlet transport and thermodynamic properties for the 

NASA-specified 500 MW(th) MHD generator have been computed with 

a zero-dimensional, equilibrium method using the STDNASP code 

from the STD THERMODYNAMICS family of codes. For tf''se 

conditions, the inlet thermodynamic and transport data are 

20 3 

electron concentration of 2.71 x 10 electrons/m and 
electrical conductivity of 6.83 S/m. 
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Even with equilibrium chemistry, the conductivity in the 
inlet plane is more realistically a function of the 
distribution of static temperatures and pressures and also a 
function of the distribution of current. Since the electrical 
conductivity is approximately an exponential function of 
temperature and the electrical current raises the electron 
temperature, both of these two important effects will raise the 
value of average *,lectrical conductivity in the inlet plane 
above that computed with the zero-dimensional, equilibrium 
code. 


For example, at the 0.813 m station (inlet plane) of the 
HPDE the average electrical conductivity from computation 
CHPQUUUJBO is 10.4 S/m, and the electrical conductivity 
calculated for the average static temperature and average 
static pressure is 9.46 S/m. Thus, for this example, the 
multidimensional, nonequilibrium code computes an average inlet 
conductivity which is 10.6% higher than would be calculated by 
a zero-dimensional method. In another example, taken from the 
inlet plane of the Mark VI-C simulation STD68830HI , the average 
conductivity is 12.4 S/m and the zero-dimensional conductivity 
at the average static pressure and average static temperature 
is 12.2 S/m, a difference of only 2.0%. 
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TABLE 6-1 


Pull Load Operating Characteristics 
For a NASA-Specified 500 MW(th) HHD Channel 

1. Channel mass flow rate: 97.5 kg/s 

2. Channel inlet average stagnation pressure: 6. 

3. Channel inlet average stagnation temperature: 

4. Channel inlet height: 0.51 m 

5. Channel inlet width: 0.51 m 

6. Channel length. 14.8 m 

7. Faraday load factor: 0.85 

8. Channel inlet average mach nvimber: 0.85 

9. Burner and nozzle heat loss: 21 MW 

10. Slag heat loss: 3.7 MW 

0-4364 


> atm 
2837 K 
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TABLE 6>2 


Thermodynamic Data to Define the Working Fluid 
For a NASA-Specified 500 MW(th) MHO Channel 


1. Coal and Ash Analysis 


A. Proximate analysis, weight percent 


B. 


Moisture 

22.7 

Volatiles 

29.4 

Fixed Carbon 

39.2 

Ash 

8.7 

Ultimate Analysis, weight 

percent 

Hydrogen 

6.0 

Carbon 

52.1 

Nitrogen 

0.79 

Oxygen 

31.5 

Sulfur 

0.85 

Ash analysis, weight percent 

SIO 2 

37.6 

AlgOg 

17.3 


5.1 

TIO2 

0.7 

P2O5 

0.4 

CaO 

11.0 

MgO 

4.0 

Na 20 

3.1 

K 2 O 

0.5 

SO 3 

17.5 (assumed 


0-4365 
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TABLE 6-2 (Continued) 


D. Heating values 

wet 20.73 MJ/kg 

dry 26.87 HJ/kg 


E. 

Coal dried to 5.0% moisture 


F. 

Coal injected at 300 K 


G . 

Coal carrier gas: none 


Oxidizer 


A. 

Air composition of 300 K, 0.0132 kg 

H 20 /kg dry air 

B. 

Oxygen enriched to 35% Og by volume 


c. 

Preheat to 922 K 


Seed 

A. 

1% K by weight of channel mass flow 

rate 

B. 

0% by weight KgSO^ 


c. 

100% by weight KgCOg 


D. 

anhydrous ^ 200 ^ 


E. 

Seed temperature at 300 K 


F. 

Seed injected downstream of ash removal 

G. 

Seed carrier fluid: none 


Combu 

F '•or 


A. 

Ash removal of 65% 



B. Oxidizer/Fuel stoichiometric ratio of 0.90 

5. Pilot 

A. None 

0-4366 
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TABLE 6-3 


List of Species Used to Define the Working Fluid 
at the Inlet of the NASA-specified 500 MW(th) Generator 


Species 

^5.3302 **4.2192 **0.0693 

SiOg 

AlgOs 

MgO 

TiOg 

CaO 

FegOs 

P 

NagO 

*^2^03 

SO3 


H2O 


Mole Ratio 

^0.0326 %.8716 1-0000 

0.0269 

0.0064 

0.0037 

0.00033 

0.0074 

0.0012 

0.00022 

0.0019 

0.0708 

0.0082 

9.6694 

5.3397 

0.542 


0 - 43‘^7 
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7.0 


SCALE-UP CONSIDERATIONS FOR MHD POWER TRAINS 


In order to put the critical phenomena and the 
performance obtained to dace in the power trains discussed 
above into the context of the overall development program for 
MHD power generation, it is necessary to consider the scale 
dependence of the performance and behavior of MHD power trains. 
In (7-1), STD Research Corporation outlined convenient 
interaction parameters for correlating the performance and 
behavior of high interaction MHD power trains. These 
interaction parameters and the appropriate measures of power 
train performance are defined in Table 7-1. Case E of (7-11 is 
the original "Revision 7-C" desi§n conditions of the U-25 
generator. Case r of (7-1) pertained to the "Revision 4", "hot 
wall" design conditions of the AEDC/HPDE. The scale dependence 
of the performance of these two devices, relative to power 
trains ranging from the smallest (UTSI) to the largest 
commercial scale power train (the "STD Benchmark Power Train"), 
are depicted in Figs. 7-1 through 7-4 and 7-6 through 7-8, 
which are reproduced from (7-1). 

In addition to the quasi-three-dimensional calculations 
discussed in (7-11, STD Research Corporation, under U.S. 
Department of Energy Contract AC-01 -79ET15501 performed a 
series of calculations for cases B, C, F, and G in order to 
elucidate the scale dependence of other critical phenomena 
besides those considered in (7-1). These phenomena included 
quantitative calculations of the scale dependence of secondary 
flows transient behavior, part load operation, and the effects 
of alternate cross-sectional shapes on channel performance. 
While the results of these additional investigations into the 
scale dependence of these phenomena did not generally yield as 
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simple and as elegant correlations as were obtained In [7-1], a 
number o£ interesting and Important relationships were observed 
which will be described in the present section. 

This section reviews the AEDC/HPDE performance to date 
in the context of the scaling parameters presented in Figs. 7-1 
through 7-8 and reviews the additional investigations carried 
out under Contract AC-01-15501. 

7.1 AEDC/HPDE Performance to Date 


Figs. 7-1, 7-2, 7-3, and 7-4 summarize the electrical 
performance obtained by the AEDC/HPDE to date. It is observed 
that, at a magnetic field of 3.27 T, the HPDE has achieved 
approximately half of its original design goals. The enthalpy 
extraction rate is approximately 8% instead of 16%. The 
isentropic efficiency is approximately 37% rather than 
approximately 60%. The specific energy extraction is 
approximately 0.45 MJ/kg versus approximately 0.86 MJ/kg. The 
interaction parameter based on pressure obtained to date is 
approximately 0.99, which is the highest interaction parameter 
based on pressure achieved in modern experiments aimed at 
commercial power generation. This interaction parameter is 
compared to the original design interaction parameter of 
approximately 1.8. 

Fig. 7-5 shows the development of the axial vel^-city 
profiles between the electrode walls and between the sidewalls 
in the Run 006-014 Q3D simulation CHPQUUUJBO. (The first three 
profiles are offset from the main channel outline due to a 
plotter error.) The development of the 15% sidewall overshoots 
at the exit of this generator are consistent with the 
interaction parameter based on velocity, S^, of approximately 
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2 . 2 , which is obtained in this experiment. (As pointed out in 
[7-1 1, supersonic generators have lower interaction parameters 
based on velocity for a given interaction parameter based on 
pressure. The design value of interaction parameter based on 
velocity is in excess of 10, because the original design was 
based on subsonic flow.) The ratio of the sidewall displace- 
ment thickness to the average electrode wall displacement 
thickness is exhibited in Fig. 7-6. It is clear that the 
AEDC/HPDE experiment has already achieved negative displacement 
thickness at the exit of the generator, according to the Q3D 
simulation. If properly instrumented with pressure probes, the 
experiment should be able to observe velocity overshoots under 
the Run 006-014 conditions. Fig. 7-7 shows the dependence of 
total channel blockage at the exit of the generator. It is 
seen that the AEDC generator performance to date has nearly 
passed the milestone of zero blockage at the exit of the 
channel. Finally, the average skin fric *.on coefficient 
throughout the generator is exhibited in Fig. 7-8 as a function 
of the interaction parameter based on velocity. It is seen 
that the sidewall friction, while greater than the electrode 
wall skin friction, has not reached the steep part of the curve 
associated with very strong velocity overshoots on the 
sidewalls. 

It is evident that the AEDC/HPDE offers an early 
opportunity to study multidimensional phenomena in high 
interaction MHD devices. Many of the effects which are 
predicted by the STD/MHD codes will become important as the 
interaction level of the AEDC/HPDE is increased from its 
present levels. 
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7.2 Optimum Mach Number Distribution and Its Effect on 
Par«--Loaci and Transient Behavior 


The Mach number distribution determines performance 
because the power density falls dramatically on either side of 
an optimum Mach number value- The optimum Mach number is 
determined by the thermodynamic and transport properties of the 
working fluid and the stagnation conditions at each generator 
station. A series of part-load and off-design calculations 
performed by STD illustrate the impact of the Mich number 
distribution on generator performance. 

Critical phenomena related to the optimum Mach number 
distribution and its effect >n part-load and transient behavior 
were deduced by analyses and evaluation of computations with 
the TRANSIENT family of codes which were completed under 
Contract AC-01-79ET15501 . A diagram which schematically 
presents the logic of the variations of the part-load and 
unsteady cases which were completed with the TRANSIENT code 
family under is given as Fig. 7-9. For each of four channels 
with thermal inputs of 20 MW, 50 MW, 300 MW, and 2000 MW, a 
nominal steady-state solution was computed. These steady-state 
solutions were the starting point of the following variations; 
(1) 5% step reduction in mass flow, (2) + 5% sine wave input in 
mass flow; (3) 30% quasi-steady mass flow reduction; and (4) 
30% instantaneous load reduction. The electrical connection of 
ejch channel was then redesigned from Faraday to diagonal or 
diagonal to Faraday and the 30% quasi-steady mass flow 
reduction cases were calculated. A list of the final 
computations performed for these cases is presented in Table 
7-2. 
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The Mach number distributions for these channels during 
mass flow rate reduction are displayed in Figs. 7-10 through 
7-17 for both the nominal electrical connections and the 
alternate electrical connections. For the 300 MW Faraday 
connected channel, the mass flow linearly decreases from 100% 
to 70% over the 300 ms time ringe of the computation, as shown 
in Fig. 7-14. For this case, the flow starts out almost 
entirely subsonic. As the mass flow rate decreases, the small 
supersonic region becomes subsonic and the Mach number at each 
station in the entire flow field decreases. 

The Mach number distribution for the 50 MW Faraday 
connected channel is shown in Fig. 7-13. This channel is 
initially completely supersonic with a shock near the end of 
the diffuser. As the mass flow rate decreases to 70% of its 
initial value, the Mach number throughout the channel 
decreases. The Mach number at the channel entrance decreases 
only 0.3% but the Mach number at the channel exit decreases by 
12.5%. 


A correlation of the normalized gross power to 
normalized mass flow rate for the four channels is shown in 
Figs. 7-18 and 7-19. The reduction in power follows very 
closely the reduction in mass flow, except in the case of the 
50 MW(ch) channel. The power output of the 50 MW(th) channel 
remains high as the mass flow decreases because its profile 
does not go far from optimum. 

As the mass flow rate varies, so does the combustor 
pressure and the required compressor power. The combustor 
pressure is given as a function of normalized mass flow rate 
for the 20. 50, 300, and 2000 MW(th) flow trains in Fig. 7-20 
for the nominal electrical hookup and in Fig. 7-21 for the 
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alternate electrical hookup. The net power is given as a 
function of normalized mass flow rate for the 20 and 50 MW(th) 
flow trains in Figs. 7-22 and 7-23. The normalized net power 
is given as a function of the normalized mass flow rate for the 
300 and 2000 MW(th) flow trains in Figs. 7-24 and 7-25. 

The relationships between gross power per unit mass flow 
rate (specific energy extraction) and the product of 
interaction parameter and electrical efficiency are shown in 
Figs. 7-26 and 7-27. For each channel, a trace is shown which 
describes the result of linearly decreasing the ^■'.ass flow from 
its nominal value to 70% of the nominal. For the two subsonic 
channels, nominally 300 HW(th) and 2000 MW(th), both the 
specific energy extraction and the interaction decrease as the 
mass flow rate is reduced. This is attributed to the departure 
of the channel Mach number distribution further from the 
optimum distribution. Conversely, for the two supersonic 
channels, nominally 20 MW(th) and 50 MW(th), t.ie specific 
energy extraction and the interaction increase with decreasing 
mass flow. This effect is attributed to the fact that tie 
Mach number distribution for the supersonic designs of the 
small channels approaches the optimum Mach number distribution 
as the mass flow rate is reduced. These results point to the 
maintenance of the optimum Mach number distribution as an 
important consideration for MHD generator load following. 

The relationship between power output and the specified 
mass flow which varies sinusoidally with a + 5% amplitude, is 
given for the nominal 20 MW(th), Faraday-connected case in Fig. 
7-28 and for the nominal 300 MW(th), Faraday-connected case in 
Fig. 7-29. These figures show both the relative amplitudes and 
the phase difference between the mass fluctuation and the 
resulting power fluctuations. 
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7.3 Effect of Alternate Cross-Sectional Shapes on Channel 

Performance 


The effects of alternate cross-sectional shapes on 
c- '1 performance were deduced by analysis and evaluation of 
cciUf nations completed under Contract AC-01-79ET15501 with 
codes from the Q3DYZ of the Q3D Family of codes. A diagram 
which represents schematically the variations of the cases 
which were completed is presented as Fig. 7-30. '’or each of 
four channels with nominal thermal inputs ol' 20 50 MW, 300 
MW, and 2000 MW, three axial stations were chosen analysis. 
These three axial stations were chosen to be representative of 
the front end, middle, and the back end of the channels. At 
each of these stations and for each channel, three different 
cross-sectional shapes, rectangular, hexagonal, and elliptical, 
were considered. Therefore, there were 36 cases completed. A 
list of these computations is given in Table 7-3. 

The Q3DYZ code has been applied to analyze Faraday, 
diagonal-conc. jcting-wal 1 , and insulating-wall-diagonal 
channels. For the computations listed in Table 7-3, the 20 MW 
and 2000 MW generators are Faraday-connected, the 50 MW(th) 
generator is a diagonal-conducting-wall generator, and the 300 
‘JW generator is an insulating-wall-diagonal. 

At each station in each channel, the load factor, 
enthalpy flux, magnetic field strength, average static 
pressure, wall temperature, anode-center boundary layer 
thickness, cathode-center boundary layer thickness, sidewali- 
center boundary layer thickness, electrode center wall 
temperature, ana sidewall-center wall temperature were the same 
for each different geometry. At each station in the 20, 50 and 
300 MW(th) channels, the cross-section .1 area and aspect ratio 
were also the same for each geometry. 
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The geometrical definitions used for the hexagonal, 
rectangular, and elliptical shapes of these computations are 
given in Fig. 7-31. For all hexagonal shaped channels the 
included electrode angle was 120 degrees. For the 20, 50, and 
300 MW(th) channels, the angles subtendec. ^y the electrodes for 
the elliptical shapes were the same angles, as measured from 
the center point in the channel, as for the corresponding 
rectangular shapes. For the 2000 MW(th) channel, the hexagonal 
cases were regular hexagons and the elliptical cases had aspect 
ratios of unity (i.e., circular). 

The data which have been compiled for the 36 cases are 
listed in Tables 7-4, 7-5, 7-6, and 7-7 for the 2'>, 50, 300, 
and 2000 MW(th) channels, respectively. The data which are 
listed in these Tables are run identifier, shape, axial 
location, magnetic field, load factor, diagonalization angle 
(for the diagonally-connected channel), aspect ratio, cross- 
sectional area, cross-section perimeter, enthalpy flow, 
enthalpy flux, average static pressure, average static 
temperature, core static temperature, average axial velocity, 
core axial velocity, average electrical conductivity, core 
electrical conductivity, average Hall parameter, electrode wall 
temperature, sidewall temperature, anode boundary layer 
thickness, caths^de boundary layer thickness, sidewall boundary 
layer thickness. Hall field, Faraday voltage, load currer. -, 
power output, internal resistance, plasma nonuniformity factor, 
and two power production parameters, Cge? and Cgg. The power 
production parameter Cgg is defined as the power divided by the 
product of load factor, one minus the l»^ad factor, co<'c 
electrical conductivity, magnetic field squared, core a.:ial 
velocity squared, and the cross-sect iona ’ urea, and this 
parameter represents the performance at the station relative to 
an ideal channel with no thermal or velocity gradients. The 
power production parameter Cgg is defined as the power divided 
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by the product of load factor, one minus the load factor, 
average electrical conductivity, magnetic field squared, 
average axial velocity squared, and the cross-sectional area; 
this parameter represents the performance at the station 
relative to a simplified one-dimensional model. 

The analyses and interpretations of the computations 
with alternate cross-sections have resulted in the following 
implications; 

(1) Some of the inefficiency caused by conductivity 
nonuni formity can be removed by appropriate choice of a channel 
cross-section. The larger the nonuniformity factor in a 
particular rectangular channel, the larger the potential gain 
which can be obtained with a nonrectangular geometry. As the 
nonuniformity factor increases with axial distance down the 
channel, the improvement in power obtained from nonrectangular 
geometry should also increase with axial distance. 

In order to obtain the potential increase in 
performance, it is necessary to choose a geometry so that the 
boundary layer voltage drop is minimized, but not a geometry 
which causes the power producing core to become shorted out. 
An incorrect choice of geometry for a particular case can 
result in none of the potential gain from being exploited. In 
most c€.jes, the potential loss from an inappropriate choice of 
geometry is much greater than the potential gain from an 
appropriate choice. 

An example of the reduction of the boundary layer 
voltage drop is given with assistance from Figs. 7-32 through 
7-35. The isopotentials and current streamlines for the 2 m 
axial station of the nominal 50 MW(th) channel are given for 
the rectangular cross-section in Figs. 7-32 and 7-33, 
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respectively. The Isopotentials and current streamlines for 
the corresponding elliptical cross-section are shown in Pigs. 
7-34 and 7-35, respectively. The region of voltage drop near 
the anode and cathode is much larger for the rectangular case. 
Fig. 7-32, than for the elliptical case. Fig. 7-34. The peak 
potential (at the zero field point) for the cathode region in 
the rectangular case is 32 V, and the peak potential for the 
cathode region in the elliptical case is 17 V. That is a 
difference of about a factor of two. 

(2) The geometry (cross-sectional shape) that 
improves the performance of one particular channel can yield 
poorer performance than the rectangular geometry for a 
different channel. 

The optimum geometric shape is a complex function of 
many parameters (boundary layer thicknesses, wall temperatures. 
Hall parameter, to name a few) and the substitution of this 
shape into a different channel with different values for these 
important parameters will yield different results. An 
illustration of this can be found with the hexagonal shapes for 
the 20 MW(th) and 300 MW(th) channels. For each of the 
stations for the 20 MW(th) channel, the hexagonal case had 
better performance than the corresponding rectangular case, but 
for each of the stations for the 300 MW(th) channel, the 
hexagonal case had worse performance than the corresponding 
rectangular case. 

(3) As shown by the results for the 50 MW(th) channel 
cases, compared to the results from the three other channels, 
the electrical performance of a completely diagonal-conducting- 
wall channel is not as strong a function of cross-sectional 
shape as for the insulating wall Faraday-connected channel. 
Since the diagonal-conducting-wall channel has all walls as 
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conductors, the current can enter the walls where it is most 
efficient for the particular geometry and distributions of 
plasma properties. However, this does not say that the DCW 
channel is better than a Faraday-connected channel. 

(4) Geometries which have less surface-to-area ratios 
and larger hydraulic diameters than the rectangular cross- 
sectional shape, but have the same area, have smaller losses 
due to heat transfer to the walls, wall shear stresses, and 
viscous effects. For particular cases, each of these three 
losses can be simultaneously decreased by over 10%. 

(5) Geometries which allow the current to enter the 
electrodes parallel to the magnetic field have a lower internal 
plasma impedance. The lower internal plasma impedance can be 
utilized to improve the performance of the channel. 

An example of the effect of lower plasma impedance is 
given with assistance from Figs. 7-36 through 7-39. The 
isopotentials and current streamlines for the 1 m station of 
the nominal 20 MW(th) channel are given for the rectangular 
cross-section case in Figs. 7-36 and 7-37, respectively. The 
isopotentials and current streamlines for the corresponding 
hexagonal cross-sections are shown in Figs. 7-38 and 7-39, 

. espectively. As shown on Pig. 7-37, the current enters the 
electrode nearly perpendicular to the magnetic field for the 
rectangular case and, as shown on Fig. 7-39, the current enters 
the electrode nearly parallel to the magnetic field for the 
hexagonal case. The internal resistance of the hexagonal case, 
illustrated by Figs. 7-38 and 7-39, is 32% less than that for 
the rectangular case. At the same Faraday load factor of 0.5 
(for maximum power production), the Faiaday voltage is 13% 
lower for the hexagon, but the load current is 28% higher, so 
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that the power output is 12% higher for the hexagon than for 
the rectangle. 

Although the cases represented in Table 7-3 and Fig. 
7-30 define a systematic approach to analysis of the effects of 
cross-sectional geometry, the detailed analyses performed to 
date have not resulted in complete and unambiguous 
interpretations which are applicable in all cases. These 
calculations of performance with variable cross-sectional 
shapes did not include sufficient parametric variation and were 
applied over too great a range of conditions to obtain regular 
trends with scale and position in the channel. It is probable 
that many of the elliptical or hexagonal cross-sections with 
lower performance than the rectangular design values could be 
optimized to outperform the rectangle when additional 
systematic parametric variations are considered. In addition, 
other cross-sectional shapes (e.g. , octagonal) may provide 
better performance than any of the three shapes considered. 
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TABLE 7-1 

Interaction Parameters and Performance Indices for 
Evaluation of MHD Generator Performauce Scaling 


where 


Interaction parameter based on pressure 


I 

P 


•o P 


dx 


Interaction parameter based 
u pU 


on velocity 


Generalized nonuniformity factor 
G * g + (g-KP^)) 
g = (»•>(</«»■> 


Electrical conversion efficiency 

T » E 


r • (u X B) 


Isentropic efficiency 
IP 

^ ^^^isentropic 


Enthalpy extraction ratio 

IP power out 

mHp thermal power in 


Specific energy extraction 


4- : power per unit mass flow rate 
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0-4354 
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List of Final TRANSIENT Computations Performed Under Contract AC-01-79ET15501 

( continued) 
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Thermal Electrical 
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TABLE 7-3 


List of Final Q3DYZ Computations Performed Under 
Contract AC-01-79ET15501 


Run 


Date 

Nominal 

Thermal 

Flux 

(MW) 

Cross- 

Sectional 

Shape 

Axial 

Location 

(m) 

BRCYCHPIFG 

17 

June 

1979 

20 

rectangle 

0.50 

BRCYCJAIFQ 

17 

June 

1979 

20 

rectangle 

1.00 

BRCYCJCIFQ 

17 

June 

1979 

20 

rectangle 

1.49 

BRCYCISIFQ 

17 

June 

1979 

20 

hexagon 

0.50 

BRCYCIUIFQ 

17 

June 

1979 

20 

hexagon 

1.00 

BRCYCJKIFQ 

17 

June 

1979 

20 

hexagon 

1.49 

BRCYCIQIFQ 

17 

June 

1979 

20 

ellipse 

0.50 

BRCYCIWIFQ 

17 

June 

1979 

20 

ellipse 

1.00 

BRCYCJFIFQ 

17 

June 

1979 

20 

ellipse 

1.49 

BRCYBIEIFQ 

17 

June 

1979 

50 

rectangle 

1.00 

BRCYBIIIFQ 

17 

June 

1979 

50 

rectangle 

2.00 

BRCYBIQIFQ 

17 

June 

1979 

50 

rectangle 

3.00 

BRCYBHFIFQ 

17 

June 

1979 

50 

hexagon 

1.00 

BRCYBILIFQ 

17 

June 

1979 

50 

hexagon 

2.00 

BRCYBIHIFQ 

17 

June 

1979 

50 

hexagon 

3.00 

BRCYBHQIFQ 

17 

June 

1979 

50 

ellipse 

1.00 

BRCYBRNIFQ 

17 

June 

1979 

50 

ellipse 

2.00 

BRCYBIQIFQ 

17 

June 

1979 

50 

ellipse 

3.00 

BRCYBGLIFQ 

17 

June 

1979 

300 

rectangle 

2.00 

BRCYBGBIFQ 

17 

June 

1979 

300 

rectangle 

3.50 

BRCYBFZIFQ 

17 

June 

1979 

300 

rectangle 

5.00 

BRCYBGHIFQ 

17 

June 

1979 

300 

hexagon 

2.00 

BRCYBGFIFQ 

17 

June 

1979 

300 

hexagon 

3.50 

BRCYBFLIFQ 

17 

June 

1979 

300 

hexagon 

3.00 

BRCYBGJIFQ 

17 

June 

1979 

300 

ellipse 

2.00 

BRCYBGDIFQ 

17 

June 

1979 

300 

ellipse 

3.50 

BRCYBFXIFQ 

17 

June 

1979 

300 

ellipse 

5.00 

CHRQBUEIF3 

29 

June 

1979 

2000 

rectangle 

4.00 

CHRQLNIIF3 

29 

June 

1979 

2000 

rectangle 

7.50 

CHRQBVFIF3 

29 

June 

1979 

2000 

rectangle 

11.00 

CHRQEHRIF3 

29 

June 

1979 

2000 

hexagon 

4.00 

BRCYEHPIF3 

29 

June 

1979 

2000 

hexagon 

7.50 

CHRQBXJIF3 

29 

June 

1979 

2000 

hexagon 

xl.OO 

CHRQDWBIF2 

28 

June 

1979 

2000 

ellipse 

4.00 

BRCYHQJIF2 

28 

June 

1979 

2000 

ellipse 

7.50 

BRCYINQIF2 

0-4368 

28 

June 

1979 

2000 

el li pse 

11.00 
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Fig. 7-1. Variation of enthalpy extraction rate with 

the interaction parameter based on pressure. 
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research corporation 



Fig. 7-2. Variation of Isentroplc effxciency with in- 
teract Mn parameter based on pressure. 
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SPECIFIC ENERGY EXTRACTION, P/m (Mj/kg) 



INTERACTION PARAMETER TIMES LOAD FACTOR, IpK 


Fig. 7-3. Variation of specific energy extraction with 

the product of interaction parameter based 

on pressure and the Faraday load factor K 

or diagonal generator electrical conversion 

efficiency n • 
e 
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AVERAGE CONDUCTIVITY NONUNIFORMITY 
FACTOR INCREMENT, (g-1) 



Fig. 7-4 . Variat ion of conductivity nonuniformity 

factor g, averaged over the channel length, 
with sensible heat rate or thermal input tc 
the MHD generator. Solid symbols denote 
arc-mode current transport , and triangles 
denote slagging walls. 
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VITY NONUNIFORMITY FACTOR 




Fig. 7-5. Development of the axial velocity profiles In the Q3D 
simulation CHPQUUUJBO of the AEDC/HPDE P.un 006-014. 
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Fig. 7-7. Variation of the MHD generator exit blockage 
with interaction parameter based on velocity. 
The STD computational Cases A-G are compared 
with the estimates used by the ETF conceptual 
design contractors and the AEDC performance 
to date. 
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Fig. 7-8. Average skin friction coefficient on the 

electrode walls and sidewalls as a function 
of the interaction parameter based on 
velocity . 
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Fiji. 7-9. Schematic duurram representing the cases which were completed with the 
TRANSIENT code under contract AC-01 -79ET15501 . 
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NORMALIZED GROSS POWER 


M ■ Ur/UA’' 
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NORMALIZED MASS FLOW RATE 


CHANNEL MASS FLOW POWER 
kg/sec MW 


20 MWth 2. AS 0.A23 

50 MWth 8.00 1.819 

300 MWth 50.0 39.91 

0-4114 2000 MWth 705.1 602.9 


Fig. 7-18 Normalized gross power vs. normalized mass flow- 
rate for channels with their nominal electrical 
hook-up . 
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NORMALIZED MASS FLOW RATE 


CHANNEL MASS FLOW POWER 

kg/sec 

20 MWth 2.49 0.502 

50 MWth 8.00 1.777 

300 MWth 50.0 36.67 

0-4111 2000 MWth 705.9 522.9 

Fl". 7-19 Normalized gross power vs. normalized mass flow rate 

for 4 channels with their alternate electrical hook-up. 
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COMBUSTOR PRESSURE (atm) 



0-4108 NORMALIZED MASS FLOW RATE 


Fig. 7-20. Combustor stagnati^J.*. pressure vs. normalized mass 
flow rate for different size generators with mass 
flow varying between 70% and 100%> at nominal design 
value. 
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Fig. 7-22 . Net power vs. normalized mass flow rate for two 
generators with mass flow varying between t 

100% of nominal design value. 
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Fig. 7-26 . Gross power/mass flow rate vs. interaction param- 
eter times electrical efficiency for four different 
size channels with mass flow varying from design 
to 70% of design. 
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Fig. 7-27 . Gross power/inass flow rate vs. Interaction param- 
eter times electrical efficiency for four different 
size channels with mass flow varying from design 
to 70% of design. 
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Fig, 7-28. Normalized gross power and normalized mass flow 
rate for an inputted ±5% sine wave in mass flow 
for a 20 MWth channel with Faraday connection. 
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Fig. 7-29. Normalized gross power and normalized mass flow 
rate for an inputted + 5% sine wave in mass flow 
for a 300 MWth channel with Faraday connection. 
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8.0 


SUMMARY AND CONCLUSIONS 


This section summarizes the major conclusions derived 
from the studies described in this report. 

8. 1 AEDC/HPDE Data Analysis and Interpretation 

Attention has been concentrated on the data from two 
power producing runs of the AEDC/HPDE (Runs 006-008 and 
006-014). Calculations have been performed with codes from the 
Q3D and TRANSIENT families in order to compare theory and 
experiment where possible. 

The comparison of theory and experiment for Run 
006-008 is incomplete because the experimental data were 
inadequate to accurately define the conditions that existed 
during the operation of that run. Nevertheless, the 
qualitative agreement is sufficiently good to presume that if 
the Hall voltage across the channel were known, the discharge 
current ccuid have been accurately determined and the 
comparison would have been excellent. 

The comparison between theory and experiment was much 
better for Run 006-014. The data from the two recording 
systems were much more consistent for this run and are more 
reliable. On the whole, the predictions are well within the 
scatter for the data. 

8.2 Critical Phenomena in the AEDC/HPDE 

The power dissipation in the thermal boundary layers 
of the HPDE is much more important than was originally assumed 
in the design of the HPDE. However, the enthalpy extraction 
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goal of the experiment can be met without providing the 
original design values of electrode wall temperatures (2100 K - 
1950 K). On the basis of the calculations performed to-date, 
provision of a uniform 1500 K wall temperature would be 
sufficient to meet the enthalpy extraction goal of the 
experiment . 

Under the operating conditions of Run 006-014, which 
represents the best-diagnosed, high-power run of the HPDE to 
date, the principal operational parameters available to 
increase the HPDE performance toward the goals of the 
experiment are the electrode surface trmperature and magnetic 
field. 


The central role of the Mach number distribution in 
achieving the performance goals of the experiment is 
demonstrated by the calculations performed to date. The 
optimum Mach number distribution in the AEDC/HPDE appears tc be 
transonic . 


Finite electrode and near electrode effects appear to 
have a four to five percent impact on the performance of the 
MHD generator under the operating conditions of Run 006-014. 
The experimental data therefore are perhaps not sufficiently 
refined to discriminate between the differences calculated with 
some models for near electrode effects (small arcs vs. diffuse 
mode transport), and electrode segmentation. The data to date 
do appear to rule out the existence of very strong arcs and 
suggest that the severe losses associated with seed 
condensation cn the elect. >des and seed depletion in the 
boundary layers are not controlling mechanisms. 
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The AEDC/HPDE computations performed to date iudicate 
that the fluid mechanical behavior of the generator is 
approaching the interesting regime of zero blockage and 
measurable sidewall velocity oversnoots. The data thus far arr 
not sufficiently complete to detect these effects. Howeve-’, 
the general agreement with the pressure predictions and the 
other electrical and gasdynamic predictions with the 
multidimensional Q3D calculations suggests that the inclusion 
of these effects in the computations may be necessary. As the 
performance of the MHD generator improves, and the interaction 
level increases, the differences between simplified ana 
sophisticated models of the flow field will become more and 
more apparent. The AEDC/HPDE should be afforded additional 
instrumentation with which to diagnose the local fluid 
phenomena which are predicted by the analyses. 

The unconventional fluid behavior will be magni.f’.ed 
when the generator is operated in the subsonic mode. It would 
be possible, under conditions similar to HPDE Run 006-014, to 
achieve an interaction parameter based on velocity sufficient 
to exhibit many of the secondary flow and velocity overshoot 
effects predicted by the STD/MHD codes, if the mass flow rate 
were reduced sufficiently to provide subsonic flow throughout 
the generator. While this may not be the optimum enthalpy 
extraction and isentropic efficiency point, it would be 
possible, neverthess, to simulate conditions in a large, 
subsonic, baseload MHD generator, if the appropriate 
interaction parameter based on velocity is provided. 

The inlet pressure loss was calculated fo one 
operating condition to be approximately 1.5% of the total 
pressure entering the generator. Such a pressure loss in the 
inlet and exit regions of the generator would not present a 
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substantial penalty to the overall performance of an MHD plant 
operating with an MHD generator of the type represented by the 
AEDC/HPDE. The exit region pressure loss remains to be 
evaluated. 

The electrode voltage drops are sensitive not only to 
the wall temperature but also to the amount of Joule 
dissipation heat deposited in the boundary layers. Second 
order improvement in performance may be achieved by tailoring 
the load factor distribut on to achieve high dissipation in the 
inlet regions of the generator. Analytical models which fail 
to account for Joul*; dissipation effects in the electrode 
boundary layers will fail to predict the electrode voltage 
drops with accuracy. 

8 . 3 Critical Phenomena in the U.S. U-25 Experiment 

The nominal conditions of the U.S. U-25 are predicted lo 
yield transonic behavior with a shock located in the channel. 
Even under these nonideal conditions, the U.S. U-25 experiment 
is predicted to De able to reach or come very close tc the goal 
of 10 MW power output. The operation of this experiment can be 
altered by changing the mass flow rate, but simultaneous 
attainment of both fully subsonic flow and at least 10 MW power 
output was not predicted under the nominal operating 
conditions. 

8.4 Electrical Conductivity Calculations i.. Coal Flames 

The calculation of electrical conductivity of the 
products of coal combustion is a subject of considerable 
uncertainty due to the lack of sufficiently accurate 
compositional and thermochemical data. Even for equM i 
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calculations, the elemental composition of the products of 
combustion is dependent upon history of the combustion, 
particularly f.r the treatment of ash and slag. In addition, 
for MHD generators the aonequi 1 ibrium effects of boundary 
layers, current transport, radiation transport, condensation, 
and chemical kinetics must be taken into account. 

Even when tne coal crmposition is defined (defined for 
computational purposes), the values of the transport and 
thermodynamic iroperties at the inlet to the generators is 
dependent upon the conditions in the combustor, separator, 
throat, nozzle, and mating sections of the flow train. It is 
particularly important to know the distributions of heat flux 
and slag removal in these sections of the flow train in order 
to define the inlet condiiions. 

For examp’ ' .ere are many ways of modeling slag 
removal: a) all of the liquid and solid slag species can be 
removed at the flame temperature, (b, the slag that is removed 
can be assumed t- have the same composition as for the ash 
analysis; (c) the condensed phase slag can be removed at the 
equiliorium composition of the gas-phase, condensed-phase 
interface temperature; or (d) any combination of these, i.e., 
some slag removed at one condition and some removed at another. 
Since condensed phase matter is removed to the walls 
continuously along each section of the flow train, the model of 
slag .emovai must be applied as a function of axial distance 
and for conditions which are functions of location in the flow 
train . 
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